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Abstract: The generalized abc-block edge transformation graph Q*°~ (G) is a graph whose 
vertex set is the union of the edges and blocks of G, in which two vertices are adjacent 
whenever corresponding edges of G are adjacent or one corresponds to an edge and other to 
a block of G are not incident with each other. In this paper, we study the girth, covering 
invariants and the domination number of Qt°~(G). We present necessary and sufficient 
conditions for Qt" (G) to be planar, outerplanar, minimally nonouterplanar and maximal 
outerplanar. Further, we establish a necessary and sufficient condition for the generalized 


abc-block edge transformation graph Q*°- (G) have crossing number one. 


Key Words: Line graph, abc-block edge transformation, generalized abc-block edge trans- 
formation graph, Smarandachely block-edge H-graph. 
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§1. Introduction 


Throughout the paper, we only consider simple graphs without isolated vertices. Definitions 
not given here may be found in [5]. A cut vertex of a connected graph is the one whose removal 
increases the number of components. A nonseparable graph is connected, nontrivial and has 
no cut vertices. A block of a graph is a maximal nonseparable subgraph. Let G = (V, E) be a 
graph with block set U(G)={B;; B; is a block of G}. If a block B € U(G) with the edge set 
{e1, €2,°+* ,€m;m > 1}, then we say that the edge e; and block B are incident with each other, 
where 1 <i<m. The girth of a graph G, denoted by g(G) , is the length of the shortest cycle 
if any in G. Let [x](|a]) denote the least (greatest) integer greater (less) than or equal to z. 
A vertex and an edge are said to cover each other if they are incident. A set of vertices in 
a graph G is a vertex covering set, which covers all the edges of G. The vertex covering number 
ao(G) of G is the minimum number of vertices in a vertex covering set of G. A set of edges 
in a graph G is an edge covering set, which covers all vertices of G. The edge covering number 
~ 1§upported by UGC-SAP DRS-III, New Delhi, India for 2016-2021: F.510/3/DRS-III/2016(SAP-I) Dated: 


29!” Feb. 2016. 
2Received December 23, 2017, Accepted May 8, 2018. 
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a1(G) of G is the minimum number of edges in an edge covering set of G. A set of vertices in 
a graph G is independent if no two of them are adjacent. The maximum number of vertices in 
such a set is called the vertex independence number of G and is denoted by {9(G). The set of 
edges in a graph G is independent if no two of them are adjacent. The maximum number of 
edges in such a set is called the edge independence number of G and is denoted by (1(G). 

The line graph L(G) of a graph G is the graph with vertex set as the edge set of G and two 
vertices of L(G) are adjacent whenever the corresponding edges in G have a vertex in common 
[5]. The plick graph P(G) of a graph G is the graph whose set of vertices is the union of the set 
of edges and blocks of G and in which two vertices are adjacent if and only if the corresponding 
edges of G are adjacent or one is corresponds to an edge and other is corresponds to a block are 
incident [8]. In [2], we generalized the concept of plick graph and were termed as generalized 
abc-block edge transformation graphs Q*°(G) of a graph G and obtained 64 kinds of graphs. 
In this paper, we consider one among those 64 graph which is defined as follows: 


Definition 1.1 The generalized abc-block edge transformation graph Q*°- (G) is a graph whose 
vertex set is the union of the edges and blocks of G, in which two vertices are adjacent whenever 
corresponding edges of G are adjacent or one corresponds to an edge and other to a block of G 


are not incident with each other. 


Generally, a Smarandachely block-edge H-graph is such a graph with vertex set E(G) ) B(G) 
and two vertices e1,e2 € E(G)U B(G) are adjacent if e1,e2 € E(H) are adjacent, or at least 
one of €), €2 not in H(H) and they are non-adjacent, or one in E(H) and other in B(G) which 
are not incident, where H is a subgraph of G with property YW. Clearly, a Smarandachely 
block-edge E(G) J) B(G)-graph is nothing else but a generalized abc-block edge transformation 
graph. 

In this paper, we study the girth, covering invariants and the domination number of 
Q*°-(G). We present necessary and sufficient conditions for Q*t°~(G) to be planar, outer- 
planar, minimally nonouterplanar and maximal outerplanar. Further, we establish a necessary 
and sufficient condition for the generalized abc-block edge transformation graph Q*° (G) have 
crossing number one. Some other graph valued functions were studied in [3, 4, 7, 8, 9, 11, 12]. 
In Figure 1, a graph G and its generalized abc-block edge transformation graph Qt°~ (G) are 


shown. 


Figure 1. Graph G and its Qt°- (G). 
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In Qt°~(G), the vertices correspond to edges of G denoted by circles and vertices corre- 
spond to blocks of G denoted by squares. The vertex e, (B,) of Q+°- (G) corresponding to edge 
e; (block B;) of G and is refereed as edge (block)-vertex. 


The following theorems will be useful in the proof of our results. 


Theorem 1.1([8]) Jf G is a nontrivial connected (p,q) graph whose vertices have degree d; and 
P 


if b; the number of blocks to which vertex v; belongs in G, then P(G) has q—p+1+ > b; 
i=1 


12 
vertices and 5 a edges. 
— 


Theorem 1.2([5]) For any nontrivial connected graph G with p vertices, 


ao(G) + Bo(G) =p = ar(G) + f1(G). 


Theorem 1.3((6]) Jf L(G) is the line graph of a nontrivial connected graph G with q edges, 
then 


§2. Basic Results on Q*°- (G) 


We start with preliminary remarks. 


Remark 2.1 L(G) is an induced subgraph of Qt°~(G). 
Remark 2.2 If G is a block, then Q*°- (G) = L(G) U Ky. 


Remark 2.3 Let G be a graph with edge set E(G) = {e1,€2,--- ,@m} and r blocks. Then 
dg+o- (aye =dgej+r—-1. 


Remark 2.4 Let G be a (p,q)-graph with block set U(G) = {Bi, Bo,--- ,B,} such that 
| E(B;) |= Nj. Then dg+o-(@)B; = qd- Gy. 


Theorem 2.1 Let G be a (p,q)—connected graph whose vertices have degree d; with r > 1 
blocks and b; (1 <i <p) the number of blocks to which vertex v;, belongs in G. Then 


p 
(1) The order of Qt°-(G) =q-—pt1+4+ >> bi; 
i=l 
p 
(2) The size of QT°-(G) = g(r — 2) + 3 OG. 
i=l 


Proof It is shown in [5] that for a connected graph G with p vertices and b; number of 
blocks to which vertex v; (1 < i < p) belongs in G. Then the number of blocks of G is given 


P 
by D(G) = 1+ 32 (b; — 1). The order of Qt°~(G) is the sum of the number of edges of G and 
i=1 
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number of blocks of G. Hence the order of Q*°- (G) 


Pp 
qt1+ 5 -(i-1)=q-p+1+)_ bi. 


i=l i=1 


The total number of edges formed by joining each of the r block-vertices to all the q 


iw 
edge-vertices is rg. The number of edges in line graph L(G) is —q + 5 S- d?. Thus, the size of 
i=l 


Qh’-(@) =rq-a-ats LY df = a(r 2) + he. 


t=1 


An immediate consequence of the above theorem is the following corollary. 


Corollary 2.2 Let G be a (p,q) graph whose vertices have degree d; with r blocks and m 
components. If b; (1 <i <p) is the number of blocks to which vertex v; belongs in G, then 


P 
(1) The order of Q*°-(G) =q-—p+m+t yO 


i=l 


P 
(2) The size of Qt°-(G) = q(r —2)+ +5> d?. 
i=l 


Theorem 2.3 Let G be a graph. The graphs Qt°~(G) and P(G) are isomorphic if and only if 
G has two blocks. 


Proof Let G be a (p,q) graph with r > 1 blocks. Suppose Q*°-(G) = P(G). Then 
|E(Q*°- (G))| = |E(P(G))|. By Theorems 1.1 and 2.1, we have 


Pp 
q(r — 2) be = se 


aes = 0. 


Since G has at least one edge and hence equality holds only when r = 2. Therefore G has 
two blocks. 


Conversely, suppose G has two blocks B, and By. Then by definitions of Q*°-(G) and 
P(G), L(G) is induced subgraph of Qt°~(G) and P(G). In Q*°-(G), block-vertex Bj is 
adjacent all the edge-vertices corresponding to edges of Bz and block-vertex Bs is adjacent to all 
the edge-vertices corresponding to edges of By. In P(G), block-vertex Bj is adjacent all the edge- 
vertices corresponding to edges of B, and block-vertex Bs is adjacent to all the edge-vertices 
corresponding to edges of By. This implies that there exist a one-to-one correspondence between 
vertices of Qt°~ (G) and P(G) which preserves adjacency. Therefore the graphs Qt°~ (G) and 
P(G) are isomorphic. 


The following theorem gives the girth of Qt° (G). 
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Theorem 2.4 For a graph G 4 2K, Ko, P3, 


37f G contains K1,3 or K3 or G=P,; n>4 or Gis union of at least 
two cycles or paths or G is union of paths and cycles, 
9(Q*°-(G)) =) 4if G=mK2,m>4, 
6if G=3ka, 
nif G=Cy,n>=> 4. 


Proof If G contains a triangle or K1,3, then the line graph L(G) of G contains triangle. 
By Remark 2.1, it follows that girth of Qt°~(G) is 3. Assume that G is triangle free and K1.3 
free. Then we have the following cases: 


Case 1. Assume G has every vertex of degree is 2. We have two subcases: 


Subcase 1.1 If G is connected, then clearly G = C, ;n > 4, we have Qt9-(G) = C, UK, 
n> 4. Therefore girth of Q*°" (G) is n. 


Subcase 1.2 If G is disconnected, then G is union of at least two cycles and Q*°- (G) 


contains at least two wheels. Therefore girth of Qt°~ (G) is 3. 
Case 2. Assume that G 4 2K2, K2 has every vertex of degree is one. It is easy to see that 


6 if G= 3Ka, 
4if G=mKo; m> 4. 


9(Q*°"(G))= 


Case 3. Assume that G # P3 has vertices of degree one or two. Then G is either union of 


paths P,, or union of paths and cycles. Therefore girth of Q*°~ (G) is 3. 


§3. Covering Invariants of Q*° (G) 


Theorem 3.1 For a connected (p,q)-graph G with r blocks, if Q*°-(G) is connected, then 
ao(Q*" (G)) = and 8o(Q* (G)) =r. 


Proof Let G be a connected (p, q)-graph. By Remark 2.1, L(G) is an induced subgraph of 
Qt°-(G). Therefore by definition of Q*°~(G), the edge-vertices covers all the edges of L(G). 
Since Q*°- (G) is connected, it follows that for each block-vertex B’ of Qt°~ (G), there exists 
a edge-vertex e’ such that e’ and B are adjacent in Q*°- (G). Therefore the vertex set of L(G) 
covers all the edges of Q+°- (G) and this is minimum covering. Hence ao(Qt°- (G)) = q. Since 
Qt°-(G) is connected. By Theorem 1.2, we have ap(Q*° (G)) + Bo(Q*° (G)) =q+r. Thus 
Bo (Q*°"(G)) = r. 


Theorem 3.2 Let G be a connected (p,q)-graph with r blocks. If Qt°~ (G) is connected, then 


r if Gis a tree, 


a1 (Q*°"(G))= 


r+[5*] otherwise. 
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and 


Bi(Qt?- (G))= q if Gis a tree, 


q—|4 |] otherwise. 


Proof Let T be the set of minimum edges covering all block-vertices of Q*°-(G). ice., 
|T| = r. Let S be the set of minimum edge cover of L(G). By Theorem 1.3, |S| = [$]. We 


consider the following two cases: 


Case 1. If G is a tree, then g = r. By the definition, T covers all block-vertices and edge- 
vertices of Qt°~(G). Thus a1(Q*°" (G)) = r. 


Case 2. If Gis not a tree, then g > r. By the definition, T covers all block-vertices and only r 
edge-vertices of Q+°- (G). Therefore there exists a set of edge-vertices F', say of Qt°- (G) such 
that no element of T is incident with any element of F in Q*°~(G). ie., |F| = q—r. Since 
each element of S covers two elements of L(G) and F c V(Q*°- (G)), it follows that we need 


[S| elements from S' to cover all elements of F. Thus a1(Qt°- (G))= r + [5*]. 


2 


r if G is a tree, 
Therefore, a1(Qt° (G))= 


r+[5*] otherwise. 
Since Qt°~(G) is connected. By Theorem 1.2, we have a1(Q*°- (G)) + 61(Q*° (G)) = 


qtr. Thus 


q if G is a tree, 


(QT (G))= 


q—[4*] otherwise. 


§4. Domination Number of Qt°-(G) 


A set D of vertices in a graph G = (V,£) is called a dominating set of G if every vertex in 
V — D is adjacent to some vertex in D. A dominating set D is called minimal dominating set 
if no proper subset of D is a dominating set. The domination number 7(G) of a graph G is the 


minimum cardinality of a dominating set in G ((10]). 


The following result is immediate from Remark 2.2. 


Theorem 4.1 If G is a block, then y(Q*°- (G)) = y(L(G)) +1. 
Theorem 4.2 If G has two blocks, then y(Qt°~(G)) = 2. 


Proof Suppose G has two blocks By and By. Then Bj dominates all the edge-vertices 
in Q*°-(G) corresponding to edges of Bz and Bs dominates all the edge-vertices in Qt°~ (G) 
corresponding to edges of B,. Therefore y(Q*°- (G))= |{ Bj, BS}| = 2 where {B, BS} is a 
minimal dominating set in Qt°~ (G). 


Theorem 4.3 For any graph G with at least three blocks, 
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2 if G contain an edge is adjacent to every other edge of its block, 
VQ (G))= 


3 otherwise. 


Proof Let G be a graph having at least three blocks. We consider following two cases: 


Case 1. If G contain an edge e is adjacent to every other edge of its block B, then block- 
vertex B’ dominates the edge-vertices corresponding to the edges not in B. And edge-vertex e’ 
dominates the block-vertices except B’ and dominates the edge-vertices corresponding to edges 
of B. Therefore y(QT°-(G)) = |{e’, B’}| = 2 where {e’, B’} is a minimal dominating set in 
gt-(a). 


Case 2. If G contain no edge is adjacent to every other edge of its block, then there exist two 
block-vertices B’, Bi and one edge-vertex e’, where e is in B in G, such that B’ dominates the 
edge-vertices corresponding to the edges not in B and edge-vertex e’ dominates all the block- 
vertices except B’ and block vertex Bi dominates the edge-vertices which are not dominated 
from e’ and B’. Therefore y(Qt°~(G))= |{e’, B’, B{}| = 3 where {e’, B’, Bi} is a minimal 
dominating set in Q*°-(G). 


§5. Planarity of Graphs Qt? (G) 


A graph is planar if it can be drawn on the plane in such a way that no two of its edges intersect. 
A planar graph is outerplanar if it can be embedded in the plane so that all its vertices lie on 
the exterior region. In [1], Kulli introduced the concept of a minimally nonouterplanar graph. 
The inner vertex number i(G) of a planar graph G' is the minimum possible number of vertices 
not belonging to the boundary of the exterior region in any embedding of G in the plane. 
Obviously G is outerplanar if and only if i(G) = 0. A graph G is minimally nonouterplanar if 
i(G) = 1. An outerplanar graph G is maximal outerplanar if no edge can be added without 
losing outerplanarity. The crossing number Cr(G) of a graph G is the minimum number of 
pairwise intersections of its edges when G is drawn in the plane. Obviously, Cr(G) = 0 if and 
only if G is planar. A cactus is a connected graph in which every block is an edge or a cycle. 
If G and H are graphs with the property that the identification of any vertex of G with an 
arbitrary vertex of H results in a unique graph, then we write G- H for this graph. 

The condition for the planar, outerplanar, minimally nonouterplanar, maximal outerplanar 
and crossing number of line graph of G and generalized abc-block edge transformation graph 
Q*°-(G) are same when G is a block. So that in this section we assume graph G under 


consideration is not a block in what follows. 


Lemma 5.1 Jf G is not a tree having more than two blocks, then Qt°~(G) is nonplanar. 


Proof Let G be not a tree having more than two blocks, i.e., G has a block B contains a cycle 
C. Then Q*°-(G) has a subgraph homeomorphic to Qt°~ (2K2 U Ks), and Cr(Qt9- (2K2 U 
K3)) =1. Therefore Q*° (G) is nonplanar. 


Theorem 5.2 Let G be a connected graph with more than one block. Then generalized abc- 
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block edge transformation graph Qt°~ (G) is planar if and only if G satisfies one of the following 


conditions: 


(1) G is a cactus having two blocks; 
(2) G is a tree of order < 5. 


Proof Suppose Qt°~ (G) is planar. Assume a connected graph G has atleast 5 blocks. We 


consider the following cases: 
Case 1. If G is not a tree, then by Lemma 5.1, Qt°~ (G) is nonplanar, a contradiction. 


Case 2. If Gisa tree, i.e., every block of G is K2, then Qt°~ (G) has a subgraph homeomorphic 
to Qt°- (5K) and Cr(Qt?- (5K2)) = 4. Therefore Qt°~ (G) is nonplanar, a contradiction. 

In either case we arrive at a contradiction. Hence G contains at most four blocks. We 
discuss two possibilities on number of blocks: 


Subcase 2.1 If G is not a cactus having two blocks, i.e., some block B of G contains a 
subgraph homeomorphic to C;, + e, then edge-vertices corresponding to edges of C,, + e and 
block-vertex corresponding to block other than B forms a subgraph with at least one crossing 
in Q*°-(G). Therefore Q*°- (G) is nonplanar, a contradiction. This proves (1). 


Subcase 2.2 If G is not a tree having 3 or 4 blocks, then by Lemma 5.1, Qt°~(G) is 
nonplanar, a contradiction. This proves (2). 

Conversely, suppose G satisfies (1) or (2). Then G=C,, - Kg or C,- Cy or Py or K1,3 or 
K,3- Kz or P3 or Ps. Therefore it is easy to check that Q*°~(G) is planar. 


Theorem 5.3 Let G be a connected graph with more than one block. Then generalized abc-block 


edge transformation graph Q*°-(G) is outerplanar if and only if G is a tree of order < 4. 


Proof Suppose Qt°~(G) is outerplanar. Then Q+°~(G) is planar. By Theorem 5.2, we 
have, G is a cactus having two blocks or G is a tree of order < 5. Assume G is a tree of 
order 5. Then Qt°- (G) has a subgraph homeomorphic to Q*°- (4K2) and i(Qt°- (4K2)) = 4. 
Therefore Q*°- (G) is nonouterplanar, a contradiction. Assume G = Cyn, - Cy or Cp» Ko. Then 
Q*°- (G) is nonouterplanar, a contradiction. In either case we arrive at a contradiction. Hence 
G is a tree of order < 4. 

Assume G is not a tree of order < 4, i.e., G has a block B contains a cycle C. Then 
edge-vertices corresponding to edges of C and a block-vertex corresponding to block other than 
B forms a subgraph wheel in Q*°~ (G). Therefore Q*°~ (G) is nonouterplanar, a contradiction. 
Hence G is a tree of order < 4. 

Conversely, suppose G is a tree of order < 4. Then G = P3 or Py or Ky 3. Therefore 
Q*°-(G) is outerplanar. 


Theorem 5.4 Let G be a connected graph with more than one block. Then generalized abc-block 
edge transformation graph Q*°- (G) is minimally nonouterplanar if and only if G = C,,- Ke. 


Proof Suppose Q*°-(G) is minimally nonouterplanar. Then Q*°~(G) is planar. By 
Theorem 5.2, we have, G is either cactus having two blocks or tree of order < 5. If G is a tree 
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of order < 4, then by Theorem 5.3, Q*°~ (G) is outerplanar, a contradiction. If G is a tree of 
order 5, then Qt°~ (G) has a subgraph homeomorphic to Q*° (4K), and i(Qt°~ (4Ka)) = 4. 
Therefore Qt°~ (G) is not minimally outerplanar, a contradiction. 

Suppose G 4 C,, - Kg is cactus having two blocks. Then G = P3 or Cy, - Cy. Therefore 
Q*°-(G) is not minimally nonouterplanar, a contradiction. Thus G is C, - Ko. 


Conversely, suppose G = C,,- Ka. Then Q*°- (G) is minimally nonouterplanar. 


Theorem 5.5 Let G be a connected graph with more than one block. Then generalized abc-block 
edge transformation graph Qt°~(G) is maximal outerplanar if and only if G = Ki 3. 


Proof Suppose Q*°~(G) is maximal outerplanar. Then Qt°~(G) is outerplanar. By 
Theorem 5.3, we have, G is a tree of order < 4. Assume G # Kj 33 is a tree of order < 4. 
Then G = P3 or Py. Therefore Q*°- (G) is not maximal outerplanar, a contradiction. Hence 
G= Ky3. 


Conversely, suppose G = K1,3. Then Qt°~(G) is maximal outerplanar. 


§6. Graphs Q*°- (G) and Crossing Number One 


Lemma 6.1 Let G be a connected graph having two blocks. Then generalized abc-block edge 
transformation graph Qt°-(G) has crossing number one if and only if G is either Cy + (Cs +e) 
with A(G) < 4 or Ko: (Cy +e). 


Proof Suppose Qt°~(G) has crossing number one. Assume G # C;:(C,+e) with A(G) < 4 
or K2-(C, +e). Then we have the following cases: 


Case 1. If G is a cactus, then by Theorem 5.2, Q*°- (G) is planar, a contradiction. 


Case 2. If G is not a cactus, then G is homeomorphic to Ky - (C; + 2e) or Ke» (K2U K3) 
or (C, + e)- (Cs +e) or Cy: (C, +e) with A(G) = 5. Therefore Cr(K2 - (Cy + 2e)) > 2, 
Cr(K2:(K2U K3)) > 2, Cr((Cy+e):(C.+e)) > 2, Cr(C;-(Cs+e)) = 2. Hence Cr(Q*®- (G)) > 
2, a contradiction. 

Conversely, suppose G is either C;- (C, + e) with A(G) < 4 or Ko: (C,; +e). Then 
Cr(Qt?-(@) = 1. 


Theorem 6.2 Let G be a connected graph with more than one block. Then generalized abc-block 
edge transformation graph Qt°~(G) has crossing number one if and only if G = C;- (Cs, +e) 
with A(G) < 4 or Ko: (Cs +e) or C,,- P3. 


Proof Suppose Qt+°~(G) has crossing number one. Assume G #4 C;-(C,+e) with A(G) < 4 
or K2-(C, +e) or C;, - P3. We consider the following cases: 


Case 1. If G is a tree, then we consider following subcases: 


Subcase 1.1 If G is a tree of order < 5, then by Theorem 5.2, Qt°~(G) is planar, a 


contradiction. 
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Subcase 1.2 If G isa tree of order at least 6, then Qt°~ (G) has a subgraph homoemorphic 
to Qt°- (5K) and Cr(Qt°- (5K2)) =4. Therefore Cr(Q*°- (G)) > 4, a contradiction. 


Case 2. If G is not a tree, then G contains at least one cycle. We consider the following 
subcases: 


Subcase 2.1 If G has more than 3 blocks, then Qt°~ (G) has a subgraph homeomorphic to 
Qt°- (3K2U K3) and Cr(Qt?- (3K2U K3)) =5. Therefore Cr(Q*° (G)) > 5, a contradiction. 


Subcase 2.2 If G has three blocks, then Qt°~(G) has a subgraph homeomorphic to 
Q*°- ((C4+e) -P3) or Gy where G; = Kj —e, e is pendant edge, and Cr(Qt°- ((Cy+e) - P3)) > 
4, Cr(Qt°-(G)) = 2. Therefore Cr(Qt°~ (G)) > 2, a contradiction. 


Subcase 2.3 If G has two blocks, then by Lemma 6.1, crossing number of Q*°- (G) is 
not equal to one, a contradiction. 

Conversely, suppose G=C; - (C; + e) with A(G) < 4 or Ky-(C, +e) or C, - P3. Then 
Qt°-(G) has crossing number one. 
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Abstract: In this study, we investigate the classical differential geometry of isotropic 
curves in the complex space C*. We examine the constant breadth of isotropic curves and 
obtain some results regarding these isotropic curves. We express some characterizations 
of these curves via the E. Cartan derivative formula. We also indicate that the isotropic 
vector of these curves and pseudo curvature satisfy a third order vector differential equation 
with variable coefficients. We study this differential equation in some special cases. We 
dene evolute and involute of the isotropic curve and express some characterizations of these 
curves in terms of E. Cartan equations. The isotropic rectifying curve and isotropic helix 
are characterized in C*. Finally, we present the conditions for an isotropic curve to be an 


isotropic helix. 


Key Words: Complex spaces, isotropic helix, isotropic curve of constant breath, Bertrand 


curves, iotropic rectifying curves. 
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§1. Introduction 


At the beginning of the nineteenth century, V. Pancelet’s isotropic curve opened a door for a 
number of new concepts. The imaginary curve in the complex space was pioneered by Cartan. 
He defined his moving frame and the Cartan equations in C3. Altimisik extended the Cartan 
apparatus of isotropic curves to C+. Furthermore, isotopic Bertrand curves and isotropic helices 
in C? were characterized, [9], [10], [16]. Also, the concept of a slant helix in the complex space 
in C* was offered by Yilmaz [13]. 

Curves of constant breadth were introduced by Euler [3]. The curves have been studied 
in different spaces by researchers. For instance, Izumiya and Takeuchi defined slant helices [5]. 
Ali and Lopez gave some characterizations of slant helices in Minkowski 3-space [1]. Yilmaz 
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studied spherical indicatrices of curves in Euclidean 4-space and Lorenzian 4-space [14], [15]. In 
[7], Magden and Yilmaz extended the well known properties of constant breadth of the curves 
in four dimensional Galilean space G*. 

Many researchers have studied involute-evolute curves in other spaces. The Frenet appa- 
13 and E* is given [4], [8]. In [12], Turgut 
and Yilmaz studied involute-evolute curve couple in Minkowksi space-time. Semin mentioned 


ratus of involute-evolute curves couple in the space 


involute-evolute isotropic curve in [11]. In Euclidean 4-space, rectifying curves are introduced 
by Ilarslan and Nesovig in [6] as space curves whose position vector always lies in its rectifying 


plane, spanned by tangent, the first binormal and second binormal vector fields 7, B, and 
74. 


Bz.The position vector of a rectifying curve a in according to chosen origin satisfies the 


equation 


a(s) = A(s)T(s) + y(s)Bi(s) + u(s)Ba(s), 


where A, y and ys are some differentiable functions of the pseudo arc-length parameter s. 

Thus, the main goal of this paper is to define some isotropic curves in the four dimensional 
complex space C*. In the present paper, we first study isotropic curves of constant breadth 
and the involute-evolute of the curve in C*. Then we introduce the Bertrand curve and present 
some characterizations of the mentioned curves in terms of E. Cartan equations. Also, we give 
a new characterization of the isotropic helix. Throughout this study some complex curves are 
characterized in the complex space C*. 


§2. Preliminaries 


To meet the requirements in the next sections, the basic elements of the theory of imaginary 
curves in the space C4 are briefly presented (a more complete elementary operation can be 
found in [11]). 


Let xp be a complex analytic function of a complex variable t. Then the vector function 


is called an imaginary curve, where t = t; + ite, x : C — C4 and k, are standard basis unit 


vectors of E*, i? = —1. An arbitrary vector x € C‘, is called an isotropic vector if and only if 
, y 


x? = 0, (x £0). In this space, the curves for which the square of the distance between any two 
points equal to zero, are called minimal or isotropic curves [11]. Let s denote pseudo arc-length 
(for details, see [10] or [11]). Then, a curve is an isotropic curve if and only if 


ds? = dx? = 0. 


The complex four dimensional space C*, is the real vector space E* endowed with the 
standard flat Euclidean metric given by 


g = dx? + 2daxidzr3 — dzi, 
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where (21, 22,23, 24) is the complex coordinate system of C’. 


The E. Cartan frame moving along the isotropic curve x in the space C4 is denoted by 
{e1,e2,e3,e4}. This frame is defined ([11]) as 


e) = x 

eo = ix” 

e = — Ex +x” 

es, = ple; A eo Aes3) (2.1) 


where 3 = (x’”)?, wis taken as +1. If y is taken as +1, the determinant of matrix [e1, e2, e3, e4], 
the E. Cartan frame becomes positively oriented. Here, the triple vector product is cross product 
expressed as in [2]. The inner products of these frame vectors are given by 


0 ift+j =1,2,3mod 4 
e,:e; = 1 ift+j=4 
—-1 if*+j7=8 


where the vectors e; and e3 are isotropic vectors; eg is real and e4 is a complex vector. E. 
Cartan derivative formulas can be expressed as follows: 


e, = —1e2 
e, = tke, + 1€3 
e, = —ikes (2.2) 
é 
e, = —€(k" + €k)e, — Eke3 4 rae 


where k(s) = $/3(s) is the pseudo curvature of the isotropic curve in the class C° and €(s) = 

4 1 _ (iv) 2 7 : J 7 

+——>——.,, where 7(s) = (x , the derivative being taken with respect to the pseudo arc 
TGTETE ya) e™) g p P 


length s. In the rest of the paper, we shall suppose pseudo curvature is non-vanishing except 
in the case of an isotropic cubic. 


An isotropic hypersphere with centre m and radius r > 0 in C* is defined as 
S° = {p = (pi, p2, p3, pa) € C* : (p— m)* =r}. 


Definition 2.1 An isotropic curve x = 2(s) in C* is called an isotropic cubic if the pseudo 


curvature k(s) =0, where s is the pseudo arc-length parameter of the curve. 


Definition 2.2 Let x= 2(s) be a complex curve in C+. If the pseudo curvature of the curve is 


constant, then x(s) is called a pseudo helix or isotropic helix in C*. 


Definition 2.3. An isotropic curve «= «(s) in C* is called an isotropic helix if inner product 
Pp 


of its tangent vector e, is constant with some fixed isotropic vector v, that is, e,.v = constant. 
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Definition 2.4 Let x= 2(s) be an isotropic curve in C4. If there exists another isotropic curve 
x* = x*(s) in C* such that principal normal vector field x* coincides with that normal vector 
field of x, then «x is called a Bertrand curve and x* is called the Bertrand mate of x and vice 


versa, where x(s) and x*(s) are opposite points of the curve. 


Definition 2.5 Let p and 6 be two unit speed complex curves in C*. If the tangent vector of 
the curve yp at the point y(so) ts orthogonal to the tangent vector of the curve 6 at the 6(so) 


then curve 6 is called the involute of the curve y as follows: 
g(€1p; €15) = 0, 


where {€iy, €2y, 63, Cap} and {e€15, €25, €35, eas} are Frenet frames of p and 6, respectively. 
Also, the curve y is called the evolute of the curve 6. This definition suffices to define this 
curve mate as 6 = pt requ. 


Definition 2.6 Let a be a complex curve in C+. A rectifying curve is defined in C* as an a 
isotropic curve whose position vector always lies in orthogonal complement ex of its principal 


normal vector field eg. 


§3. Isotropic Curves of Constant Breadth and Their Characterizations 


Let x(s) and x*(s) be isotropic curves in C*. If the tangent isotropic vector e1 of x(s) coincides 
with the tangent isotropic vector ej of x*(s) opposite directions at the corresponding points 
and the distance between these points is always constant, then x(s) is a constant breadth of the 
isotropic curve. Suppose that x(s) and x*(s) are isotropic curves of constant breadth. Then ej 
can be expressed by 


e, = —e} 


where e; and ej are inverse direction and parallel vectors. 
Let x(s) and x*(s) be isotropic curves of constant breadth in C*. Taking into account the 
Cartan equations, it can be decomposed by 


X*(s) = X(s) + mi(s)e1 + mo(s)eg + m3(s)e3 + m4(s)e4, (0 < s < 1), (3.1) 


where X(s) and X*(s) are opposite points and e1, e2,e3,e4 denote the E. Cartan frame in C4. 


Differentiating the equation (3.1) with respect to s, we get 


dX= dX” ds* _ ot ds* 
ds ds* ds ‘ds 
d d 
— eae + moaik + mam )e1 t ( mii ae ms3tk)es 
ds ds 
_ dm dm 
+(moi + ee + manz)e3 + (—= + mans )ea, (3.2) 


where m(s) = —€(k” + €k), n2(s) = —€k,3(s) = £ and k = 4/3 is a pseudo curvature of the 
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isotropic curve in the class C°. Since ej = —e1, we obtain 
pea aa cs 
dg ee ds 
dm 
—myti+ i m3ik = 0 
ds 
dm 
moi + —= — mano =0 
ds 
dma 
— =0. 3.3 
ae als (3.3) 


Puttin, s) =—-1-— &, in the equation (3.3), it can be written as 
g ds 


dm : 
tt = —maik — mam + f(s) 
d 
Ge = mii + m3ik 
ds 
dmg r 
=— —m 
as Mgt 412 
dmsz 
—_ =- : 3.4 
We Mans ( ) 


By virtue of the equation (3.4), (i.e. the fourth expression of the equation (3.4)) we have 
ma = c is constant. Rearranging the equation (3.4) we get 


BOS ope Sei ee Bs) 
ds 
a 
us = mit msik 
ds 
GEN a rr op spe (3.5) 
ds 


The following corollary is a consequence of the equations (3.4) and (3.5). 


Corollary 3.1 Let x = a(s) be an isotropic cubic. The isotropic position vector of x with 


respect to E. Cartan frame can be formed by the equations (3.5) and can be obtained as 
x(s) = x"(s)+ (/ f(s)ds + ku(s)) ey 


te (| / (/ F(s)ds) + ka (s)ds " in(s)) e: 


i (fof (asda) = kx(s) ds + tke(s) + kx(s)) es +.ce4. 


Proof Let x = x(s) be an isotropic cubic. Then, & = 0 from Definition 2.1. From equation 
(3.5); we get 44 = f(s). Integrating this expression we have, 


m= f f(s)ds+ hs, 
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where k, is a complex constant, from equations (3.4),(3.5)2 and (3.5)s3, 
ee (ff Or ee TOD ea 
fig = fof f(s)ds) + ku(s)5)ds) + tke(s) + h(s) 


and m4 = c is constant. After mj ,m2,m3 and mz, are substituted into the isotropic position 


vector x = x(s), the proof is completed. 


Theorem 3.1 Let x= 2x(s) be complex curve of constant breadth with pseudo arc-length in C+. 


If x= x(s) lies fully in the e3e4 subspace, then «= x(s) is an isotropic helix. 


Proof Let x = x(s) be the pseudo arc-length parameter of constant breadth of complex 
curve in C+. From equations (3.5), if we take m, = m3 = 0, then we have mz = c, (where c; is 


a constant). Using this expression in the third equation of (3.5), we obtain k = 47 is constant. 


From Definition 2.3), it is clear that the curve x = x(s) is an isotropic helix. 


Theorem 3.2 Let x= 2x(s) be complex curve of constant breadth with pseudo arc-length in C+. 


There is no constant breadth of isotropic curve that lies fully in the e,e2 subspace. 


Proof Let x = x(s) be the pseudo arc-length parameter of constant breadth of complex 
curve in C+. If we take m3 = m4 = 0 in equation (3.5), we get m; = 0 and mz = cki. So 


x = x(s) lies fully in the eye2 subspace. 


Theorem 3.3 Let x = a(s) be complex curve of constant breadth with pseudo arc-length in 
C*. There is no constant breadth of complex curve which lies fully in the e,e4 subspace, and 


x = a(s) is isotropic cubic. 


Proof Let x = x(s) be the pseudo arc-length parameter of constant breadth of complex 
curve in C*. From equation (3.5), we get m1 = 0, m4 =c and k= 0. So x = x(s) lies fully in 


the e;e, subspace. From Definition 2.1, x = x(s) is an isotropic cubic. 


Theorem 3.4 A pseudo arc-length isotropic x= x(s) in C* is of constant breadth if and only 


if it satisfies the following third order differential equation. 


Proof From equation (3.5)1, we get 


a ck" + c&k — f(s) + dia 
—ik ' 


m2 


Substituting into (3.5)2, this expression m3 is obtained 


dm 
d ck’ +c&k— f(s)+— ; 
ds “ak Sete 
m3 = 


ik 
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Taking the derivative of this expression, we obtain 


d_ ( ck" +e&k~f(s) +73 
dmz ad | % ( =k ) a 
ds ds k 


Substituting into equation (3.5)3, this expression, we have a differential equation of third 
order with complex variable coefficients as follows: 


d| 1d ck" + c&k — f(s) + Z d es 
ds ik ds —itk ds\ k 
-z (cx” + c&k — f (s) 4 oe + ck =0. (3.6) 


The differential equation of third order with variable coefficients in equation (3.6) is char- 
acterized for the constant breadth of isotropic curve x = x(s). 


Now, we characterize the distance between opposite points of the curves of constant breadth 
in C+. Remember the equation (3.1) 


X*(s) = X(s) + mi(s)e1 + ma(s)eo + m3(s)es 
+ma(s)e4,(0 <5 <1). 


If the distance between opposite points of (C) and (C*) is constant, then we can write that 


\|a* — x|| = m2? + 2mym3 — mz = P = constant. (3.7) 
Hence, we write 
dmg dm, dmg dma “< 
ue ds ins ds Tare ds See ds aa oe) 


from equations (3.5) since m4 = c is constant. Rearranging the equation (3.8), we obtain 


dmg dm, dmg a 
ae TB te Ea =0. (3.9) 


Considering equations (3.5), we have 


/ 
k 
m3 |u(s) — ki - =| =0. (3.10) 
m3 
We write m3 = 0 or p(s) — k?i — mack = 0, obviously, m3 4 0. Then it can be expressed 


in the following cases: 


Case 1. Let us suppose m3 = c, 4 0 constant. From equations (3.5)2 and (3.5)3 we easily 


have mz = cki,m, = —c,k. Then the isotropic position vector of y* can be written as follows: 


yp =p+ cake, + ckieg + cye3 + ces. 
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Case 2. Let us suppose that m3 is constant and ¢ is isotropic helix. Thus, the equation (3.6) 
takes the form F 
FOS) kh(s) + ck® =0, (3.11) 
ds? 
where h(s) = c&k — f(s). The solution of the equation (3.11) is 


1 
h(s) — LieV** + Loe VF + 3 = Woe (3.12) 


where L, and Lz are real numbers. 


Case 3. Let us suppose 
mck _ 0 


p(s) — k7i— (3.13) 


m3 


In this case, (C*) is transformed by the constant vector 7 = m 1e1 + m2e2 + m3e3 + m4e4 
of (C). Now, let us investigate the solution to Case 3. 


Suppose that y is an isotropic cubic. Then, we get from equation (3.13) j(s) = 0 and from 


equation (3.5) we get m1 =constant, mz = 0,m3 = —¢. 


§4. Involute and Evolute of Isotropic Curves in C* 


Theorem 4.1 Let y and 6 be complex curves and y be an evolute of 6. The Cartan apparatus of 


Pl Eig, C29; C30; €4y, ky} can be formed according to the Cartan apparatus of 6{€1, €2, €3, e4, k}. 


Proof Let ~ and 6 be complex curves and ¢ be an evolute of 6. According to the property 


of involute-evolute curve couples, we have 
p=d+ re. (4.1) 


Differentiating both sides of the equation (4.1) with respect to s, we obtain 


dp dsy dX : 
—_— .— = — (= : 4.2 
ds, ds oe ds\! ata) ee) 
Rearranging equation (4.2), we have 
dy ds, dX 
ar re (44 )e, — : 4. 
dsy ds (et a Mes (4:3) 


Similarly, based on the definition of involute and evolute curves, we can say ej,1e. 
Obviously, we get 
1+ Ne 0 (4.4) 
ag ; 


We get \ = c— 8, where c is constant. Rearranging the equation (4.1), we get 


p=6+(c—s)er. (4.5) 
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By differentiating the equation (4.5), we have the following equation 


ds 


19 


ee eye = (c— s)(—ieg). (4.6) 
Taking the norm of both sides, we get 
1p = —€2 (4.7) 
and d 
sek (c—s)i (4.8) 
ds 
Differentiating the equation (4.6) two times with respect to s, we get 
yp" = —(c— s)ke; + iep + (c— s)eg (4.9) 
and 
yp” = [-2k + (c— s)k’] e; + [—2i(c — s)k] eg — 2e3. (4.10) 
Thus, we have the following expressions for e2,,e3, and ky. 
C25 = (c coe s)kie, —eg + (c na s)ies 
e35 = [-2k + (c— s)k']e1 + i(c Ne 2k)eo — 2es (4.11) 
and 
ky = —2[-2k + (c— s)k'] + [-2(e— s)k]*. (4.12) 
Using the exterior product a(e1y A e2y A e3y), we get 
C4y = 0 [2(c — s)ik(1 + (c— 8)k)ea], (4.13) 
where o = +1. 
Since from equation (4.7), it follows that e,, is not an isotropic vector, we can state the 
following. 


Remark 4.1 Let y be an evolute of a complex curve in C*. The curve y cannot be an isotropic 


curve. 


Theorem 4.2 Let y and 6 be complex curve and vy be an evolute of 6 in C+. The evolute y 


cannot be an isotropic helix in C4. 


Proof Considering the definition of isotropic helix, we write 


€1,-V = constant, 


(4.14) 
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where v is a constant isotropic vector. From equation (4.7), we easily have 
—e2.v = constant, (4.15) 
Differentiating both sides of equation (4.15), we get 
—(ike; + ie3).v = 0. (4.16) 
Therefore vLe; and v_Le3. Let us decompose v as 
v=tyeo + toes. (4.17) 


Differentiating equation (4.17) consecutively and using Cartan equations, we have t; = 0 
and tg = 0. According to the result, we write 


v=0. (4.18) 


Equations (4.14) and (4.18) yield a contradiction. Therefore, evolute y cannot be an 


isotropic helix in space C7‘. 


§5. Bertrand Couple Curves of Isotropic Curves in C4 


Theorem 5.1 Let a* and a be Bertrand curves in complex space C*. The Cartan apparatus 
of a*{et, e, e3,e4,k*} can be formed by the Cartan apparatus of a{e1, €2, €3, e4, k}. 


Proof Suppose that {a(s),a*(s*)} is an isotropic Bertrand pair of curves. Then a*(s*) 
can be expressed by 
a*(s*) = a(s) + A(s)eo, (5.1) 


where \(s) is the non zero analytic function and s* is the pseudo arc-length parameter of a*(s*). 
Differentiating both sides of the equation (5.1) with respect to s, we get 


_ da®* ds* 
~ ds* ds 


* 


ds* dX 

= ee = (14+ Aki) ey + —eg + Mes. (5.2) 
ds ds 

The definition of Bertrand curves yields ef Leg. Multiplying both sides of equation (5.2) 


with eg we have 
dv _ 


ae 
which implies that A is constant. Using this in the equation (5.2) and taking the norm of the 


0 (5.3) 


both sides, we get 
d. 2k 
= 2(1 + Aki) Ni 
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and the tangent vector e7 is equal to 


1 ka } 
e; = =a ED 2 a oe (5.4) 
2 (1+ Aki) i 2 (1+ Aki) i 


Taking the derivative of the equation (5.2) two times with respect to s, we get 


a®” = (1+ Ak’i)e; + (—1 — Aki + Ak)eo (5.5) 


at” = (14 AKG — ki + AK? + ARP i)ey + (—1 — AWE + AR’ Jen + (-i-— AK + Aki)es. (5.6) 
Using the equation (5.6), we get the vectors e}, e} and pseudo curvature k*, as follows: 


1 


1 
= mal 1— Aki + 2(1 + ARG — ki + AK74)(—4 — AK + Aki)Jer 
+(—1— dk’ + Ak’ Jeo + [-1 — Aki + Ak]e3} 


and 


k* = 5{ 1— Aki + AR + 2(1 + ARG — ki + AK? + AKPZ)(—i — Ak + Aki)}. 


So, the pseudo curvature k*(s) is a non zero constant. 


Remark 5.1 Obviously, e7] isn’t an isotropic vector from equation (5.4). So, the Bertrand 


curve a* cannot be an isotropic curve. 


Remark 5.2 Let a* and a be Bertrand curves in C+. If one of the Bertrand curves lies fully 
in e,e2e3 subspace, then the Bertrand mate also lies fully in the same subspace of C+. 


Theorem 5.2 Let x = x(s) be an isotropic curve in C*. Then, x(s) is a pseudo isotropic helix 
if and only if the following statements are equivalent: 


(a) det (x! (s),x!"(s), x) (s)) =0; 


(b) det(e; (s), e7(s), e'(s)) = 
c) det(e(s), €3(s), €3"(s)) = 0. 


Proof Taking the 2”4, 3"¢ and 4*” derivatives of equation (2.1), we obtain 


x! =e, x" = —1e2, x!” = ke, + e3, x(t) = k'e, a ike. (5.7) 
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We calculate that 


0 -i 0 
det (x"(s),x!"(s),x(s))=| k 0 1 |=—ztk. 
ki -2ik 0 


Since x = x(s) is an isotropic helix, then det(x’(s),x!”(s),x)(s)) = 0. Conversely, let the 
statement ”a)” be true. Then det(x"(s),x!”(s),x”)(s)) = —ik’ = 0. Thus, k is constant and 
x(s) is an isotropic helix. This completes the proof ”a)”. Similarly, denoting x’ = e,,x” = 
ex!” = ef and x) = e/’, we easily see that ”a)” and ”b)” are equivalent. Also, because of 


the fact that the equations 


e; = —ikeg 
e, = ke, —ik'eg + kes 
ej’ = 3kk’e; — (2k? + ik” eo + 2k’es, 


are hold, we can calculate that 


det(e',(s), ef (s), ef/(s)) = — 8K! =0. 


Since x = x(s) is an isotropic helix, then det(e3(s), e3(s),e4’(s)) = 0. Conversely, let us 
3 3 3 


say that in the determinant, 


det(e}(s),e%(s),e4/(s)) =| k? —ik’ k | =—kk' =0, (5.8) 
3kki — (2k? + ik”) 2k! 


we get ak = 0 then k is a constant. As an immediate consequence of Definition 2.2, x = x(s) 


is an isotropic helix. 


§6. Isotropic rectifying curves in C* 


In this section, we firstly characterize the rectifying curves in C* in terms of their pseudo 
curvature. In analogy with Euclidean four dimensional case, we define the rectifying curves in 
complex space C* as a curve whose position vector always lies in the orthogonal complement 
et of its principal normal vector field eg. Hence, e+ is a three dimensional subspace of C7‘, 
spanned by vector filed e;,e3 and e,. Therefore the position vector with respect to some chosen 


origin of a rectifying curve a in C4, satisfies the equation 
a(s) = A(s)ei(s) + wes(s) + 6(s)ea(s) (6.1) 


for differentiable functions X(s), w(s) and 6(s) with pseudo arc-length parameter s. Firstly, let 
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us characterize the rectifying curve a in C+ in terms of its pseudo curvature. Let a = a(s) be a 
unit speed complex rectifying curve in C*, with non zero pseudo curvature k(s). By definition, 
the position vector of complex curve a satisfies equation (6.1) for some differentiable functions 
A(s), u(s) and 6(s). Differentiating the equation (6.1) and using Cartan derivative formulas 
(2.2), we get 
: ; a 
[V — 1— 6&(k” + Ek)]ex + [Ai — pikjeg + [A’ — 6&k]e3 + [6 + oz les =0. 


It follows that 


N'(s) — 4(8)&(s)(k"(s) + &(s)k(s)) = 1 
A(s)i — w(s)k(s)i = 0 
V(s) — (s)&(s)A(s) = 0 
/ ols)! = 
0’(s) + 6(s) gs) 0 (6.2) 
and thus 
Ns) = ef k(s)ds 
ls) = gf beyas 
d(s) = Zs)" (6.3) 


Conversely, assuming that the pseudo curvature k(s) of an arbitrary unit speed complex 


curve a in C4, satisfied the following equation 


we (ef k(s)as) east (<5 [ He)as) oes (=) ets 


Remark 6.1 (i) a@ cannot be an isotropic cubic, since Ey #0; 


(it) If @ is a helix, then a(s) = s [(ck Jer + (c)e3 + (<5) e| ; 
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§1. Introduction 


The notion of quasi Einstein manifold was introduced in a paper [8] by M.C.Chaki and R.K.Maity. 
According to them a non-flat Riemannian manifold (M",g),(n > 3) is defined to be a quasi 
Einstein manifold if its Ricci tensor S of type (0,2) satisfies the condition 


S(X,Y) =ag(X,Y) + bA(X)A(Y) (1) 
and is not identically zero, where a,b are scalars b £ 0 and A is a non-zero 1-form such that 


g(X,&1) = A(X), VX € TM, (2) 


where, €; is a unit vector field. 

In such a case a, b are called the associated scalars. A is called the associated 1-form. Such 
an n-dimensional manifold is denoted by the symbol(QE)n. 

Again, U.C.De and G.C.Ghosh defined generalized quasi Einstein manifold. A non-flat 
Riemannian manifold is called a generalized quasi Einstein manifold if its Ricci-tensor S of 


type (0,2) is non-zero and satisfies the condition 
S(X,Y) =ag(X,Y) + bA(X)A(Y) + cB(X)B(Y), (3) 
where a, b,c are non-zero scalars and A, B are two 1-forms such that 
9(X,61) = A(X) and g(X,f2) = B(X) (4) 
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with 1, €2 unit vectors which are orthogonal, i.e, 


g(€1, 2) = 0. (5) 
This type of manifold are denoted by G(QE)n. 


In [16], H.G. Nagaraja introduced the concept of N(k)-mixed quasi Einstein manifold and 
mixed quasi constant curvature.A non flat Riemannian manifold (MM, g) is called a N(k)-mixed 


quasi Einstein manifold if its Ricci tensor of type (0,2) is non zero and satisfies the condition 
S(X,Y) =ag(X,Y)+bA(X)B(Y) 4+ cB(X)A(Y), (6) 

where a, b,c are smooth functions and A, B are non zero 1-forms such that 
W(X, 61) = A(X) and g(X, 2) = B(X) VX, (7) 


with £1, € the orthogonal unit vector fields. Such aa manifold is denoted by the symbol N(k) — 
(MQE)n. 


The notion of hyper-generalized quasi Einstein manifold has been introduced by A.A.Shaikh, 
C. Ozgiir and A.Patra{17] in 2011. An n-dimensional Riemannian manifold (M",g) (n > 2) is 
called a hyper generalized quasi-Einstein manifold if its Ricci tensor of type (0,2) is non zero 
and satisfies the following condition 


S(X,Y) = ag(X,Y)+bA(X)A(Y) + A(X) B(Y) + A(Y)B(X)] 
+d[A(X)D(Y) + A(Y)D(X)] (8) 


for all X,Y € x(M), where a,b,c and d are real valued, non-zero scalars functions on (M”, g). 
A, B and D are non zero 1-forms such that 


9(X,&1) = A(X), g(X,&2) = B(X) and g(X,f3) = D(X), (9) 


wheer €1, £2,€3 are three unit vector fields mutually orthogonal to each other at every point on 
M.Here a,b,c,d are called the associated scalars, A, B,D are called the associated main and 
auxiliary 1-forms. We denote this type of manifold (HGQE),». 


§2. Preliminaries 


From (8) and (9), we get 
S(X, X) = a|X|? + b]g(X, &1)|? + 2e]9(X, 1) 9(X, €2)| + 2d]g(X, &1)g(X, &3)|, VX. (10) 


Let 6, be the angle between €; and any vector X; 02 be the angle between £2 and any 
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vector X; 43 be the angle between €3 and any vector X. Then 


cos 61 = G(X, &1) _ G(X, £1) 
V9(E1, €1) V GX, X) J Gg(X,X) 
as g(€1,€1) = 1), and 
cos 02 = 9X, &) and cos63 = GX §3) 


G(X, X) VIX, X) 


If b > 0, c > 0, we have from (10) 


(a+b+2c+2d)|X|? > alX|? + dlg(X,&)I? + 2elg(X, €1)9(X, &)| 
+2d|g(X, &1)9(X, €3)| = S(X, X). (11) 
Now, contracting (8) over X and Y, we get 
r=na, (12) 


where r is the scalar curvature. Since 1, €2 and &3 are orthogonal unit vector fields, therefore 


g(1,€1) = 1, g(a, 2) = 1, 9(€3, €3) = 1, o(£1, €2) = 0, g(f1,€3) = 0 and g(£2, £3) = 0. 


Again, putting X = Y = € in (8) we get S(£1,4;) =a+0. Putting X¥ = Y = & in (8) we 
get S(&2,2) =a. Putting X = Y = &3 in (8) we get S(&3,&3) =a. 


If X is a unit vector field, then S(X,X) is the Ricci-curvature in the direction of X. 


Notice that Q is the symmetric endomorphism of the tangent space at each point corre- 
sponding to the Ricci-tensor S, where 


g(QX,Y) =S(X,Y) V X,Y ETM. (13) 


Let I? denote the squares of the lengths of the Ricci-tensor S$. Then 


2 => S8(Qei,e%), (14) 
i=1 
where {e;}, i = 1,2,--- ,n is an orthonormal basis of the tangent space at a point of (HGQE),. 


Now from (8) we get 


S(Qei,e:) = ag(Qei, ei) + bA(Qe;) A(ex) + c[A(Qei) Ble) + Alex) B(Qei)] 
+d[A(Qe;) D(e;) + A(e;)D(Qe:)], 


? = (n—2)a? + (a+b)? + 2c? + 2d”. (15) 


These result will be used in the sequel. 
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§3. Ricci Semi-symmetric (HGQE),,(n > 3) 
Chaki and Maity proved that (QE),(n > 3) is Ricci Semi-symmetric if and only if 
A(R(X,Y)Z = 0. 
Let us suppose that (HGQE),, (n > 3) is Ricci-Semi symmetric. Then 
A(R(X,Y)Z =0. (16) 


From (16) we get 
A(Q(X)) = 0, (17) 


where Q be the symmetric endomorphism of the tangent space at each point corresponding to 
the Ricci tensor S$. Then 
WQX,Y) = S(X,Y). (18) 


Then from (8) we get 
A(Q(X)) = (a+ b|)A(X) + cB(X) + dD(X). (19) 
From (17) and (19) it follows that 
(a+ b)A(X) + cB(X) + dD(X) = 0. (20) 
Thus we can state the following. 


Theorem 3.1 If a(HGQE),, is Ricci Semi symmetric than (a+b)A(X)+cB(X)+dD(X)=0. 


§4. Sufficient Condition for a Compact Orientable (HGQE),(n>3) 
Without Boundary to be Isometric to a Sphere 


In this section we consider a compact, orientable (HGQE),, without boundary having constant 
associated scalars a,b,c and d. Then from (11) and (15), it follows that the scalar curvature is 
constant and so also is the length of the Ricci-tensor. 

We further suppose that (HGQE),, under consideration admits a non-isometric conformal 
motion generated by a vector field _X. Since /? is constant, it follows that 


£xl =0, (21) 


where £x denotes Lie differentiation with respect to X. 

Now, it is known ([2], [4], [5], [9], [12], [13], [14], [15]) that if a compact Riemannian manifold 
M of dimension n > 2 with constant scalar curvature admits an infinitesimal non-isometric 
conformal transformation X such that £x1? = 0 then M is isometric to a sphere. But a sphere 
is Einstein so that b, c and d vanish which is a contradiction.This leads to the following theorem. 
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Theorem 4.1 A compact orientable hyper generalized quasi Einstein manifold (HGQE),, (n > 


3) without boundary does not admit a non-isometric conformal vector field. 


§5. Killing Vector Field in a Compact Orientable (HGQE),(n >3) 
Without Boundary 


In this section, we consider a compact, orientable (HGQE),,(n >3) without boundary with 
a,b,c and d as associated scalars. 
It is known [4] that in such a manifold M, the following relation holds 
| [S(X, X) — |VX|? — (divX)?]du < 0 VX. (22) 
M 
If X is a killing vector field, then divX = 0 ([4]). Hence (22) takes the form 


| [S(X, X) — |VX|?]dv = 0. (23) 
M 


Let b> 0,c>0,d>0. Then by (11) 


(a+b+2c+2d)|X|? > S(X,X). (24) 
Therefore, 
(a+ 6+ 2c + 2d)|X|? —|VX|? > S(X, X)— |VX|?. (25) 
Consequently, 
| [(a + b+ 2c + 2d)|X|? — |VX|?|dv > ; [S(X, X) — |VX|?]dv, (26) 
M M 

and by (23) 

| [(a + b+ 2c 4+ 2d)|X|? — |VX|?]du > 0. (27) 
M 


Ifa+b+ 2c+ 2d <0, then 
[(a +b + 2c + 2d)|X|? — |VX|?]dv = 0. (28) 
M 
Therefore, X = 0. This leads to the following. 


Theorem 5.1 [fin a compact orientable (HGQE),(n >3) without boundary and the associated 
scalars are such that b >0,c >0,d>0 anda+b+ 2c+ 2d < 0, then there exists no non-zero 
killing vector field in this manifold. 


Corollary 5.1 If in a compact orientable (HGQE)n (n> 3) without boundary, and each of 
the associated scalars a,b,c,d, is greater than zero, then any harmonic vector field X in the 
(HGQE),, is parallel and orthogonal to one of the generators of the manifold which makes 
greatest angle with the vector X. 
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Abstract: In this paper, we are developing Quadrature Methods (numerical integration 
method) of continuous function f(z) on a compact interval [a,b] and deriving a polynomial 
Pm(a) of degree m such that integration of P(x) from a to b is equal to integration of f(x) 
from a to b. We are using least square method to fit the polynomial P(x). Also derive 
Newton-Cotes formulas and composite formula from this method, estimate errors and given 
MATLAB codes. 
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§1. Introduction 


With the advent of the modern high speed electronic digital computers, the Numerical Integra- 
tion have been successfully applied to study problems in Mathematics, Engineering, Computer 
Science and Physical Sciences. Numerical integration, also called Quadrature, is the study of 
how the numerical value of an integral can be found. The purpose of this paper is quadrature 


methods for approximate calculation of definite integrals 


b 
Ilf]= if f(a) (1.1) 


where f(x) is integrable, in the Riemann sense on [a,b]. The limit of the integration may 
be finite. Numerical integration is always carried out by mechanical quadrature and its basic 


scheme is as follows: 
fe) = 0 Afi t+ RIF], (1.2) 
where f; = f(a;) is continuous function in[a, b]. Ajand x; are called Coefficients(Weights) and 


nodes for Numerical Quadrature, respectively, and R[f] is error of Quadrature method. Once 


the coefficients and nodes are set down, the scheme (1) can be determined. 
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§2. Preliminaries 


2.1 Order of Quadrature Method 


Order of accuracy, or precision, of a Quadrature formula is the largest positive integer n such 


that the formula is exact for x", for each k = 0,1,-+- ,n. 


2.2 Error of Quadrature Method 


The integration (1.1) is approximated by a finite linear combination of value of f(x) in the form 
(1.2). The error of approximation of (1.2) is given as 


C 


Rn = aut ), (2.1) 


where € = (a,b), m > n is order of (1.2) and error constant of (1.2) is 


C= fp ors Aa: (2.2) 


2.3 Interpolation Polynomial 


Let f(x) be a continuous function defined on some interval [a,b], and be prescribed at n+ 1 
distinct tabular points 79, 71,--- ,@» such that a= % <4 <4 <--+ <4, =b. The distinct 
tabular points 2, 21,-°-- ,£n are equispaced, that is 7,41 — x7, =h, k= 0,1,2,---,n—1. The 
problem of polynomial approximation is to find a polynomial P,,(x), of degree< n, which fits 
the given data exactly, that is, 


P, (ai) = f(ai),t = 0,1,2,--+,n. (2.3) 


The polynomial P,,(x) is called the interpolating polynomial. The conditions given in 
(5)are called the interpolating conditions. 


2.4 Least Squares Interpolation Polynomial 
Let the polynomial of the m*” degree 
Pym (x) = a9 + a1% + Agu? + +++ + Am2™ 


be fitted to the data points (xj, f(z;)) i = 0,1,2,---,n, where m < n and ajs are satisfy the 


the system of equations 


(n+l)ag+a1 So aj +az Soa? +++ +am 5 0” = S° f(a), (2.4) 
i=0 i=0 i=0 i=0 
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n n n n 
2 1 
ag > tj ta1 oa +++ + Om So ort a eer 
i=0 i=0 i=0 i=0 


n n n n 
ao ) ve +ay ) BE be diy ) i ) ui” f (xi). 
i=0 i=0 i=0 i=0 
There are m+ 1 equations in m+ 1 unknowns. 


Lemma 2.1 Let P,,(x) be the least squares interpolation equation of f(x) on [a,b]. Then 


> Pmlai) & D> fla), (2.5) 
1=0 1=0 


where Lp = 4,tn, =b, 2 =atih andh= (b—a)/n. 


Proof Let the P(x) = ao + aya + agx? +-+++am2x™ is least squares interpolation 
equation of f(x) on [a,b]. Then Pp(ro) = ao + ait + a2x% +--+ + ama and P(x) = 
ag +4421 + aga? +++ +am2™, and so on Py(tn) = a9 + 1% + agx2 +-+++amz™. Adding 
all we get 


n n n n 
SO Plee) = (n+ 1)ag +01) aj+a2> af+-:-+am>_ a? 
i=0 i=0 


1=0 i=0 


apply Equation (2.4), we get 


n n 


S Pale)S SS fe). 


i=0 i=0 


Theorem 2.2 Let P,,(x) is least squares interpolation equation of the integrable function f (x) 
on finite interval [a,b] and 


JS Pla) & D> F(a), 
1=0 1=0 


where xo = a, Xn = b if and only if 


/ Py (a) da =) f(a)dz. (2.6) 
xo xO 
Proof Multiplying with h = (b— a)/n and take limit h — 0 in (2.5), we get 


a eee) Pe 


This completes the theorem. 


§3. Least Square Quadrature Method 


Consider the integral in the form (1.2) for each ¢ = 0,1,2,---,n. Now we are dividing the 


Mechanical Quadrature Methods from Fitting Least Squire Interpolation Polynomials 35 


interval [a, b] into n (finite) equal sub interval and take the nodes x's are equispaced points such 
that x; = xo+ih € [a,b], i = 0,1,2,--- ,n, where zp = a, 2, = band h = (b—a)/(n). So we have 
data points (x;, f(a;)) i = 0,1,2,--- ,n for fit a polynomial P,,(x) = ap +ayr+a2x7+-+-+am2"™. 
we have 


[- f(a)dax 


ll 
s—, 
3 
a 
Rae 
Q 
8 


2 2 3 3 m+1 m+l 
Di LE, vt, & x —— 
= do(%n — 0) + a> + ap Sa See ero Ge 


(3.1) 


This method is called L?, -Quadrature method(L?, — rule), here m is donate degree of 
polynomial and n is donate number of data points. To solve the least square Quadrature 
method we have at least m+1 points. Order of this method is greater then or equal to m, since 
it’s exact for polynomial of degree m. The error constant of (3.1) is 


Ln i gi — gi 
C= x —aot+ “ % £ a 
and error 
C 
— ~ ¢k) 


where k > m,a < € < b. Now following cases arise: 


Case 1. m=0, that is Po is a constant function. 


From (2.4) we have ao(n + 1) = 37}, f (xi) and a, = az =--++ = am = 0, substituting this 
values in (9) web get 


f(v)de = 22 = 89) ¥ p¢w,), (3.2) 


ro i=1 
Case 2. m=1, that is P, is a linear polynomial. 
From (2.4) we have 


n n n n n 
ag(n +1) +a1 9° 2% = So fiao (vita) a; = a 
i=0 i=0 i=0 i=0 


i=0 


and az = a3 =--: =Q@m = 0. Solving for a; and az we get 


(n+ 1) pa, i= (Yizo Zi)? 


(n+ 1) pte — Sao to fs . 
(2 +1) no BF — Oia 24)? 


ag = 


ay = 
After simplification we get 


2 7 me 
mG) |e ae) > fi—3(0 +n) x ifil 


ao = 
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ay = 


ainda) 
1=0 1=0 
Substituting this values in (3.1), and simplification we get 


at nh 
[fete = AY seo. 


This is same as m = 0. The method (3.2) is called L?- Quadrature method and the error 
constant of (3.2) is 


6 
nh(n + 1)(n + 2) 


7 ace. 7 ine . g shen? _— =(@n—@a)*? _ (b—a)° 
c= x dx mee = ee om 


and error of (3.2) is 
R= == a) pay (€) = ~——_ f )(€) (3.3) 
6n- 2! 12n i 
where 2p < € < ay. To solve this method, we have at least 2 data points and the order of (3.2) 


is 2. 
Case 3. m= 2, that is Pj is a polynomial of degree two. 


From (2.4) we have 


n n n 
(n+1)ao + a1 $5 ai +0250 a? Poe: = A, 
i=0 i=0 i=0 


n 


n n nm 
ao )) xj +a) >a} +a.) a3 = S "(20 +ih)f; = Aro + hB, 
i=0 i=0 i=0 


i=0 
n 


n n n 
ag) 2? +a, ) 23 +42) 24 = S\(ao + ih)? f; = Awd + 2Bhay + Ch?, 
1=0 1=0 1=0 1=0 


where A= 77) fi, B= Veo thi, and C = yy fi. we have a3 = a4 =--- = am = 0. 


Solving the three linear system of equation for ag,a, and ag by MATLAB, we get 


3 
ee (n + 1)(n3 + 4n? +n — 6)h2n 
x (3Ah?n* + 12Ahn329 — 12Bh?n3 — Ah?n? — 6Ahn?29 + 10An? x2 
+6Bh?n? — 64Bhn7x9 + 10Ch?n? — 2Ah?n — 6Ahnay — 10AnzZ 
+6Bh?n — 8Bhnao — 60BnaxZ — 10Ch?n + 60Chnap + 12Bhay + 60C' x2) 
a, = —{6(6Ahn? — 3Ahn? + 10An?2 — 32Bhn? — 3Ahn — 10Anz 


—4Bhn — 60Bnx + 30Chn + 6Bh + 60C'x)}/h?n(n? + 3n + 2)(n? + 2n — 3) 
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and 
30(An? — An — 6Bn + 6C) 


02 Fant + 5n8 + 5n® — 5n — 6) 


Substituting these values in (3.1), and simplification we get 


a _ hn(An? — An? + 6An + 30Bn — 6A — 30C) 
[/ 1eyae= Gn +ayn+ 2+) 


Substituting A,B and C we get 


(n> — n? + 6n — 6+ 30ni — 3077) f;. (3.4) 


Ln hn 
RO encencesicesp> 


This method is called L3-Quadrature method. To solve this method, we have at least 3 data 


points. 
Case 4. m=3, that is P3 is a polynomial of degree three. 


Following previous case we get the same as (3.3). The error constant of (3.4) is 


n 


“te hn 
C = i xdz — YS, (3 — 1? + 6n — 6 + 30ni — 3027) (a + th)! 
% (n—1)(n+3)(n + 2)(n +1) S| 


_ (8n? = 8n+18)n7h? _ (3n? — 8n+18)(@#n— 20)? ___ (38n? — 8n + 18)(b— a)” 


10) 


210 210n3 210n3 

The error of (3.4) is 

(3n? — 8n + 18)(b— a)? 
210n3 - 4! 

where a < € < b. The order of (3.4) is 4. 


i FOO, (3.5) 


Note 3.1 If m 20 is even number then L7, method same as L?,,, method. 


§4. Newton-Cotes Formulas from Least Square Method 


We can derive trapezoidal rule, Simpson 1-3rd rule and Simpson 3-8th rule from least square 
method. 


Taking n = 1 in (3.2) we get 
xy h 
[fear = Flo +). 


This formula is called trapezoidal rule. The error of trapezoidal rule is, from (3.3) 


—(b—a)? 


(2) 
SFO, a<E<b. 


R= 
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Taking n = 2 in (3.4) we get 


ao ee 
a f(x)dx = eer aaa \fi 


* (0fo + 40h + 10fs) = $Ufo + 4hi + 2) 


This formula is called Simpson 1-3rd rile rule. The error Simpson 1-3rd rule of is, from (3.5) 


_ =(b= 4)” a4) 
Aa (€), a@<€ <b. 


Similarly, Simpson 3-8th rule come from (3.4) with n = 3, that is 


| fla)de = (fo + Bf + 8f2 + fs) 


and error come from (3.5), R = —(3/80)h°fM(€),a<& <b. 


The weights of the integration method of (3.4) with equispaced point for n < 6 are given 
in Table 1. 


a 2/45 7 32 12327 a 11 26 31 26 11 
15 | 5/288 19 75 50 50 75 19 5/336 31 61 78 78 61 31 
f6 | 1/140 41 216 27 272 27 216 41 1/14 712151615127 


Table 1. Weight of Newton-cote rules and Weights of LY Quadrature Method 


aan y ico a f \t@ 
B (x), ii Y 1] ue j ' 
in Y \ | h 
| 7 \ y / | = 
i \ | / | 
We, \\ ; | | 
i \ F | 
a \ 


Figure 1 a,b,c 
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§5. Graphically Meaning of Least Square Integration Method 


Let the polynomial P,,(x) of degree m is fitted by least square interpolation method by using 
data points (a;, f;) i = 0,1,2---,n,. If m=1, take n is large number then the the polynomial 
P,(ax) is going to exact fit polynomial such that the area A+C=B (fib: 1(a)). That’s way the 
integration of P;(x) on [a, b] is gives exact value of integration of f(x) on [a,b]. Similarly P2(2) 
(or P3(a)) is best interpolation polynomial such that the area A+C=B, such as those shown in 
Figure 1. 


§6. Problems 


Problem 6.1 Find approximate value of 


3 
r= | sin(x)e*dx 
1 


fit a straight line y(x) such that oie y(a)dx = I. 


Solution Let f(x) = sin(x)e* and yy, be the straight line by fit n+1 data points (2;, f(2;)), 
i = 0,1,2,---,n. Now we divide the interval [1 3] into two equal subinterval, that is n = 2 
or h = 1. then 3 data points are (1, f(1)), (2, f(2)) and (3, f(3)). we fit a straight line y2 by 
normal equation (5) we get 
yo = 0.27" +3.4 


following this we get 
ya = 0.782 4+ 3.15, 


yg = 1.172 + 2.77 
yig = 1.39x% + 2.51 
y32 = 1.514% + 2.36 


and 
yea = 1.57x + 2.28. 


But we know if n — oo then fP Yn(x)dx ee f(x)dx. Therefore, I = fP(L.57e + 2.28)dx = 
10.84. 


Problem 6.2 Fit quadratic equation P2(x) such that 


1 1 
7 P(x)dxz = | sVa + ldx 
0 0 


and find approximate value of fe aVJa + ldz. 


Solution Let P2, be the quadratic equation by fit n equal space data points in [0, 1]. By 
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least square method we have 


Po,(x) = 0.378937382 + 1.03527618a + 3.61400724(E — 20), 
Po,,(x) = 0.378928452? + 1.039562852 — 0.00227848, 

Po, (x) = 0.378392732? + 1.04141576a — 0.00304322, 
P>,,,(%) = 0.3783134a? + 1.04167012 — 0.00314653. 


Let I, = is Po, (x)dx then Iz = 0.643950551, 41 = 0.643812428, I51 = 0.643795564 and 
T1091 = 0.643792992. The exact value of Hib xVx + 1dx upto five decimal is 0.64379. 


Problem 6.3 Find the approximate value of 


1 
1 
r= | dx, 
0 2+2 


using L} and LF rules with different equal subintervals. Using the exact solution, find the 


absolute errors. 


Solution Results for the L? and L3 rules to estimate the integral of f(#) = 1/(2+.2) from 
x = 0 to 1. The exact value is Ieract = Ve 1/(2+2)dx = In(x +2) = In(3) — In(2)=0.4054651. 
We get 


| s | o.4054801 | —_ 0.000150 
| 32 | 
| 64 | 
| 128 | 


6 0.0000084 
: 
z 
128 0.0000012 


co 


§7. Conclusion 


We develop this new method for easy to solve Definite Integral of finite interval with equispaced 
nodes and derived Simpson 1/3rd rule and Simpson 3/8th rule from DL} Quadrature Method. 
In this method (£4) weights are increasing from a to midpoint(i.e (a + b)/2) of interval and 
decreasing from midpoint to b. The advances is the weights of L} — method are positive 
(since(n?—n? +6n—6+30ni—30i7) > 0 for alln > 2 for alli) . We have given the MATLAB code 


also, give any continuous function f(a) on [a,b] that will be give an approximation integration 


value of f(a) from a to b. Also, we are developing this concept to high degree polynomials and 


high dimension. 
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Abstract: In this paper, we examine the motion of a robot end-effector by using the 
Blaschke approach of a ruled surface generated by a line fixed in the robot end-effector. In 
this way, we determine time dependent linear and angular differential properties of motion 
such as velocity and acceleration which play important roles in robot trajectory planning. 
Moreover, motion of a robot end-effector which can be represented by a right conoid and an 


additional parameter called spin angle is investigated as a practical example. 


Key Words: Blaschke frame, curvature theory, robot end-effector, robot trajectory plan- 


ning, ruled surface. 


AMS(2010): 53A05, 53A17, 53A25. 


§1. Introduction 


In robotics, a robot end-effector is a device at the end of a robotic arm. Robot end-effectors are 
widely used in transportation, welding industry, medical science, military and many other areas. 
Recently, they can be used in the research areas which have critical importance of accurate 
motion such as surgical operations and bomb disposal. So accurate trajectory planning of a 
robot end-effector becomes an important research area of robotics. In this area, one of the 
most interesting problems is determining time dependent differential properties of motion of a 
robot end-effector which are linear and angular velocities and accelerations. These differential 
properties play important roles in robot trajectory planning. 

As a robot end-effector moves on a specified trajectory in space, a line fixed in the end- 
effector generates a ruled surface [13]. There is an important relationship between time depen- 
dent properties of motion of the robot end-effector and differential geometry of the ruled surface. 
By using this relationship, Ryuh and Pennock proposed a method based on the curvature theory 
of a ruled surface generated by a line fixed in the end-effector to determine linear and angular 
properties of motion [12, 13, 14]. After that, this research area was also studied in Lorentzian 
space. Ekici et al. examined motion of a robot end-effector in Lorentzian space by using the 
curvature theory of timelike ruled surface with timelike ruling [7]. Ayyildiz and Turhan also 
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determined differential properties of motion of a robot end-effector whose trajectory is a null 
curve [3]. 

On the other hand, there is also an efficient relationship between directed lines and dual 
unit vectors. This relationship known as “E. Study mapping” or “transference principle” which 
can be stated as: “there exists one-to-one correspondence between the directed lines in line 
space and dual unit vectors in dual space” [11, 16]. By the aid of this correspondence, W. 
Blaschke defined a frame called Blaschke frame on a ruled surface by taking directed lines pass 
through striction curve of the ruled surface instead of real unit vectors used in Frenet frame 
of ruled surface. He also gave some invariants which characterize the shape of a ruled surface. 
Several authors used Blaschke frame in their researches concerning with kinematics, spatial 
mechanisms and many other areas [1, 2, 18]. 

In this paper, we use the relationships between kinematic, ruled surfaces and dual vector 
algebra. First, we represent motion of a robot end-effector on a specified trajectory in space 
as a ruled surface generated by a line fixed in the end-effector and an additional parameter 
called spin angle. We define a dual frame called dual tool frame on robot end-effector in order 
to obtain a relationship between Blaschke frame of the ruled surface, which is used to study 
the differential geometry of a ruled surface by means of dual quantities, and time dependent 
differential properties of robot end-effector. By using this relation, we determine time dependent 
differential properties of motion of a robot end-effector which are linear (translational) and 
angular (rotational) velocities and accelerations. These differential properties have important 
roles in robot trajectory planning. In this method, we use just a dual vector called dual 
instantaneous rotation vector of dual tool frame to determine the differential properties. So, this 
method has more advantages than traditional methods which based on matrix representations 


in terms of being simple and systematic. 


§2. Preliminaries 


In this section, we give a brief summary of basic concepts for the reader who is not familiar 
with dual numbers, dual vectors and dual space. 

As introduced by W. Clifford, a dual number can be defined as @ = a+ €a*, where a and 
a* are real numbers and called real part and dual part of dual number 4G, respectively, and ¢ is 
dual unit which satisfies the condition e? = 0, [17]. The set of all dual numbers can be denoted 
by D. Addition and multiplication of two dual numbers @ = a + ca* and b = b + eb* can be 
defined as 

a@+b=(a+5b)+<e(a* +8") 


and 
ab=ab+e(ab* + a*b) 


respectively [4, 10]. The set D is a commutative ring, not a field. A function of a dual number 


f(@) can be expanded in a Maclaurin series as 


f@) = flat ea") = f(a) + ea" f'(a), 
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where the prime indicates derivation of f(a) with respect to a [5]. 


A dual vector can also be defined as G= a+ea*, where a and a* are three dimensional 


vectors in real space and e? = 0. The set of all dual vectors is a module over the ring D and 


is called dual space or D-module, denoted by D°, [15]. Dual scalar and vector products of two 


dual vectors @= at+ea* and b= b+eb* can be defined as 


(a,b) = (a,b) +¢((a,b*) + (a*, bY) 


and 
axb=axb+e(ax b* +a* xb) 


respectively [16]. The norm of a dual vector @ can also be given by [10, 17] 


* 


(a, a*) 


I|a|| = llall + lal” 


(a # 0). 
If ||@|| = 1, then @ is called a dual unit vector. The set 

S? = {a= atea* | ||@||=1; a, a* € R*} 
is called dual unit sphere. 


Theorem 2.1([8]) The set of all directed straight lines in R® are in one-to-one correspondence 
3 


4 


with the set of all points of the dual unit sphere in 


A dual angle between two oriented lines in three dimensional real space can be defined as 
6 = 0+.e 6*, where @ and * are the real angle and the shortest distance between these lines, 
respectively, [4]. 


§3. A Robot End-Effector and its Dual Tool Frame 


In this section, we introduce tool frame of a robot end-effector which consists of three mutually 
perpendicular unit vectors described by Ryuh and Pennock [13] in detail. Then, we represent 
motion of a robot end-effector by using a ruled surface generated by a line fixed in the end- 
effector and an additional parameter called spin angle. By taking three lines instead of three 
unit vectors, we define a dual frame called dual tool frame on robot end-effector which will be 
used to study the motion. 

The tool frame consists of three orthogonal unit vectors strictly attached to robot end- 
effector. These are; orientation vector O which is a unit vector in the direction of the gripper 
motion as it opens and closes, approach vector A which is a unit vector in the direction normal 
to the palm of robot end-effector, and normal vector N which is a unit vector in the direction 
perpendicular to the plane of the gripper (see Figure 1), [12]. The origin of the tool frame 
is called tool center point and denoted by TCP. By using tool frame and tool center point, 
location and orientation of a robot end-effector can be described completely. 
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As a robot end-effector moves on a specified trajectory in space, a line called tool line 
fixed in the end-effector which passes through TCP and whose direction vector is parallel to 
the orientation vector O generates a ruled surface [12]. This ruled surface can be expressed as 


X(t,v) = a(t) + v u(t), 


where a is the specified trajectory which robot end-effector follows (directrix of the ruled sur- 
face), u is a unit vector called ruling parallel to the orientation vector O, t is the parameter of 


time, and v is an arbitrary parameter. 


During motion, the approach vector A may not be always perpendicular to the ruled 
surface. As seen in Figure 1, there may be an angle between the approach vector A and the 
surface normal vector on the directrix which is denoted by S,. This angle is called spin angle 
and denoted by 7 [12]. Thus, a robot end-effector motion which has six degrees of freedom in 
space can be completely described by a ruled surface generated by a line in robot end-effector 
which provides five independent parameters and a spin angle. 


Specified 
trajectory 


Tool center point (TCP) ' | 


Rulings of ruled surface 


Figure 1 Robot end-effector and spin angle 


Tool center point 


Figure 2 Dual tool frame of a robot end-effector 
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Now, we define dual tool frame by taking three directed lines instead of three unit vectors of 
the tool frame. These lines pass through the TCP and their direction vectors are the orientation 
vector O, the approach vector A and the normal vector N, respectively. From Theorem 2.1, 
these lines correspond to three dual unit vector which can be called dual orientation vector, 
dual approach vector and dual normal vector and can be denoted by O, A and N , respectively 
(see Figure 2). 


§4. Blaschke Approach to the Motion 


In this section, we give Blaschke frame of a ruled surface generated by a line fixed in the robot 
end-effector. By relating Blaschke frame and dual tool frame, we determine linear and angular 
differential properties of motion. Furthermore, we give corollaries for some special cases of 


motion. 


From Theorem 2.1, it can be said that a ruled surface can be represented by a dual unit 
vector based on a real parameter. So, we can consider the ruled surface generated by motion 
of robot end-effector as a dual unit vector ti(t) = u(t) + cu*(t), where wu is ruling of the ruled 
surface, u* is moment vector of u about the origin, t is the parameter of time, and <? = 0. 
The moment vector can be found as u* = c x u, where c is striction curve of the ruled surface 
satisfies the condition that (c’, u’) = 0, [6]. In this paper, we consider the case without u(t) = c1 
which means ruled surface is a cylinder and u*(t) = c; which means ruled surface is a cone, 
where c; is a constant. In order to simplify formulations, arc-length parameter of the striction 
curve denoted by s can be used instead of the parameter of time ¢ and it can be obtained as 


The Blaschke frame of a ruled surface is defined on striction curve and it consists of three 


orthogonal dual unit vectors given as follows [4]: 


ae 
a a Pe Uy ap ae i 
Uy=U, UW=—, UZ= U1 X U2, 
Pp 
where p = p+ ep* = ||a‘||, @2 and Gs are normal line and tangent line of the ruled surface on 


the striction curve, respectively, and the prime indicates the derivation with respect to s, [4]. 
The derivative formulae of Blaschke frame can be given as 


Uy 0 i) 0 uy 

d 

de U2 = —p 0 qd U2 3 (1) 
U3 0 —-q U3 


det(ti1, w’1, u”’1) 


eI” 
characterize the shape of a ruled surface. If p = 0, ruled surface is a cylinder; if p* = 0, ruled 


where Gg = qt+eq = p and g which are called the Blaschke’s invariants 
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surface is a developable ruled surface which is a surface that can be flattened onto a plane 
without distortion; if g = 0, all rulings of ruled surface are parallel to a plane; if g = 0, ruled 
surface consists of binormal vectors of a curve, [4]. 

Let 6 = y+e y* be a dual angle between dual unit vectors A and tiz, where p=nto 
is real angle, where 7 is the spin angle mentioned in Section 3 and o is an angle between two 
normal vectors of ruled surface, one is on the directrix and other is on the striction curve, and 


/ / 
ie 7 (see Figure 3). 


* is the shortest distance from striction curve to directrix, i.e., p* = 


i [lu 


Line of striction 


it, 


Specified trajectory (directrix) 


! Tool center point (TCP) 


Figure 3 Dual angle between the dual unit vectors A and U2 


By the aid of dual angle g, we can give dual tool frame relative to Blaschke frame in matrix 


form as 
O 1 0 0 Uy 
Al=]|0 cosg sing fig | - (2) 
N 0 -sing cos¢@ U3 


By differentiating equation (2) and substituting equation (1) into the result, we have 


O' 0 7) 0 iy 
A' |=| -peos@ —ésing dcos@ tia |; 
N’ psing —dcos@ —dsing tg 


where 6 = g' + q. By using equation (2), derivative formulas of the dual tool frame can be 


obtained in terms of itself in matrix form as 


O! 0 pcos@ —psing O 
A’ | =| -—pcos@ 0 5 A 
N’ psng —6 0 N 


From the above matrix equality, dual instantaneous rotation vector of the dual tool frame 
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which plays an important role to determine both linear and angular differential properties of 
motion of a robot end-effector can be obtained as 


tio =5 O+psing A+ peos@ N. 


By using equation (2), the dual instantaneous rotation vector can also be expressed in 
terms of the Blaschke frame as 
Wo = 6 ty 4p ti. (3) 
This dual vector is similar to dual Pfaff vector in terms of playing role in motion. The dual 
Pfaff vector is considered as dual velocity vector in dual spherical motion (see ref. [9]). So, 
we can consider the dual instantaneous rotation vector of dual tool frame wo as dual velocity 
vector of the motion of robot end-effector. 
w 
The dual tool frame attached to robot end-effector moves along unit direction —? | with 


|| Wo| 
dual angle ||to||. This dual motion contains both rotational and translational motion in real 


space. The real and dual parts of the dual vector wo correspond to instantaneous angular 
velocity and instantaneous linear velocity, respectively. By separating equation (3) into the real 
and dual parts, these velocity vectors can be found as follows 


wo = du, +p Us, (4) 


and 
Wo = Ou + O*u, + pus + p*us. (5) 


In order to find dual acceleration vector of the motion, we should differentiate dual velocity 
vector. By differentiating equation (3) and using equation (1), the dual acceleration vector can 
be obtained in terms of the Blaschke frame as 


Wo = 6 t+ PP te + pis, (6) 


where the prime indicates differentiation with respect to s. By separating equation (6) into 
the real and dual parts, instantaneous angular acceleration vector and instantaneous linear 


acceleration vector can be found as 
Wo = 6 ut ¢'p us + p'ug (7) 


and 
wo =O upto urty'’p us + (y'p* + ¢% p)us + p'us +p” us, (8) 


respectively. Thus, linear and angular velocities and accelerations which are important dif- 
ferential properties of motion of a robot end-effector are found in terms of the parameter s 
which is the arc-length parameter of striction curve of the generating ruled surface. In order 
to determine time dependent differential properties, the vectors given in equations (4), (5), (7) 
and (8) should be related to the parameter of time. Now, we give time dependent linear and 
angular differential properties of motion of a robot end-effector as corollaries. 
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Corollary 4.1 Let the motion of a robot end-effector be represented by a ruled surface X(t, v) = 
a(t) +v u(t) and a spin angle n, where a is specified trajectory of robot end-effector, u is a unit 
vector parallel to the orientation vector O, and t is the parameter of time. Angular and linear 


velocities of robot end-effector can be given, respectively, as 

VA = Wo § (9) 
and 

UL = WO 5, (10) 


where wo and we are given by equations (4) and (5), respectively, and the dot indicates differ- 
ds 


dt. 
Corollary 4.2 Let the motion of a robot end-effector be represented by a ruled surface X(t, v) = 


entiation with respect to time, t.e., § = 


a(t) +v u(t) and a spin angle n, where a is specified trajectory of robot end-effector, u is a unit 
vector parallel to the orientation vector O, and t is the parameter of time. Angular and linear 
accelerations of the robot end-effector can be given, respectively, as 


aa = wo §+wo8? (11) 
and 

an = Wo 3+wi 8, (12) 
where Wo and ws are as given by equations (7) and (8), respectively. 


Now, we consider some special cases of motion of a robot end-effector and give some 
corollaries about these cases. 


Case 1. As a robot end-effector moves on a specified trajectory in real space, spin angle 7 
may be constant. Then, the derivative of the spin angle is equal to zero. For this case, by 
substituting the value of spin angle into equations (4), (5), (7), and (8), and by rearranging 
these equations, we can give time dependent linear and angular differential properties of the 
motion of a robot end-effector as in the following corollaries. 


Corollary 4.3 Let the motion of a robot end-effector be represented by a ruled surface X(t, v) = 
a(t) +u u(t) and a spin angle n, where a is specified trajectory of the robot end-effector, u is 
a unit vector parallel to the orientation vector O, and t is the parameter of time. If the spin 


angle n is a constant, then angular and linear velocities of robot end-effector can be given as 
va = ((o' +4) 1 +p us) § 


and 


ut = ((o' + q)uj + O° ur + pus + p*us) 8, 
respectively. 


Corollary 4.4 Let the motion of a robot end-effector be represented by a ruled surface X(t, v) = 
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a(t) + v u(t) and a spin angle n, where a is the specified trajectory of the robot end-effector, u 
is a unit vector parallel to the orientation vector O, and t is the parameter of time. If the spin 
angle n is constant, then angular and linear accelerations of robot end-effector can be respectively 
given as 


aa = ((o'+ q)u1+pus) &+((o" +4’) uw + o'pus + p'us) 8°, 


ap = ((o'+ qui + our + pus + p*us) § 


* 


H((o" +q') ut +. 5% uy + o'p ud + (o'p* + yp pus + p'us + p* us)3?. 


Case 2. A specified trajectory which robot end-effector follows may be striction curve of ruled 
surface generated by a line fixed in the robot end-effector. Namely, directrix and striction curve 
of generating ruled surface may be the same curve. Then, the angle o which is the angle between 
two normal vectors on directrix and striction curve and the distance between these curves are 
equal to zero. For this case, by rearranging equations (4), (5), (7), and (8), we can give time 
dependent linear and angular differential properties of the motion of a robot end-effector as in 
the following corollaries. 


Corollary 4.5 Let the motion of a robot end-effector be represented by a ruled surface X(t, v) = 
a(t) +u u(t) and a spin angle n, where a is specified trajectory of the robot end-effector, u is a 
unit vector parallel to the orientation vector O, and t is the parameter of time. If the specified 
trajectory is also the striction curve of the ruled surface, then angular and linear velocities of 


robot end-effector can be given as 
va = ((1/ +9) m1 +p us) § 


and 


respectively. 


Corollary 4.6 Let the motion of a robot end-effector be represented by a ruled surface X(t, v) = 
a(t) + v u(t) and a spin angle n, where a is specified trajectory of robot end-effector, u is a 
unit vector parallel to the orientation vector O, and t is the parameter of time. If the specified 
trajectory is also the striction curve of the ruled surface, then angular and linear accelerations 


of robot end-effector can be given as 


2 


tq’) uy + 17'puz + p'us) $ 
and 


ap = (( +q)up t+ q*ur + pug + p*ug) 8 
H((n” +4’) ut tq ur +7/p ud + n/p*us + p'u§ + p* us) 8?, 


respectively. 
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Case 3. Ruled surface generated by a line fixed in a robot end-effector may be a developable 
ruled surface (except a cylinder and a cone). So, the dual part of Blaschke’s invariant p is equal 
to zero, i.e., p* = 0. For this case, by making the necessary arrangement in equations (4), (5), 
(7), and (8), we can give time dependent linear and angular differential properties of the motion 
of a robot end-effector as in the following corollaries. 


Corollary 4.7 Let the motion of a robot end-effector be represented by a ruled surface X(t, v) = 
a(t) +v u(t) and a spin angle n, where a is specified trajectory of robot end-effector, u is a unit 
vector parallel to the orientation vector O, and t is the parameter of time. If the ruled surface 


is developable, then angular and linear velocities of robot end-effector can be given as 
va = ((n/ +4) u1 +p us) 8 


and 
vt = (7! + qhuy + O° ur + pus) 8, 


respectively. 
Corollary 4.8 Let the motion of a robot end-effector be represented by a ruled surface X(t, v) = 
a(t) +v u(t) and a spin angle n, where a is specified trajectory of robot end-effector, u is a unit 


vector parallel to the orientation vector O, and t is the parameter of time. If the ruled surface 


is a developable, then angular and linear accelerations of the robot end-effector can be given as 


2 


aa = ((n' + q) ui +p us) &+((n" +’) ui +7 puz + p'us) 8 
and 


v1 


tq’) uy + OY uy + np us + yp ug + p'us)s”, 


a = ((n! + q)uy +O" us + pus) §+((n 


respectively. 


§5. An Example 


Let the motion of a robot end-effector be represented a right conoid given by the equation 
X(t, v) = (vcost, vsint, 2sint) and aspin angle 7, where t is the parameter of time (see Figure 
4). Directrix and ruling of the right conoid are a(t) = (0, 0, 2sint) and u(t) = (cost, sint, 0), 
respectively. Since (a’,u’) = 0, directrix and striction curve of the ruled surface are the same 
curve, i.e., c= a. The right conoid can be expressed as a dual unit vector 


ai(s) = u(s) + eu*(s) = (cost, sint, 0) + e(—2sin?t, sin 2t, 0) 


where s is the arc-length parameter of striction curve. The first dual unit vector of Blaschke 
frame is ti(s) = &(s). The second and third dual unit vectors of Blaschke frame can be found 
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as 


tig(s) = (—sint, cost, 0) + e(—sin2t, —2sint, 0) 


and 
tu3(s) = (0, 0, 1), 


respectively. The Blaschke’s invariants can be obtained as p = p+ ep* = 1+ €2cost and 
G=q+teq*X =0+e0. Let 6 = y+ ey* be a dual angle between dual unit vectors A and Zi, 
where y and y” are the real angle and the shortest distance between the lines correspond to the 
dual vectors A and éz, respectively. Since directrix is also striction curve, the distance between 
these curves equals to zero, i.e., y* = 0, and the angle between two normal vectors on directrix 
and on striction curve equals to zero, i.e., g = 0. Thus, we have 6 = 7+¢0. Dual instantaneous 
rotation vector of dual tool frame can be found as 


Wo = wo + ews = (n' coss, 7’ sins, 1) + e(—2n' sin? s, ’ sin 2s, 2coss). 


Angular and linear velocities of the robot end-effector can be obtained by substituting wo and 
wo into equations (9) and (10), respectively. By differentiating the dual instantaneous rotation 
vector, we get 


~ , , . . 
Wo =Wotewg = (7 coss—7'sins, 7 sins +17 coss) 


+e(—2n" sin? s — 2n/ sin 2s, 7" sin 2s + 2n’ cos2s, —2sins). 


Angular and linear accelerations of the robot end-effector can also be obtained by substi- 
tuting wo and w%, into equations (11) and (12), respectively. 


Figure 4 Motion of a robot end-effector 
which can be represented by a right conoid and a spin angle 7 


Blaschke Approach to the Motion of a Robot End-Effector 53 


§6. Conclusions 


In this paper, time dependent differential properties which are linear and angular velocities and 
accelerations of the motion of a robot end-effector are determined by using Blaschke approach 
of a ruled surface generated by a line fixed in the end-effector. These differential properties are 
important information in robot trajectory planning. By the aid of Blaschke approach which 
uses dual numbers and dual vectors as basic tool, both linear and angular differential properties 
can be determined. This is achieved only by using a dual vector which is dual instantaneous 
rotation vector of dual tool frame. Thus, Blaschke approach presents a simple and systematic 
method to study motion of a robot end-effector without redundant parameter. This paper 
does not contain a computer program which compares Blaschke approach and conventional 
method of scalar curvature theory of ruled surfaces in real space. This is the subject of ongoing 
research works. However, it is believed that the presented method based on Blaschke approach 
will reduce computation time in computer programming for determining differential properties 


of motion and contribute to research area of robot trajectory planning. 
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§1. Introduction 


Let G(V,£) be a simple finite graph. The order of G is the number of vertices of G. A set 
S C V is a dominating set if every vertex v € V — S is adjacent to at least one vertex in S. 
The domination number of G, denoted by 7(G), is the minimum cardinality of the dominating 
sets in G. Generally, a dominating set S is said to be a Smarandachely k-dominating set if each 
vertex of G is dominated by at least k vertices of S. Let D(G,i) be the family of dominating 
sets of G with cardinality 7 and let d(G,i) = |D(G,)|. The polynomial 


IV(G)| 


D(G,t)= S> dG, ia" 


1=7(G) 


is defined as domination polynomial of G. For more information on this polynomial the reader 
may refer to [8]. A root of D(G,x) is called a domination root of G. It is easy to see that 
the domination polynomial is monic with no constant term. Consequently, 0 is a root of every 
domination polynomial (in fact, 0 is a root whose multiplicity is the domination number of the 
graph). 


§2. d-Number 


In this section we mainly focus on the number of real domination roots of some specific graphs. 


So we introduce a new definition as follows. 


1Received March 8, 2018, Accepted August 11, 2018. 
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Definition 2.1 Let G be a graph. The number of distinct real domination roots of the graph 
G is called d-number of G and is denoted by d(G). 


Theorem 2.1 For any graph G, d(G) > 1. 


Proof It follows from the fact that 0 is a domination root of any graph. 


Theorem 2.1 Jf a graph G consists of m components G1, G2,...,Gm, then 


d(G)< 4 a(Gi)—m+1. 


uM: 


t=1 


Proof It follows from the fact that D(G,x) = [] D(Gi, x). 


i=. 


Theorem 2.3 If G and H are isomorphic, then a(G) = d(H). 


Proof It follows from the fact that if G and H are isomorphic, then D(G,x) = D(H, 2). 


Theorem 2.4 If G has exactly two distinct domination roots, then d(G) = 2. 


Proof It follows from the fact that 0 is a domination root and complex roots occurs in 


conjugate pairs. 


Theorem 2.5 Let G be a graph without pendent vertices. If G has exactly three distinct 


domination roots, then a(G) = 1. 


Proof It follows from the fact that with the given conditions in theorem, Z(D(G,)) C 


{0, -2+ ivV2, =34iV3) (8). 


Theorem 2.6 For all n we have the following : 


a(Ky) 1; ifn is odd, 
‘ 2 ; ifn is even. 


Proof We have known the domination polynomial of K,, is 


D(Kn, x) = (14+ a2)" -1. (1 


The result follows from the transformation y = 1+ x in equation (1). 


Theorem 2.7 For any graph G, d(Go K,) =2. 


Proof Notice that D(Go ky,x) = x"(a + 2)” ([8]), where n is the order of G. Therefore 
d(G fe) ky) =. 


Theorem 2.8 For any graph G, a(Go Kz) =3. 
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Proof Notice that D(G o Ko,x) = x? (x? + 3x + 1)? ((8]), where n is the order of G. 
Therefore Z(D(G,x)) = {0, =3v5) This implies that d(G o Kz) = 3. 


Theorem 2.9 For all n the d-number of the star graph Sy, is 


a(S,) = 2 ; ifn is odd, 


3. ; of n is even. 


Proof We have known the domination polynomial of S;, is 


D(Sn,2) =2(1 +2)" +2". (2) 


Therefore it suffices to prove that f(x) = (1+ 2)" +2"~? has exactly one real root if n is 
odd and two real roots if n is even. But the number of real roots of f(x) is equal to the number 
of real roots of g(x) = (1+ 1)" + 4. Again the number of real roots of g(x) is equal to the 
number of real roots of g(+) = (1+a)" +2. Consider (5-7) =y"+y-—1, we find the number 
of real roots of h(y) = y” +y— 1. We have h(0) = —1 < 0 and A(1) = 1 > 0. Therefore by 
the intermediate value theorem, h(y) has at least one real root in (0,1). Also by De Gua’s rule 
[11] for imaginary roots, there are at least n — 1 complex roots for odd n and there are at least 
n — 2 complex roots for even n. Therefore we can conclude that h(y) has exactly one real root 
for odd n and two real roots for even n. It remains to show that all the real roots of f(a) are 
distinct. Suppose a € R is a double root of f(x). Whence, 


(1+a)"+a”'=0, (3) 
n(1 +a)""1 + (n—1)a"* =0. (4) 
From equation (3) we get 
qt 
eg NSS aoe ( since a # —1). (5) 


Putting the value of (1+ .a)"~! in (4) and simplify, we obtain a = n—1. Which is a 


contradiction since a < 0. 
Theorem 2.10 For all n the d-number of Kanan ts 1. 
Proof Notice that the domination polynomial of Kenan is 
D(Konan,t) = (1+ 2)?” — 1)? + 20”, (6) 


Suppose for a € R, ((1 +a)?” — 1)? + 2a?” = 0, then ((1 +a)?" 1)° = —2a?". But this 
is true only if a = 0, hence d(Kan2n) = 1. 


Theorem 2.11 The d-number of Kan+1,2n+1 is greater than or equal to 3 for all n. 
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Proof We have known the domination polynomial of Kon+1,2n+1 is 
D(Kansiantis) = (+2)? = 1)? + 2041, (7) 


It is easy to verify that 


1 1 1 
D (Hensransn 5) Spit Go - 1) > 0 


D (Konsinti,-l) =-1<0 
D (Kon41,2n41; -2) = 2?(1 — 27") < 0 
D (Kan41,2n+1,—3) = (277+! + 1)? -2 x 3rt1 > 0 


Therefore by the intermediate value theorem, Kon+41,2n+1 has at least one real domination 
root in (—1, —3) and at least one in (—3,—2), hence d(Kan+1,2n41) > 3. 


The Dutch-Windmill graph G3 is the graph obtained by selecting one vertex in each of n 
triangles and identifying them. 


Theorem 2.12 For n> 2 the domination polynomial of the Dutch-Windmill graph G3 is 


D(G8,2) = 2(1 +2)?" + (22 + 2)". 


Proof Let v be the center vertex of G3. It is clear that {v} is the only dominating 
set of cardinality 1. Therefore 7(G3) = 1 and d(G¥,1) = 1. The number of ways of selecting 
dominating set of cardinality which containing the center is (ee Also there are 2” dominating 
sets of cardinality n which does not contain the center vertex v. Similarly there are (") 2s 
ways to select a dominating set of cardinality n +7 which does not contain the center vertex v. 
Therefore D(G%, x) = a(1+ a)?" + (2u + 27)”, 


Theorem 2.13 For all n the d-number of the Dutch windmill graph Gat is 1. 
Proof We have known the domination polynomial of the Dutch windmill graph G3,,,, is 
D(G3"*" 2) =a2(1 +2)*"*? + (Qe + 2?)?"41, 


Suppose there is a number a € R with a £4 0 such that a(1 + a)4"*? + (2a + a7)?"*! = 0. 


Then we have a < 0 and by a simple calculation we have 


+--(-a) ° 


Suppose —2 < a < 0, then the left side of the equation (7) is negative but the right side is 
positive, a contradiction. Now suppose a < —2. Then the left side of the equation (7) is less 
than or equal to —2 but the right side is greater than —1, a contradiction. Therefore there is 


no nonzero real domination root for G3"! and hence d(G3"*") = 1. 
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Theorem 2.14 The d-number of G3” is greater than or equal to 3 for all n. 


Proof Notice that the domination polynomial of the Dutch windmill graph G3” is 
D(G?", 2) = a(1 +2)" + (Qe +27)*". 


It is easy to verify that D(G3",—1) > 0 and D(G2",—2) < 0. Also if a is a negative real 
number near to 0, then D(G3",a) < 0. Therefore by the intermediate value theorem, we have 
G2” has a real domination root in (—2,—1) and a real domination root in (—1,0) and hence 
a(Ge 3. 


The lollipop graph L,,,; is the graph obtained by joining a complete graph K,, to a path 
P,, with a bridge. 


Theorem 2.15 For n> 2 the domination polynomial of the lollipop graph Ly, is 


D(Lyn1,0) =a ((1 +a)" +(1+2)""*-1). 


Proof Let {v1,v2,--:,Un} be the vertices of the complete graph K,, and v be the path 
P, and let v is adjacent to v1. Clearly, y(£n,1) = 1 and d(£y1,1) = 1. For 2<i<n-1, 
the only non dominating sets of i vertices of Ly,1 are the subset of {v2,v3,--- , Un}. Therefore 
d(Ln.1,4) = ("5") — ("7"). Also d(Zn1,n) =n+ land d(Lni,nt+1)=1. 


Theorem 2.16 For all n > 2 the d-number of the lollipop graph Ly, is 


a(Ln.1) 2 ifn is odd, 
nig li) = 
3 ifn is even. 


Proof By Theorem 2.15 it is enough to prove that f(y) = y” + y”~! —1 has only one real 
root if n is odd and has exactly two real roots if n is even. By De Gua’s rule for imaginary 
roots, there are at least n — 1 complex roots if n is odd and there are at least n — 2 complex 
roots if n is even. Now, f(0) = —1 < 0 and f(1) = 2 > 0 for all n and f(—1) = —1 < 0 and 
f(—2) = 2”-!—1> 0 for all even n. Therefore by the intermediate value theorem, we have the 


result. 


The generalized barbell graph By,»,1 is the simple graph obtained by connecting two 
complete graphs K,, and K,, by a path P,. 


Theorem 2.17 Form <n, the domination polynomial of generalized barbell graph Bm ni 18 


Proof Let V = {v1,v2,-++ ,Um} and U = {uy, u2,--- , Un} be the vertices of By nj such 
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that if 2 A 7 every vertices V are adjacent, every vertices U are adjacent and v,, and uy, is 
adjacent. There is no one element dominating set and {v;, u;} is a dominating set of cardinality 
2 of Brin. Therefore ¥(Bmn,1) = 2 and d(Byn.1,2) = mn. Also observe that for 2 <i < 2n, 
a subset S' of vertices Bm »,1 of cardinality i is not a dominating set if either SC V or SCU. 
Therefore d(Bmni,t) = 7”) - (") - ee for 2 <4 < m, d(Bn»ni,i) = CP) - (2): for 
mt+1<i<nand d(Byna,i) = (ear forn+1<1i< 2n. This implies that D(Bmni,2) = 
(1+a)™"—1)[+a2)"-1]. 


Theorem 2.18 For all m,n the d-number of the generalized barbell graph Bm.n,1 1s 


1 if both m and n are odd, 
A( Brinn) —= 
2 otherwise. 


Proof The result follows from the transformation y = 1+ in the domination polynomial 
of Bmjn- 


The n— barbell graph B,,; is the simple graph obtained by connecting two copies of 
complete graph K,, by a bridge. 


Corollary 2.19 The domination polynomial of the n-barbell graph By. ts 


D(Bn1,£) = ((14+2)"—1)?. 


Proof It follows from the fact that the n-barbell graph B,,; and the generalized barbell 
graph B,,»,1 are isomorphic. 


Corollary 2.20 For all n, the d-number of the n-barbell graph By ts 


1 ifn is odd, 
d(Bn,1) = 
2 ifn is even. 


A bi-star graph Bin) is a tree obtained from the graph K2 with two vertices u and v by 
attaching m pendant edges in u and n pendant edges in v. 


Theorem 2.21 The domination polynomial of the bi-star graph Bom, n) 18 
D(Bonn)s%) =o” + a072(1t a)" part (1+ a)" +a (142). 
Proof Let {u,v}, U = {ur, u2,-++ , tn} and V = {v1, v2,--+ ,Um} be the vertices of Bm n 


such that u and v are adjacent, every vertices U are adjacent to u and every vertices V are 
adjacent to v. Clearly there is no one element dominating set. The set {u,v} is the only 
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dominating set of cardinality 2 of By». Therefore y(Byn) = 2 and d(Bmn,2) = 1. For 
3 <i<™m, the i—element dominating set of Bm, must contain {u,v}, and the i — 2 elements 


can have Cas) choice. For m+1<i<n, there are Cs i—element dominating set of Bm in 


containing {u,v} and (,_,”_,) i-element dominating set of By,» containing V{u}. For n + 


t—m—1 
1<i<m-+n-1, there are (eee i—element dominating set of By,» containing {u, v}, (,_,”_,) 


i—element dominating set of Bm n containing VU{u} and (,_™ 


ued) i—element dominating set 


of Bm» containing UU{v}. For i = m+n, there are (*;’) (m +n)—element dominating 
set of By» containing {u,v}, n (m+n)—element dominating set of Bm» containing V U{u}, 
m (m+ n)—element dominating set of Bm» containing UU{v} and one (m + n)—element 
dominating set of By.» not containing {u,v}. Also d(Bmn,m+n+1) = m+n-+2 and 


d(Bnn,m+n+2) =1. That is, 


1 if¢=2,m+n+2 
eae) if3<i<m 

I Bie) = (2) a pe 24) ‘ m+1 i an 
(" Sia io eran cule Say ifn+1l<icgm+n—-1 
(Vy) +tnt+mt+1 ifi=m+n 
m+n+2 ifi=m+n41 


Hence 


D(Bin) = 27? +221 ta)? + e™ 1 ta)? +27 (1 t+a)™. 


Corollary 2.22 The domination polynomial of the bi-star graph Bim) 1s 


D(Bian), 2) = (el +2)" +2"). 


Theorem 2.23 For the bi-star graph Bum), m An we have the following : 


3. if both m and n are odd, 
d(Bomm)) = 4 5 if both m and n are even, 
4 ifm and n have opposite parity. 


Proof By Theorem 2.21 we have, 


D Beavis ae ha ayers CL ea et CL 
=a? (amr 2 4 (Lt ayer pam ta) +ar-1(14+2)™) 
=a (eee) ea) (Cea 2es8)) 
=a? ((l+2)™+2™1) ((1+a2)"4+2"71). 


62 Shyama M.P. and Anil Kumar V. 


We have known that there is no real number satisfying both the equations (1+a2)™+a™~1 = 
0 and (1 +2)" +2"! = 0 simultaneously. Therefore it suffices to prove that (1+ a2)" +a™7! 


has exactly one real root for odd m and two real roots for even m. The remaining proof is 


similar to the proof of Theorem 2.9. 


Theorem 2.24 For bi-star graph Bryn), we have the following : 


2 ifn is odd, 
d( Brin) = 


3 ifn is even. 


Proof The proof similar to the proof of Theorem 2.9. 


The corona H o G of two graphs H and G is the graph formed from one copy of H and 
|V(H)| copies of G, where the i*” vertex of H is adjacent to every vertex in the i'” copy of G. 


Lemma 2.25([9]) Let G and H be nonempty graphs of order m and n respectively. Then 


D(Go H,2x) = (x(1+2)"+ D(H,z))”. 


Theorem 2.26 If Ky, and Ky, be the complete graphs with m and n vertices respectively. Then 


the domination polynomial of Ky, 0° Ky is 


m 


D(Km 0° Kn, x) = ((1+2)"t! —1) 


Theorem 2.27 For the corona Ky,° Ky, we have the following : 


2 ifn is odd, 
d(Kmo Kn) = f 
1 ifn is even. 


Proof It follows from the transformation y = 1+. in the domination polynomial D(K,, © 
Ky, 2). 


Consider the graph K,, and m copies of K,, . The graph Q(m, n) is obtained by identifying 


each vertex of K,, with a vertex of a unique Ky. 


Corollary 2.28 For m > 2, the domination polynomial of Q(m,n) is 


D(Q(m,n), x) = ((1+a)"—1)™. 


Proof It follows from the fact that Q(m,n) and Ky, 0 Ky_-1 are isomorphic. 
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Corollary 2.29 For the graph Q(m,n), we have the following : 


1 ifn is odd, 


2 ifn is even. 


§3. Domination Stable Graph 


In this section we introduce d-stable and d-unstable graphs. We obtained some examples of 
d-stable and d-unstable graphs. 


Definition 3.1 Let G = (V(G), E(G)) be a graph. The graph G is said to be a domination 
stable graph or simply d-stable graph if all the nonzero domination roots of G lie in the left open 
half-plane, that is, if real part of the nonzero domination roots is negative. If G is not d-stable 


graph, then G is said to be a domination unstable graph or simply d-unstable graph. 
Theorem 3.1 IfG and H are isomorphic graphs, then G is d-stable if and only if H is d-stable. 


Proof It follows from the fact that if G and H are isomorphic graphs then D(G,xz) = 
D(H, x). 


Corollary 3.2 If G and H are isomorphic graphs then G is d-unstable if and only if H is 
d-unstable. 


Theorem 3.3 If a graph G consists of m components G1,G2,...,Gm, then G is d-stable if 
and only if each G; is d-stable. 


Proof It follows from the fact that 


m 


D(G,x) = |] D(Gi, 2). 


i=l 


Corollary 3.4 If a graph G consists of m components G1,G2,...,Gm, then G is d-unstable 
if and only if one of the G; is d-unstable. 


Theorem 3.5 Let G be a connected graph of order n > 2 without pendent vertices. If G is 
d-stable, then 
n<1+2d(G,n- 3). 


Proof Suppose G is d-stable. Then by Routh-Hurwitz criteria, we have Routh-Hurwitz 
matrix Hj > 0. This implies that 


a(G,n — 1)d(G,n — 3) — d(G,n — 2) > 0. 
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Since G is connected and without pendent vertices we have 


ad(G,n—1) =n and d(G,n—- 2) = sn(n— 1). 


This completes the proof. 


Theorem 3.6 The complete graph K,, is d-stable graph for all n. 
Proof The domination polynomial of Kk, is 
D(Ky, 2) = (1+ 2)" - 1. 


Therefore 


Z(D(Kn,2)) = {exp (=) —1Jk=0,1,....n— i} 


Clearly, real part of all the roots are non-positive. This implies that K,, is d-stable for all n. 


Theorem 3.7 The complement of the complete graph K,, is d-stable graph for all n. 


Proof It follows from the fact that the graph K,, has no nonzero domination roots. 


We use the following definitions and results to prove some graphs which are d-unstable. 
These definitions and theorems are taken from [10]. 


Definition 3.2 If f,(x) is a family of complex polynomials, we say that a number z € C is a 
limit of roots of fr(x) if either fr(z) = 0 for all sufficiently large n or z is a limit point of the 
set Z(fn(x)), Z(fn(a)) is the set of the roots of the family fr(x). 


Now, a family f,(a) of polynomials is a recursive family of polynomials if f,,(a) satisfy a 


homogeneous linear recurrence 
k 
fa(t) = >) ai(2) fail), (8) 
i=1 


where the a;(a#) are fixed polynomials, with a,(a) # 0. The number k is the order of the 


recurrence. We can form from equation (8) its associated characteristic equation 
dM — ay (x)AP-1 — agdP-? — --- — ag (2) = 0 (9) 


whose roots A = X(x) are algebraic functions, and there are exactly k of them counting multi- 
plicity. 

If these roots, say A1(z), A2(x),--- , Ax (x), are distinct, then the general solution to equation 
(8) is known to be 


fn(a) = oy a(x) Ai (@)" (10) 
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with the usual variant if some of the A;(#) are repeated. The functions 
ai (2), a(x), mre , A (2) 


are determined from the initial conditions, that is, the k linear equations in the a; obtained 
by letting n = 0,1,--- ,&—1 in equation (10) or its variant. The details are available in [10]. 
Beraha, Kahane and Weiss [10] proved the following results on recursive families of polynomials 
and their roots. 


Theorem 3.8 If f,(x) is a recursive family of polynomials, then a complex number z is a limit 
of roots of fn(x) if and only if there is a sequence (Zn) in C such that fn(Zn) = 0 for all n and 


zn 7 Z2AS8 NO. 


Theorem 3.9 Under the non-degeneracy requirements that in equation (10) no a;(ax) is iden- 
tically zero and that for no pair i £ j is it true that A;(x) = wA;(x) for some complex number 


w of unit modulus, then z € C is a limit of roots of fr(x) if and only if either 


(1) two or more of the r;(z) are of equal modulus, and strictly greater (in modulus) than 


the others; or 


(2) for some j, A; (z) has modulus strictly greater than all the other r;(z), and a;(z) =0. 
Corollary 3.10((6]) Suppose fp(a) is a family of polynomials such that 
fr(a) = a1 (@)A1(@)” + a2(%)A2(x)" +... + an(x)AR(z)”, (11) 


where the a;(x) and the A;(x) are fixed non-zero polynomials, such that for no pair i # j is 
Ai (xz) = wdA;(x) for some w € C of unit modulus. Then the limits of roots of fn(x) are exactly 
those z satisfying (1) or (2) of Theorem 3.9. 


Theorem 3.11 The generalized barbell graph Bm n,1 18 d-stable for all m,n. 
Proof We have known by Theorem 2.17 that the domination polynomial of By n,1 is 
D(Bmnairt) = [1 +2)" — [i +2)" — I]. 


Therefore 


Qkri 
Z(D(Bmnir®)) = {exp ( =) a= 0,--.,m—1} 
m 


U{ ew (**) —1|k=0,--- aaah. 


Clearly, real part of all the roots are non-positive. This implies that the generalized barbell 


graph By in, is d-stable for all m,n. 


The domination roots of the generalized barbell graph By, 1 for 1 < m,n < 30 are shown 
in Figure 1. 
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Figure 1 Domination roots of By» for 1< m,n < 30. 


Corollary 3.12 The n-barbell graph By, is d-stable for all n. 


Proof It follows from the fact that the n-barbell graph B,,; and the generalized barbell 


graph B,,y,1 are isomorphic. 


The domination roots of the n-barbell graph 6, for 1 <n < 60 are shown in Figure 2. 


t+ 


Figure 2 Domination roots of B,,; for 1 <n < 60. 


Theorem 3.13 The corona Ky», 0 Ky, is d-stable for all m,n. 


Proof Notice that the domination polynomial of Ky, 0 Ky is 


m 


D KO hye) = ((1+a)"" =1) 


Therefore 


Z(D(Km 0 Kn,2)) = {exp (=) Thi ew nh. 
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Clearly, real part of all the roots are non-positive. This implies that the corona K,, o Ky is 
d-stable for all m,n. 


Corollary 3.14 The graph Q(m,n) is d-stable for all m,n. 


Proof It follows from the fact that the graph Q(m,n) and K,, 0 Ky -1 are isomorphic. 


Theorem 3.15 Let G be a connected graph of ordern and D(G, x) be its domination polynomial. 
If D(G,x) has exactly two distinct domination roots, then G is d-stable for all n. 


Proof It follows from the fact that the two distinct roots are real. 


Theorem 3.16 Let G be a graph of order n, then the corona Go Ky is d-stable for all n. 


Proof We have known the domination polynomial of Go Ky is 


D(Go ki,2) =a2"(x +2)”. 


Therefore Z(D(G o Ky, x) = {0,—2}, that is, Go Ky is d-stable for all n. 


Theorem 3.17 Let G be a graph of order n, then the corona Go Ky is d-stable for all n. 


Proof Notice that the domination polynomial of Go KQ is 


D(G 0 Ko, 2) = £3 (x? +3241)3. 


Therefore Z(D(G 0 Ko, 2) = {0, aval , That is, Go Kg is d-stable for all n. 


Theorem 3.18 Let G be a graph without pendent vertices and let D(G,x) be its domination 
polynomial. If D(G,x) has exactly three distinct roots, then G is d-stable. 


Proof Notice that 


Z(D(G,2)) c {0.24 iv2, a 


This implies that G is d-stable. 


Theorem 3.19 Any graph G with three distinct domination roots is d-stable 
Proof Notice that 


Z(D(G, 2) ¢ {20 2$8 2 if, 3* A, 


This implies that G is d-stable. 


Theorem 3.20 The Dutch windmill graph G3 is not d-stable graph for all but finite values of 
n. 
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Proof Using maple, we find that the Dutch windmill graph G¥ is d-stable for n < 6. Notice 
that D(G?,z) = z(14+ 2)?” + (22 +27)". We rewrite fn(x) = D(G%, x) as 


f(a) = 2 ((1+2)?)” + (1)(2e +02)” = aa dP + a2, 


where, a1 = 2, Ay = (1 +2)”, ag =1, A. = 2a +27. 

Clearly, 1 and x are not identically zero and Ay #4 wA2 for any complex number w of 
modulus 1. Therefore the initial conditions of Theorem .19 are satisfied. Now, for z= a-+ibe 
C, |Ai(z)| = |Ag(z)| holds if and only if |(1 + z)?| = |2z + 2?|. That is, |(1 + a + ib)?| = 
|2(a + ib) + (a + ib)?|. By a simple calculation we have (a + 1)? + 6? = $. Therefore 0 and the 

1 


complex numbers z such that (1+ R(z))? + (Z(z))? = § are limits of domination roots of GY. 


This implies that the domination roots of G3} have unbounded positive real part. Therefore the 


Dutch windmill graph G3 is not d-stable for all but finite values of n. 


The domination roots of the Dutch windmill graph G¥ for 1 <n < 6 and for 1 < n < 30 
are shown in Figures 3 and 4, respectively. 
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Figure 4 Domination roots of G3 for 1 <n < 30. 
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Figure 6 Domination roots of B, for 1 <n < 30. 
The domination roots of the bipartite cocktail party graph B, for 1 < n < 9 and for 


1 <n < 30 are shown in Figures 5 and 6, respectively. 


Remark 3.21 The domination polynomial of $;, is 


DS) a” +o(1 +o)” 
1(a)" +a(14+ 2)" 


ayAp + a2A5, 


l| 


l| 


where q3 = 1, \y = %, Q2 = & and Ag = 14+ a. Clearly 1 and & are not identically zero and 
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Ay # wg for any complex number w of modulus 1. Therefore the initial conditions of Theorem 
3.9 are satisfied. Now, |Ai| = |A2| holds if and only if |a — 0] = |a — (—1)|, that is, if and only if 
x is equidistant from 0 and —1. This holds if and only if real part of x is —4. Also a, is never 
0 and ag = 0 if and only if x = 0 and in this case |A2(0)| = 1 > 0 = |Ai(0)|. By these arguments 
we have 0 and the complex numbers z such that R(z) = —3 are the limits of roots of D(S;,, 2). 
Therefore we think that that there is no complex number z with positive real part is a root of 
D(S;,,x). We conjectured that the star graph S,, is d-stable graph for all n. 


The domination roots of the star graph S,, for 1 < n < 60 are shown in Figure 7. 


Figure 7 Domination roots of S, for 1 <n < 60. 


Remark 3.22 The domination polynomial of Ky, is 
D(Kmyn,t) = (1 +2)" —1)((l4+2)"-1)+a"+4+a2". 


Let m be fixed and rewrite D( Kinin, 2) as: 


I 


D(Kmn,t) =(1+a)™—1)(1+a)" + ((1+a™— (1 +a)™)) 1)” + 1(a)” 


aA} + agr5 + aga, 


where ay = (L4+a)"—1,A, =1+2,a9 =14+a™—-(14+2)™, 2 = 1, a3 = 1 and A3 = a. 
Clearly a1,a2 and a3 are not identically zero and A; 4 wA; for i A j and any complex number 
w of modulus 1. Therefore the initial conditions of Theorem 3.9 are satisfied. Now, applying 
part(i) of Theorem 3.9, we consider the following four different cases: 


4) |Aa| = |Aal = |Asl, 
it) |Aa| = |Aal| > |Asl, 
iit) |Aa| = |As] > |Aal, 


( 
( 
( 
(iv) [Aa] = [As] > [Ail 


Domination Stable Graphs 71 


Case 1. Assume that |1 + 2] = |1| = |a|. Then |a — (—1)| = |” — 0| implies that x lies on the 


vertical line z = —4, |x — (—1)| = 1 implies that 2 lies on the unit circle centered at (—1,0) 
and 1 = |a — 0| implies that x lies on the unit circle centered at the origin. Therefore the two 
points of intersection, s =e ig are limits of roots. 


Case 2. Assume that |1 +2] = |1] > |a|. Then |a — (—1)| = 1 implies that x lies on the unit 
circle centered at (—1,0), |e — (—1)| > |x — 0| implies that x lies to the right of the vertical 
line z = —4. Therefore the complex numbers x that satisfy |a — (—1)| = 1 and R(x) > —4 are 
limits of roots. 


Case 3. Assume that |1 + 2| = |2| > |1]. Then |x — (—1)| = |x — 0] implies that x lies on 


the vertical line z = —4 and |x — 0| > 1 implies that 2 lies outside the unit circle centered at 
the origin. Therefore the complex numbers « that satisfy |x| > 1 and R(x) > —3 are limits of 


roots. 


Case 4. Assume that |1| = |z| > |1+<2|. Then 1 = |x—0| implies that z lies on the unit circle 
centered at the origin and |a — 0| > |a — (—1)| implies that x lies to the left of the vertical line 
x = —5. Therefore the complex numbers x that satisfy |z| = 1 and R(x) < —$ are limits of 
roots. 


Also there may be some additional isolated limits of roots, being roots of a2 inside |1+2] = 1 
and |x| = 1. The union of the curves and points above yield that for m fixed, the limits of roots 
of the domination polynomial of the complete bipartite graph Ky, consists of the part of the 
circle |z| = 1 with real part at most -4, the part of the circle |z + 1| = 1 with real part at least 
—+ and the part of the line R(z) = —4 with modulus at least 1. So we conjectured that the 
complete bipartite graph Ky, is d-stable for all m,n. 


The domination roots of the complete bipartite graphs Ky, for 1 <m<15,1<n < 30 


and Ky,» for 1 <n < 30 are respectively shown in Figures 8 and 9. 


Figure 8 Domination roots of Ky, for 1<m<15 and 1<n< 30. 
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Figure 9 Domination roots of Ky, for 1 <n < 30. 


Remark 3.23 We have that D(Bimmn),v) = 2? ((l+2)™+a2™~1) (1+2)" +21). Let m 
be fixed, we rewrite D(Bomn), 2) as fn(a) : 


fa(z) = (e™1 4+-272(14+2)™) 1+)" 4+ (e™4+2(1+2)™) 2" 


= aA? + azAJ, 
where 
ay = (eth + 27?(1+2)"), A, =142, a2 =(a"+a(1+2)™) and A2 =z. 


Clearly (x™t! + 2?(1+a)™) and (x™ + 2(1+2)”) are not identically zero and A; # wAz for 
any complex number w of modulus 1. Therefore the initial conditions of Theorem 3.9 are 
satisfied. Now, |Ai| = |A2| holds if and only if |~ — (—1)| = |x — OJ, that is, if and only if x is 
equidistant from —1 and 0. The latter holds if and only if R(x) = —4. Notice that a1(0) = 0 
and a;(0) = 1+0=1 has modulus strictly greater than \2(0) = 0. 


Note that there may be some additional limits of roots, being roots of a, and a2. But from 
the Remark 3.21, we can conclude that a, and a2 have no roots in the right-half plane. By 
these arguments we have 0 and the complex numbers z that satisfy R(z) = —4 are the limits of 
roots of D(Bim,n),«). So we conjectured that the bi-star graph Bm) is d-stable for all m,n. 


The domination roots of the bi-star graph Bi,,,) for 1 <n < 50 are shown in Figure 10. 
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Figure 10 Domination roots of bi-star graph By, ,) for 1 <n < 50. 
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Abstract: Let Z, be the commutative ring of residue classes modulo n, PG(Z,) be the 
prime graph of a ring over a ring Z,. In this paper we study Energy and Wiener index of 
PG(Z,,) and give some results of line graph of prime graph of a ring over a ring Zn, denote 
it by L(PG(Zn)). 
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§1. Introduction 


Prime graph of a ring first introduced by Satyanarayana et al. [3]. Prime graph of a ring is 
defined as a graph whose vertices are all elements of the ring and any two distinct vertices 
x,y € Rare adjacent if and only if eRy = 0 or yRx = 0. This graph is denoted by PG(R). The 
concept of energy and Wiener index of zero divisor graph was introduced by Mohammad Reza 
and Reza Jahani in [4]. Motivated from the article in [4] in Section 2 of this paper we discuss 
energy of prime graph of a ring and give general MATLAB code for our calculation. In section 
3, We calculate Wiener index of PG(Z,), for n = p, n = p* and n = p®. In last section of 
paper, we introduce Line Graph of Prime Graph of a Ring denoted by L(PG(Z,,)) and discuss 
Planerity, Girth and degree of all vertices in L(PG(Z,,)). Also, we find center, eccentricity, 
point covering number, independence number, Energy, Wiener index and Chromatic number 
of L(PG(Z,,)), where n = p, p prime. Here, we also discuss complement of line graph of prime 
graph of a ring over a ring Z,,, denote it by L(PG(Z,,))°. We study Girth of L(PG(Z,,))*° and 
also find Eulerianity and degree of all vertices in L(PG(Z,,))°, where n = p, p prime. 
For more preliminary definitions and Notations the reader is referred to [5]-[8]. 


§2. Energy of Prime Graph of a Ring 


In this section we give some examples and calculate the Energy of prime graph of a ring. 


lReceived April 19, 2018, Accepted August 12, 2018. 
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Definition 2.1 The energy of the prime graph of a ring PG(Z,,) is defined as the sum of the 
absolute values of all the eigen values of its adjacency matrix M(PG[R]). ae. if A1,A2,°++ 5 An 
are n eigen values of M(PG[R]), then the energy of PG(Z,) is - 


E(PGIR)) = >i. 


Example 2.2 For p = 2, the adjacency matrix of PG(Z2) is 


1 
M(PG|Z3}) = 
0 
The characteristic polynomial is \? — 1.The eigen values are A; = 1,2 = —1. Therefore, 


E(PG[Zg]) = 2. 


Example 2.3 For p = 3, the adjancency matrix of PG(Zs) is 


0 1 1 
M(PG|Zs3]) = }1 0 0 
1 0 0 
The characteristic polynomial is 8 — 2. The eigen values are Ay = —1.4142,r5 = 


1.4142, 43 = 0. Therefore, E(PG[Zs3]) = 2.8284. 


Example 2.4 For p = 4, the adjancency matrix of PG(Z,) is 


M(PG{Za]) = 


Ee Fe FE oO 
ot 
Or (OS. Or 
Oo. (OOS oe 


The characteristic polynomial is \4 — 3\?. The eigen values are \y = 1.7321,\2 = 
—1.7321, A3 = 0,A4 = 0. Therefore, E(PG[Z,4]) = 3.4641. 


Example 2.5 For p=5, the adjancency matrix of PG(Zs) is 


M(PG[Z5]) = 


Ee PF EF EF oS 
(om > > 
(oe > 
SOS -O 2 
on > eo 


The characteristic polynomial is A° — 4A°. The eigen values are \y = 2,2 = —2,A3 = 
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0,A4 = 0, As = 0. Therefore, E(PG[Zs]) = 4. 


From the above Discussion we conclude the following theorem. 
Theorem 2.6 If p is a prime number then energy of PG(Z,) is 2\/p— T. 


General MATLAB code to find Energy of a Graph 


syms A To create Symbolic Variables; 

A=[--+3+++3+-+5---] To create a matrix that has multiple rows, separate the rows with semicolons; 
charpoly(A, A) Returns the characteristic polynomial of A in terms of variable ; 

p=| | To input the coefficients of characteristic polynomial; 

r = roots(p) Gives the eigen Values of matrix A; 


s = sum(abs(r)) Gives the energy of a graph. 


The values of E(PG[Z,,]) for n = 2,3,4,5,6,9 and 10 are given in table below. 


| Sr.No. n Characteristic Polynomial Energy | 


2 7-1 2 

3 r3 — 2X 2.8284 
4 At — 3? 3.4641 
5 
6 


pee 4 
TP AN Ay? 6.6858 
9 d® — 9rA7 — 208 + 645 7.4641 
10 «AM — 138 —8A7 +164 —- 9.2058 


§3. Wiener Index of Prime Graph of a Ring 


In this section, We calculate Wiener index of PG(Z,), for n = p, n = p? and n= p®. 


Definition 3.1 Let PG(R) be a Prime Graph of a Ring with vertex set V. We denote the length 
of the shortest path between every pair of vertices x,y € V with d(x,y). Then the Wiener index 
of PG(R) is the sum of the distances between all pair of vertices of PG(R), i.e. 


W(PG[R])= S© d(a,y). 


zyeV 


The following results can be easily verified. 
Theorem 3.2 W(PG[Z,]) = (p— 1)? if p is a prime. 


Theorem 3.3 W(PG[Z,2]) = ee) - (2p? — 2p +1] if p is a prime. 


Theorem 3.4 W(PG[Z,3]) = pip) [2p* + 2p? — 2p — 3] if p is a prime. 
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§4. Line Graph of Prime Graph of a Ring 


In this section we define line graph of prime graph of a ring, presented some examples and give 


some results. 


Definition 4.1 The line graph L(PG(Z,,)) of the graph PG(Z,) is defined to the graph whose 
set of vertices constitutes of the edges of PG(Z,), where two vertices are adjacent if the corre- 


sponding edges have a common vertex in PG(Z,,). 


Consider Z,, the ring of integers modulo n. 
Example 4.2 L(PG(Zz)) is a single vertex graph, there is no edge in L(PG(Z2)). 


Example 4.3 In L(PG(Z3)), there is an edge between the vertices [0,1] to [0,2], as shown in 


figure below. 


Figure 1 


Example 4.4 In L(PG(Z,)), there is an edge between the vertices [0,1] to [0,2], [0,2] to [0,3] 
and [0,3] to [0,1] as shown in figure below. 


Figure 2 
i.e. L(PG(Z4)) is a complete graph ks. 


Example 4.5 In L(PG(Zs)), there is an edge between the vertices [0,1] to [0,2], [0,2] to [0,3], 
[0,3] to [0,4], [0,4] to [0,1], [0,1] to [0,3] and [0,2] to [0,4] as shown in figure below. 


Figure 3 


i.e. L(PG(Zs)) is a complete graph kq. 
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Example 4.6 Let us construct L(PG(Ze)). 


Figure 4 


i.e. L(PG(Ze)) contains a complete subgraph ks. 


Example 4.7 Let us construct L(PG(Z7)). 


Figure 5 


i.e. L(PG(Zz)) is a complete graph keg. 
Observations 4.8 Every L(PG(Z,,)) contains a complete subgraph on n — 1 vertices. 
Observations 4.9 If Z, is a prime ring then L(PG(Z,)) is a regular graph. 


Observations 4.10 If m = p, a prime number then PG(Z,,) is a star graph. So, its line 
graph L(PG(Z,,)) is a complete graph and hence its eccentricity e(v) = 1,Vuv € V(L(PG(Z,))). 
Therefore, centre is L(PG(Z,)). 


Theorem 4.11 The graph L(PG(Z,)) is Hamiltonian if and only if n = p, a prime number 
andn> 4. 


Proof When n = 2, L(PG(Z,)) is a single vertex graph, hence there is no cycle. For 
n = 3, L(PG(Z,)) is a single edge graph, hence there is no cycle exist. For n = 4, L(PG(Z,)) 
is a triangle graph and there exist a cycle which containing every vertex. So, L(PG(Z,)) is a 
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Hamiltonian graph. Now, for n = p, a prime number then L(PG(Z,,)) is Hamiltonian graph 
because there exist a cycle containing every vertex. Hence, the graph L(PG(Z,,)) is Hamiltonian 


if and only if n = p, a prime number and n > 4. 


Theorem 4.12 Let L(PG(Z,,)) be a line graph of prime graph of a ring, where n = p and p 
is an odd prime number then point covering number and independence number of L(PG(Zn)) 


both are one. 


Proof When n = p, PG(Z,,) is a star graph. So, there is a common vertex which is 
adjacent to all other vertices and that vertex is called center of the graph. When we draw the 
line graph of PG(Z,,), for n = p, and let a; = 0 be the common vertex of PG(Z,,) which is 
the end point of every edge of PG(Z,,). Then a; appears in every vertex of the line graph. 
[a1, ui] € V(L(PG(Z,,))), where i = 1,2,3,--- ,(p — 1) forms a complete line graph of PG(Z,,) 
and here, [{a1, v1] is adjacent with all other vertices of line graph. In other words, we can say 
that single vertex cover all other vertices of line graph of PG(Z,,). Thus, the point cover is one 


and from that vertex an independence number is also one. 


The following results can be immediately verified. 


Theorem 4.13 The general formula for degree of vertex in L(PG(Z,)) is: 


deg|u, v] = gced(u, n) + ged(v, n) — 2, if u2 £0 and v? £0 
= gcd(u, n) + ged(v, n) — 3, if either u2 =0, v? =0 
= gcd(u, n) + gced(v, n) — 4, if u2 =0 and v? =0 


Theorem 4.14 L(PG(Z,,)) is planer if and only if n = 2,3,4,5 and is non-planer for n > 6. 


Theorem 4.15 The girth gr(L(PG(Z,))) = 3 if and only ifn > 4. If n = 2,3 then 
gr(L(PG(Z,,))) = oo. 


Theorem 4.16 The chromatic number x(L(PG(Z,))) =p—1 for p = 2,3,5,---. 
Theorem 4.17 The chromatic number x(L(PG(Zp»))) = p” — 1, p prime. 


Theorem 4.18 The energy E(L(PG(Z,))) = 2p —4, for p= 3,5,--- andn=4. 


Theorem 4.19 The Wiener index W(L(PG(Z,))) = P=) | for p=3,5,:-- andn=4. 


Theorem 4.20 The graph L(PG(Z,,))° is Eulerian if and only if n =p, a prime number and 
n> A. 


Proof When n = 2, there is no graph, as there is no edge between the vertices 0 and 1 in 
(PG(Zn))°. For n = 3, L(PG(Z,))°is a single vertex graph. For n = 4, L(PG(Z,,))° is triangle 
graph and every vertex is of even degree. Now, For n = p, a prime number, every vertex of 
L(PG(Z,,))° have even degree. Hence, the graph L(PG(Z,,))° is Eulerian if and only if n = p, 
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a prime number and n > 4. 


Theorem 4.21 The general formula for degree of vertex in L(PG(Z,))°, where n = p a prime 
number and n> 5 is: 


deglu, v] =n + d(n) —5 


Theorem 4.22 The girth gr(L(PG(Zn))°) = 3 if and only ifn > 4. If n = 2,3 then 
gr(L(PG(Z,,))°) = 0. 
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Abstract: A tree T contained in graph G is a Steiner tree with respect to W C V(G) if 
T is a tree of minimum order with W C V(T). The set S(W) consists of all the vertices 
of G which lie on some Steiner tree with respect to W. The set W is a Steiner set for G if 
S(W) = V(G). The minimum cardinality among the Steiner sets of G is the Steiner number 
of G, denoted as s(G). The set W is called Steiner dominating set if W is both a Steiner set 
and a dominating set. The minimum cardinality among such sets is a Steiner domination 
number, denoted as ys(G). We investigate Steiner domination number of some splitting and 
degree splitting graphs. 

Key Words: Steiner distance, Steiner set, Steiner number, domination number, Steiner 


domination number. 


AMS(2010): 05C69, 05C76. 


§1. Introduction 


We consider simple, finite, connected and undirected graph G with vertex set V and edge set 
E. For the standard graph theoretic terminology and notation we follow Chatrand and Lesniak 
[2] while the terms related to the theory of domination are used in the sense of Haynes et al. 
[6]. 


Definition 1.1 The distance d(u,v) between two vertices u and v in a connected graph G is 
the length of the shortest u—v path in G. 


Definition 1.2 The Steiner distance sd(W) of a subset W of vertices of a connected graph 
G is the minimum number of edges in a connected subgraph of G that contains W. If H is a 
subgraph of minimum size that contains a set W, then H is necessarily a tree, called a Steiner 


tree for W or a Steiner W-tree. 
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Chartrand et al. have introduced a generalization of distance in [3]. The sharp upper 
and lower bounds for the Steiner k-diameter of G and G are given by Mao [9] while the same 
author have identified some graph classes attaining these bounds. Let n be an integer such that 
2<n<|V(G)|, then the n diameter of G, diam,,(G), is defined to be the maximum Steiner 
distance of any n-subset(subset with n elements) of vertices of G. If G be any graph of order p 
with minimum degree 6 > 2 and 2<n<p then diam,(G) < a +2n—5, is proved by Ali 
et al. [1]. 


Definition 1.3 The set of all vertices of G that lie on some Steiner W -tree is denoted by S(W). 
If S(W) =V(G), then W is called a Steiner set for G. A Steiner set of minimum cardinality 


is a minimum Steiner set and this cardinality is the Steiner number s(G). 


The concept of Steiner number was introduced by Chartrand and Zhang [4]. In the same 
paper authors have proved many results on this newly defined concept. This concept was further 
studied by Santhakumaran and John [8]. For the graph G of Figure 1, there are three Steiner 
trees related to W = {w1, we} which are shown in the same figure. Since S(W) = V(G), W is 
a Steiner set of G. 

W 


W W WwW, 


i, 


W), W 
Figure 1 The graph G and its Steiner trees 


Definition 1.4 A set S C V of vertices in a graph G = (V,E) is called a dominating set if 
every verter vu € V is either an element of S or is adjacent to an element of S. A dominating set 
S is a minimal dominating set if no proper subset S’ C S is a dominating set. The domination 


number ¥(G) of a graph G is the minimum cardinality of a dominating set in graph G. 


Definition 1.5 Let G be a connected graph with vertex set V(G). A set of vertices W in G is 
called a Steiner dominating set if W is both a Steiner set and a dominating set. The minimum 


cardinality of a Steiner dominating set of G is called its Steiner domination number, denoted 
by ys(G). 


The concept of Steiner domination number was introduced by John et al. [7]. It is very in- 
teresting to investigate Steiner domination number of graph or graph families as it is known only 
for handful number of graphs. Vaidya and Mehta [11] have investigated the Steiner domination 
number of W,,, H,, and Fl, and the same authors [12] have established some characterizations 
for Steiner domination in graphs while Steiner domination number for S’(P,,), S’(Cn), M(Pr), 
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M(C,,) and F;, are obtained by Vaidya and Karkar [10]. 
For the graph G of Figure 1, W = {wi,w2} is a Steiner dominating set of minimum 
cardinality. Therefore, y;(G) = 2. 


Definition 1.6 A vertex v is an extreme vertex of a graph G if the subgraph induced by neighbors 
of v is complete. 


Definition 1.7([5]) A systematic visit of each vertex of a tree is called a tree traversal. 


Definition 1.8 The bistar By» is the graph obtained by joining the center(apex) vertices of 
Kim and Ky» by an edge. 


Definition 1.9 Let G be a graph with V(G) = S,; US_U $3 U---S,;UT where each S; is a set 


t 

of all vertices of the same degree with at least two elements and T =V(G) \ LU S;. The degree 
i=1 

splitting of G denoted by DS(G) is obtained from G by adding vertices w1, w2,W3,+++ , we and 

joining w; to each vertex of S; for1<i<t. Note that if V(G) = u S; then T = 0. 


Definition 1.10 For a graph G the splitting graph S'(G) of a graph G is obtained by adding a 
new vertex v' corresponding to each vertex v of G such that N(v) = N(v’). 


Definition 1.11 A friendship graph F,, is a one point union of n copies of cycle C3. 


§2. Main Results 


Observation 2.1 7(Bmn) = m+n. 


Theorem 2.2 y¥5(S’(Bmn)) =m+n+42. 


Proof Let u, ur, U2,+++ Um, VU, U1, V2,°°* Un bem+n+2 vertices of Bm» and u’, ul, us,--- ul 


m? 

vu’, U{,U5,-°+,v), be the corresponding vertices which are added to obtain S’(Bm»). Then 
Ul / / i. / / / / / / / 
vis (Bm,.n)) = {u, U1, U2, srt 5 Um, U1, U2,° °° Un, U,U Uy, Ug, °° Um, U 5 U1, V25°° Un}. Now wu}, Ug, 


+ Us, U4, U5,°°°U;, are extreme vertices as the subgraph induced by their neighbors is com- 
plete, namely, the complete graph K,. Therefore, they must be in Steiner dominating set 
W. If uj,us,---u, ul, ud,---u,, © W then uj,us,--- ul, up, vs,---u,,u,v € S(W). Now 
there some trees between u’ and v’ which include remaining vertices uy, U2,+ ++ Um, U1, V2,°°*Un- 
So if u’,ui,us,--- uy, U',U{,U5,---u,, € W then there are four Steiner W-trees which in- 
clude all the vertices of the graph. That is, if u’,uj,us,---ul,,v’,u},05,-°-u,, © W then 

/ / if / / / / / / / / 
U1, U2,° ++ Um; V1, V2,°° Un, Ww, Uy, Ug, UU’, Uy, Us,--- ul, © S(W). Therefore, W = {u’, ui}, us, 
--ul_,v',v1,U5,°°-U;,} becomes a Steiner set of minimum cardinality m+ n+ 2 and it is also 


a dominating set. Hence 


s(S"(Bmn)) =m+n+2. 


Theorem 2.3 y5(DS(Bmjn)) = 2. 


Proof Let u, uy, U2,++* Um, V, U1, 02,°°*Un be M+n+ 2 vertices of By» and 71, 2 be the 
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corresponding vertices which are added in order to obtain DS(By,»). Then, V(DS(Bm.n)) = 
{u, U1, U2,°+* Um, VU; U1, V2,°°*Un, 21,22}. Now if G is a connected graph of order n > 2 then 
2< S(G) <n. Without loss of generality let 71,72 € W then there are four Steiner W-tree 
traversal between x; and x2 which together include all the vertices of DS(Bm,»). Therefore, 
W = {21,22} becomes a Steiner set of minimum cardinality and it is also a dominating set. 
Therefore, W = {21,22} becomes a Steiner dominating set of minimum cardinality. Hence 


Ys(DS(Bmn)) = 2. 


Lemma 2.4 S(DS(P,)) =n—5, n> 7. 


Proof Consider P, with V(P,) = {v1,v2,---Un} with partition Vi = {v2,v3,--+Un—1} 
and V2 = {v1, Un}. Now in order to obtain DS(P,,) from P, we add x; and x2 correspond- 
ing to Vj and Va. Thus, V(DS(P,)) = {21,%2,01,V2,°++Un}. Let 21,04 € W then there 
are some Steiner W-trees which include the vertices x1, v1, v2, v3, U4, £2. So, if 71,v4 © W then 
L1, U1, V2, U3, U4, LQ © S(W). Let x1, v4, Vn—3 € W then 21, v1, v2, 3, U4, £2, Un—3; Un—2; Un—1; Un E 
S(W). Then, there does not exist tree traversal containing £1, V4, Un—3 which includes vs, v6, --- , 
Un—a. The vertices v5, Ug,°+* ,Un—a must be included in W to obtain Steiner tree of minimum 
size which include vs, v6,--- ,Un—a. Therefore, if 1, v4, U5,°++Un—4,Un—-3 € W. Then there are 


following four Steiner W- trees as listed below: 


(1) 
(2) £101 V2%2U4U5U6 * + Un—3, 
(3) 


(4) L1 Un Un—1L2Un—2Un—3Un—4° °° U5, U4, 


1 U1 U2U3 °° * Un—4Un-3; 


L1UnVn—1Un—2Un—3 °° * U5V4, 


which include all the vertices of the graph. Thus W = {21, v4, Us,-++Un—4a, Un—3} becomes a 
Steiner set of minimum size which include n — 6 vertices of P, and a vertex x;. Hence 


S(DS(P,)) =n—5. 


Theorem 2.5 y;(DS(P,)) =n-—3, n> 7. 


Proof From the Theorem 2.4 W = {a1,v4,U5,°-+-Un—4,Un—3} is a Steiner set of mini- 
mum cardinality. But it is not a dominating set as vg and vy_1 are not dominated by any 
of the vertices. Therefore, these two vertices must be in Steiner dominating set W. So, 


{21, V2, U4, U5,°°* Un—4, Un—3, Un—1} is a Steiner dominating set of minimum cardinality. Hence 


9y5(DS(P,)) =n—3. 


Proposition 2.6([7]) Ys(Kmn) = min{m,n} if m,n > 2. 
Theorem 2.7 ¥5(S'(Kmn)) =m+n. 


Proof Let v1, v2,°++Um,U1,U2,°+* Un be m+n vertices of Km mn. Now v1, v$,-+-v1,, uy, ud, 


- ,ul, be the corresponding vertices which are added in order to obtain S’(Kym,,) with parti- 
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tions W = {v},v},---ul,,uy,ug,--- ul} and X = {v1, v2,-++ Um, U1, U2,°+* Un}. It is very clear 
that W is a Steiner set as there are max{m,n} number of Steiner trees which include all the 
vertices of the graph. Here W dominates all the vertices of the graph. Therefore, it is also a 
dominating set. Thus, W is a Steiner dominating set. We claim that W is a Steiner dominating 
set with minimum cardinality. If possible let U be any Steiner set such that |U| < |W| and 
U CW. Then, there exists a vertex vi € W such that vj ¢ U. But as the vertices of W are 
mutually non adjacent, the Steiner U-tree containing vj and vu, (7 Ai, k Ai, 1<7,k <n) will 
not contain v;. Therefore, U is not Steiner set. If VU C X then some vertices of W and some 
vertices of X which are not included in U are not in any Steiner U- trees. Therefore, U is not 
Steiner set. Let U CWUX such that U contain at least one vertex from each of W and X 
then some vertices of W and X do not lie on any Steiner U-tree. Thus, U is not a Steiner set. 


So, W is a Steiner dominating set of minimum cardinality m+n. Hence 


eS Kea) =m+n. 


Proposition 2.8([4]) Let G be a connected graph of order p> 2. Then y,(G) = 2 if and only 
if there exists a Steiner dominating set S = {u,v} of G such that d(u,v) < 3. 


Theorem 2.9 y5(DS(Kmn)) =2, mAn, m,n > 2. 


Proof Let v1, v2,+++Um,Ui1,U2,***Un be m+n vertices of Km» with partitions W = 
{v1, V2,°++Um} and X = {uj,u2,--+Un}. In order to construct DS(Kinn) we add wy; and we. 
If we consider w, and wz in Steiner set W then S(W) = V(G) and W is also a dominating set. 
Therefore W becomes a Steiner dominating set and d(w1, w2) = 3. Hence by Proposition 2.8, 


Ys(DS(Kmn)) = 2. 


Proposition 2.10([7]) Each extreme vertex of a connected graph G belongs to every Steiner 


dominating set of G. 
Theorem 2.11 y5(S'(Fn)) = 2n+ 1. 


Proof Let vo, v1,V2,---Un; Un+1;°°*V2n be the 2n + 1 vertices of F;, where vp is the apex 
vertex. Now v9, U1, U9,--- Uy, Un 415° ** Von, be the vertices which are added to obtain S’(F;,). The 
vertices Up, Vj, Ug, ---U,,Un41,°** Von must be in Steiner dominating set W as they are extreme 
vertices. But W = {v9, vj, Vg, --- U2, Un4is° ++ Von} is not a Steiner dominating set as it is neither 
a Steiner set nor a dominating set. Therefore, we must include some more vertices to obtain a 
Steiner dominating set. Let vp € W then S(W) = V(S’(F;,)) and 

Ul 


ON GTS BP / / 
W= {5 U1, V2," "Uns Un41)"* "Von } 


is a dominating set of minimum cardinality. Hence 


s(8"(Fh)) = 2n +1. 
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Concluding Remarks 


The Steiner domination in graphs is one of the interesting domination models. It is always 


challenging to investigate Steiner domination number of a graph. We have obtained Steiner 


domination number of larger graphs which are obtained by means of various graph operations. 
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Abstract: Coloring the vertices of a graph G according to certain conditions can be con- 
sidered as a random experiment and a discrete random variable X can be defined as the 
number of vertices having a particular color in the proper coloring of G. In this paper, we 
extend the concepts of mean and variance, two important statistical measures, to the theory 
of graph coloring and determine the values of these parameters for a number of standard 


graphs. 


Key Words: Graph coloring, Smarandachely A-coloring, coloring sum of graphs, coloring 


mean, coloring variance, y-chromatic mean, yt-chromatic. 


AMS(2010): 05C15, 62A01. 


§1. Introduction 


Investigations on graph coloring problems have attracted wide interest among researchers since 
its introduction in the second half of the nineteenth century. The vertex coloring or simply a 
coloring of a graph is an assignment of colors or labels to the vertices of a graph subject to 
certain conditions. For example, Smarandachely A-coloring of graph G by colors in @ such that 
yp(u) £ y(v) if w and v are elements of a subgraph isomorphic to graph A in G. In a proper 
coloring of a graph, its vertices are colored in such a way that no two adjacent vertices in that 
graph have the same color. 

Different types of graph colorings have been introduced during several subsequent studies. 
Many practical and real life situations paved paths to different graph coloring problems. 

Several researchers have also introduced various parameters related to different types of 
graph coloring and studied their properties extensively. The first and the most important 
parameter in the theory of graph coloring is the chromatic number of graphs which is defined 
as the minimum number of colors required in a proper coloring of the given graph. The concept 
of chromatic number has been extended to almost all types of graph colorings. 

The notion of chromatic sums of graphs related to various graph colorings have been 
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introduced and studied extensively. Some of these studies can be found in [9, 10, 11]. The 


notion of a general coloring sum of a graph has been explained in [9] as follows: 


Let C = {c1,¢2,c3,--: ,c.} be a particular type of proper k-coloring of a given graph 
G and @(c;) denotes the number of times a particular color c; is assigned to vertices of G. 
Then, the coloring sum of a coloring C of a given graph G, denoted by we(G), is defined to be 
k 
we(G) = >> i O(c). 
i=1 
Motivated by the studies on different types of graph coloring problems, corresponding 
parameters and their applications, we discuss the concepts of mean and variance, two important 


statistical parameters, to the theory of graph coloring in this paper. 


For all terms and definitions, not defined specifically in this paper, we refer to [2, 3, 4, 
6, 15, 16] and for the terminology of graph coloring, we refer to [5, 7, 8]. For the concepts in 
Statistics, please see [12, 13]. Unless mentioned otherwise, all graphs considered in this paper 
are simple, finite, connected and non-trivial. 


§2. Coloring Mean and Variance of Graphs 


We can identify the coloring of the vertices of a given graph G with a random experiment. Let 
C = {c1,¢2,¢3,:-: ,ce} be a proper k-coloring of G and let X be the random variable (r.v) 
which denotes the color of an arbitrarily chosen vertex in G. Since the sum of all weights of 
colors of G is the order of G, the real valued function f(i) defined by 


ee a 
f(i) = |V(G)|? tw 1, 2,3, k 
0; elsewhere 


is the probability mass function (p.m.f)of the r.v X. If the context is clear, we can also say 
that f(i) is the p.m.f of the graph G with respect to the given coloring C. 


Hence, analogous to the definitions of the mean and variance of random variables, the mean 


and variance of a graph G, with respect to a general coloring of G can be defined as follows. 


Definition 2.1 Let C = {c1,c2,¢3,:-- ,ch} be a certain type of proper k-coloring of a given 
graph G and 0(c;) denotes the number of times a particular color c;, is assigned to vertices of 
G. Then, the coloring mean of a coloring C of a given graph G, denoted by ~c(G), is defined 
to be 


Definition 2.2 For a positive integer r, the r-th moment of the coloring C is denoted by tucr (G) 
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and is defined as 


> ae O(c; ) 
==, 
Definition 2.3 The coloring variance of a coloring C of a given graph G, denoted by o2(G), 
is defined to be 


2.1 yx-Chromatic Mean and Variance of Graphs 


Coloring mean and variance corresponding to a particular type of minimal proper coloring of 
the vertices of G are defined as follows. 


Definition 2.4 A coloring mean of a graph G, with respect to a proper coloring C is said to be 
a x-chromatic mean of G, if C is the minimum proper coloring of G and the coloring sum we 


is also minimum. The x-chromatic mean of a graph G is denoted by py. 


Definition 2.5 The y-chromatic variance of G, denoted by 0%(G), is a coloring variance of G 


with respect to a minimal proper coloring C of G which yields the minimum coloring sum. 


Let us now determine the x-chromatic mean and variance of certain standard graph classes. 


The following result discusses the y-chromatic mean and variance of a complete graph Ky. 


Proposition 2.6 The x-chromatic mean of a complete graph Ky, is net and its x-chromatic 
n?—1 


Variance 1S i) 


Proof Note that all vertices of a complete graph K,, must have different colors as they are 
all adjacent to each other. That is, 6(c;) = 1 for color c;, 1 <i<n. Therefore, 


and 


6 2 12 


ae m+1\?_ (n+1)Qn+1) (n+1)? _ n?-1 
(Ka) ==) #- ( 2 ) = : 
i=1 
The following theorem gives the probability distribution of a proper coloring of a complete 
graph. 
Theorem 2.7 Any proper coloring of a complete graph Ky, has the discrete uniform distribution 
on {1,2,--- ,k}(DU(k)). 


Proof Let X be the r.v representing the number of colors in a proper k-coloring of a 
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complete graph K,. For any proper k-coloring C of the complete graph K,,, 0(c;) = 1 and 
k =n. Hence, the corresponding p.m.f is 


; i; RH Wy2, 34545, 
fO= 


0; elsewhere 


which is that of the discrete uniform distribution on {1,2,--- ,k}. Hence, X ~ DU(k). 


The following result determines the y-chromatic mean and variance for a path P,. 
Proposition 2.8 The y-chromatic mean of a path P,, is 


; if n is even, 
Lx (Pn) = : = : : 
t—<; if n is odd, 


and the x-chromatic variance of Pp, is 


if n is even 
o% (Pn) = ae ee Os , 
2—=; ifn is odd. 


Proof Consider a path P, on n vertices. Being a bipartite graph, the vertices of P,, can 
be colored using two colors, say c; and cp. Then, we have the following cases. 


(i) If n is even, exactly 7 vertices of P, have color c, and 5 vertices have color cg. Then, 


the p.m.f of the corresponding r.v X is 


i= 1,2, 


; elsewhere. 


and the y-chromatic variance is 


(P=? h—-wyr=8-(3) =f. 


i=l 


(it) If n is odd, then the p.m.f of the corresponding r.v X is 


f=); i=2, 


0; elsewhere. 
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Then, the y-chromatic mean of P,, is 


n+1 n—-1 3n—-1 
P,) =1- . = 
Hy ( ) 2n Ea 2n 2n 


and its x-chromatic variance is 


2 2 
Py — 1 . 
ox ) 2n 


The following result determines the values of these parameters for a cycle Cy. 


Proposition 2.9 The yx-chromatic mean of a cycle Cy is 


mir a if n is even, 
3nt+3- ifn is odd, 


2n ? 


and the x-chromatic variance of Cy is 


+ if n is even, 


2 ‘ 2 
at if n is odd. 


Proof Consider a cycle C, on n vertices. Then, we have the following cases. 

(i) If n is even, then C;, is bipartite and is 2-colorable. Then, exactly } vertices of C;, also 
have color c; and cz each. Then, as explained in the first part of previous theorem, we have 
by (Cn) = 3 and 02 (Cn) = F. 

(it) If n is odd, then C;, is 3-colorable. Let C = {c1, c2,c3} be the minimal proper coloring 


of C,. Then, the p.m.f of the r.v X is given by 


a i = 1,2, 


2n ? 
fO™=44; if 1=3, 
0; elsewhere. 


Then, the y-chromatic mean of G is 


n—-1 n—-1 1 3n+3 
Gy) le 2- 3-- = 
H( ) 2n ¥ 2n - 2n 


and the y-chromatic variance of C’, is 


2 2 
9 _ f42 no-1 n—1 2 1\ /38n+3\"_ ni —8n+9 
o(Cn) ~ (1 2n % 2n +3 x) ( 2n ‘ An? , 


97. 


In the following theorem, we determine the x-chromatic mean and variance of a wheel 


graph W,, = Ki + Cy-1. 
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Proposition 2.10 The x-chromatic mean of a wheel graph W,, is 


Snts if n is odd, 


Ly(Wr) = an+1 ; ; 
sta; «| Cf n is even, 
and the x-chromatic variance of W,, is 

2 . . 

a if n is odd, 
oy (Wn) = ; 

n~+32n—64, - = 
a eT if n is even. 


Proof Note that the wheel graph W,, is 3-colorable, when n is odd and 4-colorable when 
is even. Then, we have the following cases. 


(i) First, assume that n is an odd integer. Then, the outer cycle C,-1 of W, is an even 
n—1 


cycle. Hence, 4 vertices of Cp—1 have color ¢;, not vertices of C,—, have color cg and the 
central vertex of W,, has color c3. Hence the corresponding p.m.f for W,, is given by 


ml. if ¢=1,2, 
Fiat ap pes, 


n? 


0; elsewhere. 


Hence, the corresponding y-chromatic mean is 


n—-1 n—-1 1 38n+3 


n=l: : : : 
Hx (Wn) 2n 5 2n n 2n 


Now, the x-chromatic variance is 


02 Wa) = 242) 2 Naat Gusting)? = (SAD 4 2) (Set 8) tend 


(ii) Next, assume that n is an even integer. Then, the outer cycle C,,-; of W,, is an odd 
cycle. Hence, noe vertices of the outer cycle C;,-, have color cy, nee vertices of C,,_1 have 


color cp and one vertex of C,_; has color c3 and the central vertex of W,, has the cy. Hence, 
the p.m.f for W,, is given by 


oe. ae pe 1 


f@=<2, if 1=3,4 


n? 


0; elsewhere. 


Hence, the corresponding y-chromatic mean is 


—2 —2 1 1 3 8 
i: pi oe Bie yore, eee nn 
2n 2n n n 2n 
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and the y-chromatic variance is 


—2 1 
o2(W,) = (12427). * 4 (374.4). = - (a (W,))? 
2n n 
_ (3m=2) P+" 3n+8\* _ n?+32n—64 
— Qn n 2n a 4n?2 . 


Remark 2.1 From the above discussion, we observe that the minimum proper coloring of 
bipartite graph follows a two-point distribution. In general, for a bipartite graph G(Vi, Vo, £), 
with |Vi| = m1 > |V2| = m2,m1 + m2 =n, the p.m.f can be defined as 

a; if i=1, 
f=); if 1-2, 


0; elsewhere. 


Hence, we have p,(G) = ™#m™ = 14 ™ and o2(G) = matting _ (1+ m2)? - 
— [(n — 1)m, + 2(2n — 1)mzy]. 


n 
2 
each partition and hence with respect to a minimal proper coloring, exactly 4 vertices having 


Remark 2.2 If G is a regular bipartite graph on n vertices, then there will be 4 vertices in 


the colors cy and cg each. Hence the p.m.f is 


saya g® ihe 
0; 


elsewhere. 


Hence, py(G) = 3 and o2%(G) = | as mentioned in Proposition 2.9. 


2.2 y+t-Chromatic Mean and Variance of Graphs 


Coloring mean and variance corresponding to another type of a minimal proper coloring of the 
vertices of G are defined as follows. 


Definition 2.11 A coloring mean of a graph G, with respect to a proper coloring C is said to be 
a x*-chromatic mean of G, ifC is a minimum proper coloring of G such that the corresponding 


coloring sum wg is maximum. The x+-chromatic number of a graph G is denoted by j1.+(G). 


Definition 2.12 The yt-chromatic variance of G, denoted by o%4(G), is a coloring variance 
of G with respect to a minimal proper coloring C of G such that the corresponding coloring sum 


1s Maximum. 


Invoking the definitions of two types of chromatic means and variances mentioned above, 
we can infer the following. 


Remark 2.3 For any arbitrary minimal proper coloring C of a given graph G, we have ju, (G) < 
lic(G) < py+(G) and o3(G) < 08(G) < o%,(G). 
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Remark 2.4 Since all vertices of a complete graph have different colors, the y-chromatic mean 
and the y*-chromatic mean are equal and the y-chromatic variance and the yt-chromatic 


variance are equal. 


Let us now discuss the y*+-chromatic mean and variance of the graph classes mentioned in 


the previous section. 


Proposition 2.13 The x*-chromatic mean of a path P,, is 


3; if n is even, 


Leas snl. ifn is odd, 


and the y*-chromatic variance of Py, is 


if n is even, 


xt mos - 
“Tne ifn is odd. 


Proof As in Proposition 2.8, we consider the following cases. 


(i) If n is even, as mentioned in Proposition 2.8, exactly $ vertices of P, have color c; and 


5 vertices have color cp. Then, the p.m.f of the corresponding r.v X is also as defined there. 
Hence, the x+-chromatic mean is j1,+(P,) = = and the y*+-chromatic variance is Car (2) ee 


(ii) If n is odd, x*-coloring assigns color c, to not vertices and color cg to the remaining 


net vertices. Then the p.m.f is 


n—1., i=1, 


2n? 
f=) SB; 1=2, 
0; elsewhere. 


Then, the y*-chromatic mean of P, is given by 


n—-1 n+l 38n+1 
+2: = 


P,) =1- 
H+ ( ) 2n 2n 2n 


and its y*-chromatic variance is given by 


n—-1 g nti Bn4+1\*  n?41 
3 +2°. _ = —. 
2n 2n 2n 4n? 


The following proposition discusses the y*-chromatic mean and variance of a cycle on n 


vertices. 
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Proposition 2.14 The x*-chromatic mean of a cycle Cy is 


nO = if n is even, 
a: sn=3. if n is odd, 


2n ? 


and the y*-chromatic variance of Py, is 


1 z : 
ai afm ts even 
oe (Cp) = 4? if 9 


2 . . 
ats. if n is odd. 


Proof Here, we have to consider the following two cases. 


(i) If n is even, as mentioned in Proposition 2.13, exactly $ vertices of C;, have color cy 


and color cz each. Then, exactly as explained there, we have, p.+(C,) = 2 and Os (Ca= = 
(ii) If n is odd, y*-coloring assigns color c; to one vertex, color cz to not vertices and 
color cg to the remaining not vertices of the cycle C,,. Then the p.m.f is 


1; =A, 
fi =4 4; 1=2,3 
0; elsewhere. 
Then, the y*+-chromatic mean of C,, is 
1 n—-1 n—-1 5n— 3 
n)b= _— 2 . — 
Beene Sop ce ne re oe, Qn 


and its x+-chromatic variance is 
xX 


The following proposition discusses the y+-chromatic mean and variance of a wheel graph 
on n vertices. 


Proposition 2.15 The x*-chromatic mean of a wheel graph W,, is 


5n—3 : ; 
; ifn is odd. 
on”? 7 
Hyt (Wr) = as ee, 
Qn--2? 


if n is even, 


and the x*-chromatic variance of Wy, is 


n-+30n—31. ‘i : 
aw ae; «= (ifn is odd, 
KE n2432n—64, 


Tz >~—COif n 4s even. 
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Proof As mentioned in Proposition 1.10, the wheel graph W,, is 3-colorable, when n is odd 


and 4-colorable when is even. Then, we have to consider the following cases. 
(i) First, assume that n is an odd integer. Then, the outer cycle C,_-1 of W,, is an even 
n—1 


cycle. Hence, we can assign color c, to the central vertex of W,, color cz to “= vertices of 
Cy~-1 and color cz to the remaining net vertices of C,-1. Hence the corresponding p.m.f for 


W,, is given by 


4; if «=1, 
f@)=% 5; if 1=2,3, 
0; elsewhere. 
Hence, the yt+-chromatic mean is 
1 n-1 n-1 5n-3 
pn ee eee ae eae 
noe (Wn) = 1 n : 2n ee 2n 2n 
and the y+-chromatic variance is 
2 2 1 2 2) noi 2 
Oy+(Wn) = W-— + (2° 43°). — (Heni(Wn)) 
_ (13(n-1) 1 5n—3\*  n?+30n—31 
a 2n n 2n - 4n? : 


(ii) Let n be an even integer. Then, the outer cycle C,1 of W,, is an odd cycle. Hence, we 
can assign color c; to the central vertex of W,,, color cz to one vertex of the outer cycle C,_1, 


color cz to 2=* vertices of C;,_1 and color cy to the remaining 24 vertices of C,,_1. Therefore 
3 2 2 ? 


the corresponding p.m.f for W, is given by 


i; if 7~=1,2 
f@)= 4 =; if 1=3,4, 
0; elsewhere. 


Hence, the corresponding y*-chromatic mean is 


n—2 n-2 T-—8 


1 1 
my el —-4+2-5 2 ee 
H+ (Wa) n * n +3 2n - 2n 2n 


and the y+-chromatic variance is 
xX 


1 n—-2 
2 2 2 2 2 2 
Ox+(Wn) = (1° +2°)-— + (38° +4") 3 — (ox (Wn) 
1 n—2 7n—8\" n24+32n—64 
= .2 495. = eos See 
(s mg 2 —) ( Qn ) 4n2 


2.3. Some Interpretations 


A block graph or clique tree G is an undirected graph in which every biconnected component 
(block) is a clique. By Theorem 2.7, minimum proper coloring of every component of G follows 
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uniform distribution. Hence, we have 


Theorem 2.16 The probability distribution of a block graph G is mixture of discrete uniform 


distributions. 


An n-partite graph is a graph whose set of vertices can be partitioned in to n subsets such 


that no two vertices in the same partitions are adjacent. Then, we have the following result. 


Theorem 2.17 Let G be a regular k-partite graph on vertices. Then, any minimal proper 


coloring of G follows uniform distribution (in each partition). 


proof Any minimal proper coloring of a k-partite graph contains k-colors. Let G be an 
r-regular k-partite graph. Then, rk = n. Then, the p.m.f of G is 


i=1,2,3,...,k, 


elsewhere. 


Ay, 
fsx © 
0; 


which is that of the DU(k) distribution. 


Corollary 2.18 Let G be a k-partite graph. Then, the x-chromatic mean (and yt -chromatic 


. . + . . . 2 
mean) of G is se and the x-chromatic variance (and x*-chromatic variance) of G is a 


Proof The proof follows immediately from the fact that the minimal proper coloring of a 


k-partite graph follows uniform distribution. 


Certain areas where these notions can be made use of are: nodes in communication and 


traffic networks. 


§3. Scope for Further Studies 


In this paper, we have extended the notions of mean and variance to the theory of graph coloring 
and determined their values for certain graphs and graph classes. More problems in this area 
are still open. 

The y-chromatic mean and variance of many other graph classes are yet to be studied. 
Determining the sum, mean and variance corresponding to the coloring of certain generalized 
graphs like generalized Petersen graphs, fullerence graphs etc. are some of the promising open 
problems. Studies on the sum, mean and variance corresponding to different types of edge 
colorings, map colorings, total colorings etc. of graphs also offer much for future studies. 

We can associate many other parameters to graph coloring and other notions like covering, 


matching etc. All these facts highlight a wide scope for future studies in this area. 


Acknowledgement 


The first author of this article dedicates this paper to the memory Prof. (Dr.) D. Balakrishnan, 


98 


N.K. Sudev, K.P. Chithra, S. Satheesh and Johan Kok 


Founder Academic Director, Vidya Academy of Science and Technology, Thrissur, India., who 


had been his mentor, the philosopher and the role model in teaching and research. 


References 


[1] M. Batsyn and V. Kalyagin, An analytical expression for the distribution of the sum of 


Oo OND Oe 


10 


11 


12 


13 


14 


random variables with a mixed uniform density and mass function, In Models, Algorithms, 
and Technologies for Network Analysis (Editors: B. I. Goldengorin, P. M. Prdalos and V. 
Kalyagin), 51-63, Springer, 2012. 

J. A. Bondy and U.S. R. Murty, Graph theory with application, North-Holland, New York, 
1982. 

A. Brandstaédt, V. B. Le and J. P. Spinrad, Graph Classes: A survey, SIAM, Philadelphia, 
1999. 

G. Chartrand and P. Zhang, Introduction to Graph Theory, McGraw-Hill Inc., 2005. 

G. Chartrand and P. Zhang, Chromatic graph theory, CRC Press, 2009. 

F. Harary, Graph theory, Addison-Wesley Pub. Co. Inc., Philippines, 1969. 

T. R. Jensen and B. Toft, Graph Coloring Problems, John Wiley & Sons, 1995. 

M. Kubale, Graph Colorings, American Mathematical Society, 2004. 

J. Kok, N. K. Sudev and K. P. Chithra, General coloring sums of graphs, Cogent Math., 
3(1)(2016), 1-11, DOT: 10.1080/23311835.2016.1140002. 

E. Kubicka and A. J. Schwenk, An introduction to chromatic sums, Proc. ACM Computer 
Science Conference, Louisville (Kentucky), 3945(1989). 

E. Kubicka, The chromatic sum of a graph: History and recent developments, Int. J. Math. 
Math. Sci., 30,(2004), 1563-1573. 

V. K. Rohatgi, A. K. Md. E. Saleh, An Introduction to Probability and Statistics, Wiley, 
New York, 2001. 

S. M. Ross, Introduction to Probability and Statistics for Engineers and Scientists, Aca- 
demic Press, 2004. 

N. K. Sudev, K. P. Chithra and J. Kok, Certain chromatic sums of some cycle related graph 
classes, Discrete Math. Algorithms Appl., 8(3)(2016), 1-24, DOI: 10.1142/S1793830916500- 
506. 


[15] E. W. Weisstein, CRC Concise Encyclopedia of Mathematics, CRC Press, 2011. 
[16] D. B. West, Introduction to Graph Theory, Pearson Education Inc., 2001. 


International J.Math. Combin. Vol.3(2018), 99-107 


On Status Indices of Some Graphs 


Sudhir R.Jog 


(Department of Mathematics, Gogte Institute of Technology, Udyambag Belagavi, Karnataka, 590008, India) 


Shrinath L. Patil 


(Department of Mathematics, Hirasugar Institute of Technology, Nidasoshi, Karnataka, 591236, India) 


Email: sudhir@git.edu, shrinathlpatil@gmail.com 
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Key Words: Status of vertex, first status connectivity index, second status connectivity 


index, harmonic status index. 


AMS(2010): 05C12. 


§1. Introduction 


There are several molecular structural graph descriptors such as Weiner Index, Zagreb Index, 
Hosoya Index etc which strongly correlate studies in graph theory with chemistry. Most of these 
indices are based on the distance between vertices in a graph.Motivated by harmonic mean we 
have harmonic index of a graph defined by Fajtlowicz [5]. For more work one can refer [6]. 


Further motivated by the same, Ramane and Yalnaik introduced the harmonic status index of 


graphs [4]. 


Definition 1.1((1]) The status of a verter u € V(G) is defined as the sum of its distance from 
every other vertex in V(G) and is denoted by o(u). That is 


a(u) = Ss d(u, v). 


ueV(G) 


Definition 1.2 The first status connectivity index S\(G)and second status connectivity index 


S(G) of a connected graph G are defined respectively as 


S(G)= S© [o(u) +o(v)] and $2(G)= YS > [o(u)o(v)]. 


uv€ E(G) uv€ E(G) 
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Similarly the first and second status connectivity coindices of a connected graph G are 
defined as 
Si(@)= So [o(u)+o(v)| and 82(@)= S° [o(ujo(r)].- 


uv¢ E(G) uv¢ E(G) 


Definition 1.3([5]) The Harmonic index of a graph G is defined as 
H(G) = : — 
ee d(u) + d(v) 


The harmonic status index of a connected graph G as ([4]) 


2 


a ae) ew). 


uve E(G) 


Similarly the harmonic status coindex of a connected graph G is defined as 


§2. Status Connectivity Indices and Coindices of Some Graphs 


In what follows, we consider a class of graphs constructed by first joining a path of length 
I(> 1) to each vertex of G and then attaching k pendent vertices to each end vertex of the 
path attached. Such a graph can be called / level thorn graph denoted by G'(+*), The usual 


thorny graph Gt* can be regarded as 0 level thorn graph. If / = 1 we get first level thorn graph 
GUC+F) | 


Example 2.1 A graph G and it’s first level thorn graph G*!(+®) are as shown below. 


Figure 1 


First we evaluate the status connectivity index and coindex of 0 level thorn graphs denoted 
by Gt*.To obtain the harmonic status connectivity index and coindex of this graph we need 
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to calculate status of each vertex and number of pairs of adjacent vertices and pairs of non 
adjacent vertices in Gt*.If G is a r regular graph then,with respect to degree there are two 
types of vertices in G+" nk pendent vertices (external),n vertices of degree 'r+nk’ we call them 
as internal. 


Theorem 2.1 The first and second status connectivity index of thorn graph K;*" are given by 


cg) 

F 

5 
I 


n(n — 1)(Q2nk +n—-—k-1)+nk(5nk + 3n — 2k — 4) 
(nk)C2 x (8nk + 2n — k — 3)? + nC x (8nk + 2n—k—3)(2nk+n—k-1) 


we 
e 

ee 
| 


Proof The graph K;* is of diameter3and there are two types of vertices in it. A set of ‘nk’ 


pendent vertices and n vertices of degree /n + 1’. Let uj, 7 = 1,2,---,nk denote the pendent 
vertices and v;,i = 1,2,--- ,n denote the vertices of degree ’n + 1’. Then the status of pendent 
vertex is 


o(ui) = 1+ 2(k — 1) + 2(n— 1) + 3k(n — 1) = 8nk + 2nk —3 
and the status of the internal vertex vu; is 
o(u;) =1(n—1) +k + 2k(n—1) = (Qnk+n—k-—1). 


Now in K;* there are nina) adjacent pairs internal vertices and nk pairs of vertices 


forming edges formed by one internal and one external vertex. Hence by definition the status 


connectivity index of K;t" is 


~1)2(2nk +n—k-1 
Sy(Kt*) = MPO RD nk(3nk + 2n— k- 3+ 2nk+m—k—I) 


= n(n—1)(Qnk+n—k—-1)+nk(5nk + 3n — 2k — 4). 


Also in K;** there are (nk)C2 pairs of nonadjacent pendent vertices and nk(n — 1) pairs of 
nonadjacent pairs of vertices formed by one pendant and one internal vertex. So that status 
connectivity coindex of K;*" is 


So(K,t*) = (nk)Co x (8nk + 2n — k — 3)? + nCg x (8nk + 2n —k —3)(Qnk+n—k-1). 


Theorem 2.2. The harmonic status index and coindex of thorn graph K** are given by 


—1) 1 2 
H Ktk = n(n _—————— _——— 
S(Kit*) >. (pena eo Garber oe a), 
— 1 2 
a en a eo | nn ry 
She) oe OR to, oR) aR on SE) 


Proof The graph K;** is of diameter 4 and there are two types of vertices in it. A set 
of ‘nk’ pendent vertices and n vertices of degree 'n + 1’. Let uj,7 = 1,2,---,nk denote the 
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pendent vertices and v;,i = 1,2,--- ,n denote the vertices of degree n+ 1’. Then the status of 
pendent vertex is 
o(uj) =14+ 2(k — 1) + 2(n — 1) + 8k(n — 1) = 8nk 4+ 2n—-—k—-3 
and the status of the internal vertex v; is 


o(u;) = 1(n—-1) +k + 2k(n—- 1) = 2nk+n—-—k-1. 


Now in K,** there are mind) adjacent pairs internal vertices and nk pairs of vertices forming 


edges formed by one internal and one external vertex. Hence by definition the harmonic status 
index of K;** is 


=i) 2 2 

HS(kt*) = mm-1 eee’! 

ous) 2 2Qnkan—k=1) | Gnktn—k—S342nktn—k=D 
a n(n — 1) 1 Li 2 


2 (Qnk+n—k-—1) (5nk + 3n — 2k — 4) 


Also in K;** there are (nk)C2 pairs of nonadjacent pendent vertices and nk(n— 1) pairs of 
nonadjacent pairs of vertices formed by one pendant and one internal vertex. So that harmonic 
status coindex of K;t* is 

eC) ——— —— a ————— 
™ 2(38nk + 2n — 2k — 4) (8nk + 2n —k -—3+4+2nk+n—k-—1) 
1 2 


=: tks ——— 4. nk — 1) ___*~* 
Ca Ron Deaw. Vigne en k=) 


Now, we discuss the status connectivity indices and the coindices of regular graphs with 
diameter 2. 


Theorem 2.3 If G is ‘r’ regular graph of diameter 2 then the first and second status connectivity 


index of Gt* are given by 


Si(Gt*) = nr(Qn4+ 2kr+k—r—2)+nk(5n + 5kr + 3k — 3r — 6), 
S(Gt*) = (an t Qkr + k—r— 2)? + nk(3n+3kr + 2k — 47 — 2). 


Proof The proof follows by direct counting. 


Theorem 2.4 If G is ‘r’ regular graph of diameter 2 then the first and second status connectivity 


co index of G*® are given by 


Si(@t*) = 2(3n + 3kr + 2k —r — 4) + nk(n — 1)(5n + 5kr + 3k — 3r — 6) 
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= nk(nk — 1)(8n+ 3kr + 2k — r — 4) + nk(n — 1)(5n + 5kr + 3k — 3r — 6) 


nmr 


+(nC2 — =)(4n + 4kr + 2k — 2r — 4), 
_— k(nk — 1 
Sx(Gtk) = ONO) an + 3k + 2k — r — 4)? 
tnk(n — 1)(8n + 3kr + 2k — 2r — 4)(2n + 2kr +k — 1 — 2) 


+(nC2 — >)(2n +2kr+k—r—2)?. 


Proof The proof follows by direct counting. 


Theorem 2.5 If G is ‘r’ regular graph of diameter 2 then the harmonic status index of Gt* is 


1 2 
H SS rs k———____ > 
Me) = S. Onb ahr k=r ao) ' Gat bere Gh oP =O) 


Proof First, we observe that if G has diameter 2 then Gt* has diameter 4. Hence from 


the structure we have the status of each internal vertex v; as 
o(u;) = 1(k +r) + 2kr+2.(n-1—r) =2n+4+ 2kr+k—r—2. 
Also the status of each pendant vertex u; as 


o(uj) =14+2.r4+ 2(k-1)4+3(n-—1—-r) =3n+4+ 3rk4 2k-r—-4. 


nr 


There are + internal edges giving harmonic status contribution 


nr 2 nr 1 


“2 20Qn+2%kr+k—r—2) 2 Qn+22kr+k—r—2) 
Similarly the pendent ‘nk’ vertices adjacent to ‘n’ internal vertices contribute, 


k 2 
rk—_—_——_—__——_.. 
(5n + 5kr + 3k — 2r — 6) 


Hence the harmonic status index of G** is 


1 2 
ti. ————— ee 
Oo One Dare), Greene bak Sor a0) 


Theorem 2.6 If G is ‘r’ regular graph of diameter 2 then the harmonic status coindex of Gt* 
18 
1 2 nr 


ee oy, 
(Qn +3kr+2k—-7r—-4) 1” Ces) ngs ) 


FB(G**) = (nkC2) GhaOhs howe) OOF 


Proof We note that there are n(k + 1)C2 — (4 + nk) non adjacent pairs of vertices in 
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Gt, There are (nk)C2 pendent nonadjacent pendent vertices, nk(n — 1) pairs of nonadjacent 


nr 


vertices combining one pendant and one internal vertex and finally (nC2 — 3%) nonadjacent 


2 


internal vertices. Taking contribution from each of them we have status connectivity coindex 


of Gt* as 


HS(G**) 


2 2 


| 
WW) a ee ap ape eee gee 


nr 2 
Ces) somos ESD 
1 2 
eS Ye, | ae 
Grae onaay TY eae aka) 
nr 1 


TMCS) ap Oke: RP Oy 


§3. Status Connectivity Indices and Coindices of First Level Thorn Graphs 


Now we discuss the harmonic status index and coindex of first level thorn graphs.We need to 


calculate status of each vertex and number of pairs of adjacent vertices and pairs of non adjacent 


vertices in G!(+*) With respect to degree there are three types of vertices in G*!(+") nk 


pendent vertices,n vertices of degree ’k + 1’ we call them as internal and lastly ‘n’ vertices 


having degree sequence added by 1.We call them external, in particular if G is ‘r’ regular their 


degrees will become ‘r + 1’. 


Theorem 3.1 The first and second status connectivity index and coindex of first level thorn 


graph of a‘r’ regular graph of order ‘n' and diameter 2 are given by 


Si(G C+) y = 


Sx(GIGH) = 


Proof The proof follows by direct counting. 


ZX (on A4nk — rk — 2r — 2k — 4) + n(12n + 9nk — rk — 4r — 4k — 10) 


4+2n? x k(7n + 5nk — rk —2r — 3k — 6), 

— x (5n+4nk — rk — 2r — 2k — 4)? 

n(5n + 4nk — rk — 2r — 2k — 4)(7n + 5nk — 2r — 4k — rk — 6), 
(nk)C2 x 2(9n + 6nk — rk — 4k — 2r — 8) 

+(nCz — =) x 2(5n + 4nk — rk — 2r — 2) 


+2(nC2)(7n + 5nk — rk — 2r — 3k — 6) + n?k(7n + 5nk — rk — 2r — 3k — 6) 
+2n(n — 1)(12n + 9nk — 2rk — 4r — 6k — 10) 
(nk)Cz x (9n + 6nk — rk — 4k — 2r — 8)? 


+(nC2 — S)(6n + Ank — rk — 2r — 2)? 


4+nCo(7n + 5nk — rk — 2r — 3k — 6)? + n?k(7n + 5nk — rk — 2r — 3k — 6)? 
+n(n— 1)(12n + 9nk — 2rk — 4r — 6k — 10). 
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Theorem 3.2 If G is ‘r’ regular graph of diameter 2 then the harmonic status connectivity 
index of G+) is 


1 
H Al(tk)y 
a ) 2 (5n + 4nk — rk — 2r — 2k — 4) 
2 
"Gn + 4nk — rk — 2r — 2k —4+4 7n + 5nk — 2r —4k —rk —6) 
2 
+n7k 


(5n + 4nk — rk — 2r — 2k —44 9n+4 6nk — rk — 4k — 2r — 8) 


Proof First, we observe that if G has diameter2 then G!(+") has diameter 6. Hence from 


the structure we have the status of each internal vertex v; as 


o(u;) = lrt+1)4+2(r+k)4+2.n-—1l—-1r)+38rk+3(n—1—r)+4k(n-1-71) 
= 5n+4nk—rk—2r—2k—4. 


Also the status of each external vertex u; as 


l 


1(k+1)+ 2r+ 3r+3(n—1—r)+4(n-—1-—r)+4rk + 5k(n-—1-1r) 
= n+ 5nk— 2r—4k—rk—6. 


o(ui) 


Finally the pendent vertices being the only vertices on the diametrical path have the status 


14+2x14+2(k—-1)+3r4+4(n—-—1—1r)+4r4+5kr+5(n—-1-r)4+6k(n—-1—-r) 
= 9n+6nk —rk — 4k — 2r—8. 


o(wi) 


In G*'(+*) there are “ pairs of internal adjacent vertices, n pair of adjacent vertices 
7 P J J 


formed of one internal and one external vertex and finally n?k pairs of adjacent vertices formed 
of one internal and one pendant vertex.Hence the harmonic status index of G!(+") is given by 


Al(+k nr 1 
HS(GN) "2 (5n + 4nk — rk — 2r — 2k — 4) 
2 
"Gn + 4nk — rk — 2r — 2k —4+4 Tn + 5nk — 2r — 4k — rk — 6) 
2p, 2 
7” “Gn + 4nk — rk — 2r — 2k —4 + 9n + Onk — rk — 4k — or — 8) 
_ or 1 
~ 2 (5n + 4nk — rk — 2r — 2k — 4) 
+n g +n7k 2 


(12n + 9nk — 2rk — 4r — 4k — 10) (14n + 10nk — 2rk — 4r — 6k — 12) 

_onr 1 dg 2 

~ 2 (5n + 4nk — rk — 2r — 2k — 4) (12n + 9nk — 2rk — 4r — 4k — 10) 
1 

(7n + 5nk — rk — 2r — 3k — 6) 


+n7k 
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Theorem 3.3. The harmonic status coindex of G\(+*) is given by 


SiGe) 


i 
BY Cpa ee et 
WON)C3 one Ga Sek Aa Ope 
nr 1 
Sc a (rman Po RG | 
1 . 1 
a a a ae 
OC nk ee ae =O) | in Pane k= = R=) 
+n(n —1) = 


(12n + 9nk — 2rk — 4r — 6k — 10) 


Proof In G+") there are (nk)C2 pairs of nonadjacent pendent vertices, 


(nCz — 4) pairs of nonadjacent vertices formed by internal vertices, nC2 pairs of nonadjacent 


vertices formed by external vertices, n?k nonadjacent pair of vertices formed by one pendant 


and one internal vertex and finally n(n—1) pairs of nonadjacent vertices formed by one internal 


and one external vertex.Hence the harmonic status connectivity coindex is given by 


HS(G*(+4)) 


§4. Conclusion 


2 
k ee ea 
(nk) C25 + Gnk rk 4h — 2 8) 
nr Ds 
tine <*>) 5G hank Sho aD 
2 
STC 2507 Sak oh =a = 6) 
2 
2 fp 
+?" Gn + 4nk— rk — Or —2k—449n + Onk— rk —4k— Or —8) 
+n(n—1) : 
mn — a0 Wo 5500 
(5n + 4nk — rk — 2r — 2k —4+4 7n+ 5nk — 2r — 4k — rk — 6) 
1 
k a eae 
nl) CP oa Gak Tk aR oF 8) 
nr 1 
TNE =o ea ek ODED) 
1 1 
_— — App 
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2 
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(12n + 9nk — 2rk — 4r — 6k — 10)’ 


We considered general | level thorn graphs and obtained in particular, status connectivity 


indices and coindices as well as Harmonic status indices and coindices of 0 level and first level 


thorn graphs for some class of graphs. 
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Abstract: Representing a subset of vertices in a graph by means of a matrix was introduced 
by £. Sampath Kumar. Let G(V,E) be a graph and S C V be a set of vertices. We can 
represent the set S by means of a matrix as follows, in the adjacency matrix A(G) of G 
replace the a;; element by 1 if and only if, v; € S. In this paper we study the set S being 


dominating set and corresponding domination energy of some class of graphs. 


Key Words: Adjacency matrix, Smarandachely k-dominating set, eigenvalues, energy of 


graph, distance energy, Laplacian energy. 
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§1. Introduction 


A set D C V of G is said to be a Smarandachely k-dominating set if each vertex of G is 
dominated by at least k vertices of S and the Smarandachely k-domination number y;,(G) of G 
is the minimum cardinality of a Smarandachely k-dominating set of G. Particularly, if k = 1, 
such a set is called a dominating set of G and the Smarandachely 1-domination number of G is 
called the domination number of G and denoted by y(G) in general. 

The concept of graph energy arose in theoretical chemistry where certain numerical quan- 
tities like the heat of formation of a hydrocarbon are related to total 7 electron energy that 
can be calculated as the energy of corresponding molecular graph. The molecular graph is a 
representation of the molecular structure of a hydrocarbon whose vertices are the position of 
carbon atoms and two vertices are adjacent if there is a bond connecting them. 

Eigen values and eigenvectors provide insight into the geometry of the associated linear 
transformation. The energy of a graph is the sum of the absolute values of the Eigen values of 
its adjacency matrix. From the pioneering work of Coulson []1 there exists a continuous interest 
towards the general mathematical properties of the total 7 electron energy ¢€ as calculated within 
the framework of the Huckel Molecular Orbital (HMO) model. These efforts enabled one to 
get an insight into the dependence of ¢ on molecular structure. The properties of ¢(G) are 
discussed in detail in [2, 3, 4]. 

The importance of Eigen values is not only used in theoretical chemistry but also in ana- 


lyzing structures. Car designers analyze Eigen values in order to damp out the noise to reduce 
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the vibration of the car due to music. Eigen values can be used to test for cracks or deformities 
in a solid. Oil companies frequently use Eigen value analysis to explore land for oil. Eigen 


values are also used to discover new and better designs for the future. 


§2. Definitions and Notations 


Representation of a subset of vertices of a graph by means of a matrix was first introduced 
by E.Sampath Kumar [5]. Let G(V, £) be a graph and S C V be a set of vertices. We can 
represent the set S by means of a matrix as follows: 


In the adjacency matrix A(G) of G replace the a;; element by 1 if and only if vu; € S. 
The matrix thus obtained from the adjacency matrix can be taken as the matrix of the set S 
denoted by Ag(G). The energy E(G) obtained from the matrix Ag(G) is called the set energy 
denoted by E's(G). In this paper we consider the set S as dominating set and the corresponding 
matrix as domination matrix denoted by A,(G) of G. Thus the energy E(G) obtained from 
the domination matrix A,(G) is defined as domination energy denoted by E,(G). 


Let the vertices of G be labeled as v1, v2,v3,-++,Un. The domination matrix of G is 
defined to be the square matrix A,(G) corresponding to the dominating set of G. The Eigen 
values of the domination matrix denoted by 1, K2,K3,°-*,Kn are said to be the A, Eigen 
values of G. Since the A, matrix is symmetric, its Eigen values are real and can be ordered 


Ky > kg > Kk3 S++: > Ky. Therefore, the domination energy 
EE AG)= >. mil (1) 


This equation has been chosen so as to be fully analogous to the definition of graph energy 


([2]). : 
E = E(G) = dla, (2) 


where 41 > Ap > Az S--: S An are the Eigen values of the adjacency matrix A(G). Recall 
that in the last few years, the graph energy E(G) and domination energy [9,10] or covering 
energy ({6]) has been extensively studied in the mathematics ([6,7]) and mathematic-chemical 
literature ([{8,12]). 


Definition 2.1(Minimal domination energy) A dominating set D in G is a minimal dominating 
set if no proper subset of D is a dominating set. The domination energy E(G) obtained for a 


minimal dominating set is called the minimal domination energy denoted by Ey—min(G). 


Definition 2.2(Maximal domination energy) A dominating set D inG is a maximal dominating 
set if D contains all the vertices of G. The domination energy E,(G) obtained for a maximal 


dominating set is called the maximal domination energy denoted by Ey—max(G). 


Similarly to domination energy of graph G, distance domination energy, Laplacian domi- 
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nation energy and Laplacian distance domination energy can also be defined as follows. 


Let the vertices of G be labeled as v1, v2, 3,-++ ,Un. The distance matrix of G, denoted by 
D(G) is defined to be the square matrix D(G) = [d,;], where d;; is the shortest distance between 
the vertex vu; and v; in G. The Eigen values of the distance matrix denoted by p11, f2, 13,°°* 5 bn 
are said to be the D Eigen values of G. Since the D(G) matrix is symmetric, its Eigen values 
are real and can be ordered fu, > pg > pg > +--+ S fn. Therefore, the distance energy 


Ep = Ep(G) = D [pil (3) 


In the distance matrix D(G) of G replace the a;; element by 1 if and only if v; € S. The 
matrix thus obtained from the distance matrix can be considered as the distance matrix of the 
set S denoted by Ds(G). The energy E(G) obtained from the matrix Ds(G) is called the 
distance set energy denoted by Dg(G). In this paper we consider the set S as dominating set 
and the corresponding matrix is distance domination matrix denoted by D,(G) of G. Thus 
the energy E(G) obtained from the distance domination matrix D,(G) is defined as distance 
domination energy denoted by Ep,(G). 


The distance domination matrix of G is defined to be the square matrix D.,,(G) correspond- 
ing to the dominating set of G. The Eigen values of the distance domination matrix denoted 
by 01,02, 03,°*+ ,On are said to be the D, Eigen values of G. Since the D,(G) matrix is sym- 
metric, its D-Eigen values are real and can be ordered as 01 > 02 > 03 >-+:: > On. Therefore, 


the distance domination energy 


Epy = Ep,(G) = De lil. (4) 


Definition 2.3(Minimal distance domination energy) A dominating set D in G is a minimal 
dominating set if no proper subset of D is a dominating set. The distance domination energy 
Ep,(G) obtained for a minimal dominating set is called the minimal domination energy denoted 
by Epa min (@). 


Definition 2.4(Maximal distance domination energy) A dominating set D in G is a maximal 
dominating set if D contains all the vertices of G. The distance domination energy Ep+(G) 
obtained for a maximal dominating set 1s called the maximal domination energy denoted by 


Epy—max(G) : 


Let the vertices of G be labeled as v1, v2, U3,--+ , Un. The Laplacian matrix of G is denoted 
by L(G) is defined to be the square matrix L(G) = d(G) — A(G), where A(G) and d(G) are 
the adjacency matrix and diagonal matrix with vertex degree of G on the principal diagonal 
element respectively. The Eigen values of the Laplacian matrix denoted by W, W2, w3,°-: Wn 
are said to be the L Eigen values of G. Since the L(G) matrix is symmetric, its Eigen values 
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are real and can be ordered yy, > w2 > v3 S++: S Wn. Therefore, the Laplacian energy 
Ex = Ex(G) => |vil. (5) 
i=1 


The energy E,,(G) obtained from the matrix L5(G) = d(G)—Ags(G) is called the Laplacian 
set energy denoted by Ls(G). In this paper we consider the set S as dominating set and the 
corresponding matrix is Laplacian domination matriz denoted by L(G) of G. Thus the energy 
E(G) obtained from the Laplacian domination matrix L(G) is defined as Laplacian domination 
energy denoted by E,,(G). 

The Laplacian domination matrix of G is defined to be the square matrix L(G) corre- 
sponding to the dominating set of G. The Eigen values of the Laplacian domination matrix 
denoted by a1, a2,03,:-+ ,Qm are said to be the L., Eigen values of G. Since the L(G) matrix 
is symmetric, its L-Eigen values are real and can be ordered as ay > a2 > a3 > °-: S Qn. 
Therefore, the Laplacian domination energy 


EL, = Ey(G) = Dy la]. (6) 


Definition 2.5(Minimal lapalcian domination energy) A dominating set D in G is a minimal 
dominating set if no proper subset of D is a dominating set. The Laplacian domination energy 
Ez,+(G) obtained for a minimal dominating set is called the minimal domination energy denoted 
by Fis min(G): 


Definition 2.6(Maximal lapalcian domination energy) A dominating set D in G is a maximal 
dominating set if D contains all the vertices of G. The Laplacian domination energy E,,(G) 
obtained for a maximal dominating set is called the maximal domination energy denoted by 


Ea tax (G) : 


The energy Eypy(G) obtained from the matrix LDs(G) = d(G) — Ds(G) is called the 
Laplacian distance set energy denoted by LDg(G). In this paper we consider the set S as 
dominating set and the corresponding matrix is Laplacian distance domination matrix denoted 
by LD,(G) of G. Thus the energy E(G) obtained from the Laplacian distance domination 
matrix L.D,(G) is defined as Laplacian distance domination energy denoted by Ey p+(G). 

The Laplacian distance domination matrix of G is defined to be the square matrix LD,(G) 
corresponding to the dominating set of G. The Eigen values of the Laplacian distance domi- 
nation matrix denoted by (1, G2, 83,--- , Gn are said to be the LD, Eigen values of G. Since 
the L.D,(G) matrix is symmetric, its L-Eigen values are real and can be ordered as 6; > G2 > 
B3 > +++ > By. Therefore, the Laplacian distance domination energy 


Expy = Expy(G) = Sa |i]. (7) 
i=l 


Definition 2.7(Minimal Lapalcian distance domination energy) A dominating set D in G is 


a minimal dominating set if no proper subset of D is a dominating set. The Laplacian dis- 
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tance domination energy E,p (G) obtained for a minimal dominating set is called the minimal 


domination energy denoted by E_py—min(G). 
Definition 2.8(Maximal Lapalcian distance domination energy) A dominating set D in G is a 
maximal dominating set if D contains all the vertices of G. The Laplacian distance domination 


energy E,py(G) obtained for a maximal dominating set is called the maximal domination energy 
denoted by Erpy—max(G). 


§3. Various Domination Energies 


Definition 3.1 A book graph (B,,) consists of m quadrilaterals sharing a common edge. That 
is, it is a Cartesian product Siyj41 and Pz, where Sm is a star graph and P, is the path graph 


on two nodes. Some book graphs are shown in Figure 1. 


B3 B4 Bs Be 


PRK IK 


Figure 1 Book graph By, 3<m<6 


Theorem 3.1 For m > 3, the minimum dominating energy of a book graph (Bm) is 


2(Vv4m+1+m-—1). 


Proof Calculation enables one to find the characteristic polynomial of B,, for m > 3 


directly. 
For m = 3, Bs is a book graph with 8 vertices. The minimum dominating set is S = 
{U1, V2}- V3 U4 
U1 v2 
U5 U6 U7 U8 


Figure 2 Book graph Bs 


Calculation shows that the domination matrix and the characteristic polynomial of B3 are 
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respectively given by 


1 11031 041 +0 
1 103101 041 
100 1 0 0 0 0 
0 1 1 00 0 0 0 
Ay (G) = 
100 0 0 1 0 0 
0 1 00 1 0 0 0 
100 0 0 0 0 1 
0 1 00 0 0 1 0 


and K8—2K7—96°412K°+ 1844-1843 13K2+86+3 = (& — 1)? (+1)? (62—-3«—1)(K2+K—3). 


And calculation shows that the domination matrix and the characteristic polynomial of By 
are respectively given by 


oOo KF CO fF CO FF OF KF KF F&F 
- CO fF CO FF OO FF CO FF & 
oo OoOcUCOOUUCUCOUUCOOUCUCrRRMDLUCOUCOCOCUCUr 
(oe > > 
oO oO.cUCDUCUOUOUC rh rMrlLUC COUCOCOCCOUO 
Co CocmUCcrUmUCCcCTCOlUCUCUlrMDm Cc OUCOUrULUCOD 
[co I SI aR > SE > ORY E> A CD J CRED TR 
oOo oO 0 RPFrlLUCOUCOCUCUCOUOUCUCOUNUC RLU 
Peoooooeoeoo oO KF 
oOoroeoooeoorF ©& 


and «9 —269— 1248+ 16K? +38%°— 364° — 526443263 +33K2—106—8 = (« — 1)? (« +.1)° (x? 
3k — 2)(K? + — 4). 


Similarly, the domination matrix and the characteristic polynomial of Bs are respectively 


given by 
1 1103101 0 1 0 1 ~0 
1 10310 1 03101 0 21 
100 10 00 0 0 0 0 0 
0 11 000 00 0 0 0 0 
100 00 1 0 0 0 0 0 0 
cers 0 10 01 0 00 0 0 0 0 
1 00 00 00 1 0 0 0 0 
0 10 00 0 1 0 0 0 0 0 
100 00 00 0 0 1 0 0 
0 10 00 0 00 1 0 0 0 
100 0 0 00 0 00 0 1 
0 10 00 0 00 0 0 1 +0 
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and (« — 1)* («+ 1)4 (Kk? — 3« — 3)(K2 + « — 5), respectively. 


And the characteristic polynomial of Bg is given by 


(«% — 1)? («& +1)? (a? — 8% — 4)(n? + & — 6) 
Generally, the characteristic polynomial of B,, using domination adjacency matrix is 


(B= 4)" ely = SS (MR 2)? +e Om): 


Solving the equation we get 


(« —1)""* =0, or (h +1)” * =0, or (Kk? — 36 — (m— 2)) = 0 or (x2 +4 —m) =0. So 
& =1,1,1,---,1 ((m— 1)times), or « = —1,—-1,—-1,--- ,—1((m — 1)times). 


By (kK? — 3k — (m — 2)) = 0, we get 


K1 = 5 (3- Vim F) and 
ko = 5 (34+ Vim F1) here m > 3. 


By (Kk? + —m) =0 we know that 


K3 = (-1- 4m +1) and 


(-1+ V4m + 1) 


K4 = 


Dl eR ple 


Hence, 


Ey—min = ymin (G) = 3 || 
t=1 


‘ 


Therefore, 
Ey—min = y—min (Bm) = 2(V4m+1+m-—1). 


This completes the proof. 


Theorem 3.2 For m > 3, the minimum distance domination energy of a book graph (By) is 


A(m — 1) + V25m? — 24m + 36 + \/mV/m + 4. 


Proof Calculation enables one to find the characteristic polynomial of B,, for m > 3 
directly. 
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For m = 3, B3 is a book graph with 8 vertices. The minimum dominating set is S = 
{v1, v2}. Calculation shows that the distance domination matrix and the characteristic poly- 
nomial of Bs are respectively given by 


ihe de, ay, 22a De uly 22 
1121212 21 
12012 3 2 3 
21103 2 3 2 
D,(G) = 
122 3 01 2 3 
213 2 1 0 3 2 
122 3 2 3 0 1 
2 AL. 33) 28 3. 2 TO 


and 08 —207—1110%—5120°—54504+50403+2400? = 0? (a + 4)* (o? — 130 — 5) (0? + 30 — 3). 
Similarly, calculation shows that the distance domination matrix and the characteristic 


polynomial of By, are respectively given by 


111212 12 «1 2 
Te Wide Dee We De NY SPT A 
12012 3 2 3 2 3 
2 de tO? Bh De 3) De Be °2 
TP! 22 2? 28e CO Dh B28 2e 3 
D,(G) as 
213 2 10 3 2 3 2 
122 3 23 01 2 3 
213 2 3 2 10 3 2 
122 3 23 23 0 1 
213 2 3 2 3 2 «21 «0 


and o!° — 20° — 20008 — 151207 — 40480 — 224005 + 435204 + 102403 = 03 (0 +4)? (0? 
180 — 4) (o? + 40 — 4). 


And the characteristic polynomial of Bs and Bg are respectively given by 


Generally, the characteristic polynomial of B,, using the distance domination matrix is 
eo oaaye [o? — (5m — 2)o + (m— 8)| (0? + mo —m) =0. 


Solving the equation we get 
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o”—-! = 0, or (o +4)” | =0, or (0? — (5m—2)o+(m—8)) =0, or (0? +ma—m) =0. So 
¢ =.0;0,.0;- we — 1)times), or o = —4,—4, —4,--- , -4 ((m — 1)times), and (0? — (5m — 


o1 = 5 (5m— Dam 25m? — 24m + 36) and 
on = 5 (5m—24 Vm? — 2m +36) here m > 3, 
Bee 
93 = 5 (-m—Vmvm-+4) and 
oa = 5 (-m+ vinvin Fa) 


E'p+—min a Ep+—min (G) 


=a(m—1) +|5 (2 Bim — Bam +30) | + E (2Vmv/m + 4)}. 
Therefore, 


Epy—min = Epy—min(G) = 4(m — 1) + 25m? — 24m + 36 + VmvVm + 4. 


This completes the proof. 


Theorem 3.3 Form > 3, the minimum Laplacian domination energy of a book graph (Bm) is 


5m + Vm? +4. 


Proof Calculation enables one to find the characteristic polynomial of B,, for m > 3 
directly. 


For m = 3, Bs is a book graph with 8 vertices. The minimum dominating set is S = 
{vi, v2}. The Laplacian domination matrix and the characteristic polynomial of B3 are respec- 


tively calculated by 


-1 3 0 -1 0 -1 0 -1 
-1 0 2 -1 0 0 0 0 
L(G) = Oo -1l -l 2 0 0 0 O 
-1 0 0 0 2 -1 0. 0 
0 -l 0 0 -1 2 0 O 
-1 0 0 0 0 0 2 -1 
0 -l 0 0 0 0 -1 2 
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and a® — 18a7 + 131a® — 496a° + 103804 — 11540* + 5430? + 36a —81 = (a — 1)? (a — 3)? (a? — 
7a + 9)(a? — 3a —1). 


Similarly, the Laplacian domination matrix and the characteristic polynomial of By, are 
respectively given by 


-1 4 0 -1 0 -1 0 -1 0 -!I 
-1 0 2 -1 0 0 0 0 0 0 
O -1 -1 2 0 0 0 0 0 0 
ie = -1 0 0 0 2 -1 0 0 0 0 
0 -1 0 0 -1 2 0 0 0 0 
-1 0 0 0 0 0 2 -1 0 0 
0 -1 0 0 0 0 -1 2 0 0 
-1 0 0 0 0 0 0 2 --1 
0 -1 0 0 0 0 0 -1 2 


and a! — 2409 +243a8 — 136007 +4618a° — 97920° + 1277404 — 95200? + 3141a?+216a—297 = 
(a — 1)° (a — 3)° (a? — 8a + 11)(a? — 4a — 1) 


And the characteristic polynomial of Bs and Bg is given by (a — 1)* (a — 3)* (a? — 9a + 
13)(a? — 5a — 1), (a — 1)” (a — 3) (a? — 10a + 15)(a? — 6a — 1), respectively. 


Generally, the characteristic polynomial of B,, using the Laplacian domination matrix is 
(a —1)”"* (a— 3)” * (a? — (m+ 4)a + (2m + 3))(a? — ma — 1) =0. 


solving the equation we get 


(a —1)"7! =0, or (a — 3)! =0, or (a? —(m4+2)a+(2m+3)) = 0, or (a2—ma-—1) = 0. 
So a=1,1,1,---,1 ((m—1)times), or a = 3,3,3,--- ,3 ((m— 1)times), and (a? — (m+2)a+ 
(2m +3) = 0, 


a= 5(m 4— f/m? 28) and 
a2 = 5(m 44+ /m? 28) here m > 3, 
(a? — ma — 1) =0, 
1 
a3 = 5 (m— vim? +4) and 
1 
aa = 5 (m+ vine +4), 
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ELy—min = Ery—min(G) = ye |ax| 
i=1 


= (m= 1) 436m — 1) +|5 (2Vne FA) | + | (2G +4) 


Therefore, Ezy—-min = Exy—min(G) = 5m+ Vm? +4 4. This completes the proof. 


Theorem 3.4 For m > 3, the minimum Laplacian distance domination energy of a Book 


Graph (Bm) is 10m — 5 + 36m? — 48m + 49. 


Proof The characteristic polynomial of B,, for m > 3 can be found directly. 


For m = 3, B3 is a book graph with 8 vertices. The minimum dominating set is S = 
{vi,v2}. The Laplacian distance domination matrix and the characteristic polynomial of B3 
are respectively calculated by 


ow wv wo 


and 6* — 1847 + 296° + 16126° — 1662964 + 755366? — 1810326? + 2223366 — 110160 = 
(3 — 2)? (6 -— 6)? (6? — 98 +.17)(G? + 78 — 45). 

Similarly, calculation shows that the Laplacian distance domination matrix and the char- 
acteristic polynomial of By are respectively given by 


4 1 1 2 1 2 1 2 1 2 

-1l 4 2 i 2 i 2 1 2 1 

-1 -2 2 1 2 3 2 3 2 3 

—2 -1 -l 2 3 2 3 2 3 2 

LD,(G) = 1 2 2 3. 2 1 2 3 2 3 
2 1 3 2 1 2 3 2 3 2 

1 2 2 3 2 30°20 SoS 2 33 

2 1 3 2 3 2 1 2 -3 -2 

1 2 2 3 2 3 2 3.2 -i1 

2 1 3 2 3 2 3 2 ee 
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G'° — 246° + 5568 + 420887 — 661926° + 4942723° — 217865644 + 59343364° — 98012164? + 
89856008 — 3504384 = (9 — 2)° (9 — 6)° (62 — 116 + 26)(62 + 116 — 78). 

The characteristic polynomial of Bs is given by (6 — 2)! (B- 6)! (6? —136+37)(67+15G 
121), and the characteristic polynomial of Bg is given by (3 — 2)° (3 — 6)” (6? — 158 +50)(8? 4 
193 — 174). 

Generally, the characteristic polynomial of B,, using the Laplacian distance domination 
matrix is 


(@—2)™"* (8 —6)™"* (G? — (2m + 3)8 + (m+ 1) + 1)(6? + (4m — 5) 8 — (5m? — 2m + 6)). 


Solving the equation we get (@ — 2)" = 0, or (@ — 6)" = 0, or 6? — (2m +3)8+ (m4 


1)? +1 =0, or 6? + (4m — 5) — (5m? — 2m +6) = 0. So 6 = 2,2,2,--- ,2 ((m — 1)times), or 


% 9 


B = 6,6,6,--- ,6 ((m— 1)times), and 6? — (2m+3)8+(m+4+1)?+1=0, 


By = 5 (2m +3—V7V4m+ 1) and 
Bo = 5 (2m +34 V4m+ 1) here m 2 3, 


B? + (4m — 5)G — (5m? — 2m 4+ 6) = 0, 
1 


pes (—V36m? = Asin: 40 — dent 5) and 


3 
fa = ; (36m? 48m + 49 — 4m 4 5), 


EL py—min ia EL py—min (G) 


n 1 1 
= 5° |6i| = 8(m—-1) + plan +9) - E (2 5m? — in +49) | 
i=1 


Whence, Eppy—-min = E_py—min(G) = 10m — 5+ V36m? — 48m + 49. This completes the 


proof. 
Ws We W, 
Wo Wio Wi 


Figure 3 Wheel graph W,,4<n< 11 


W4 


Definition 3.2 A wheel graph W,, of order n, sometimes simply called an n-wheel, is a graph 
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that contains a cycle of ordern—1, and for which every graph vertex in the cycle is connected 
to one other graph vertex (which is known as the hub). The edges of a wheel which include the 
hub are called spokes. The wheel W,, can be defined as the graph K, + Cn_1, where Ky is the 


singleton graph and C,, is the cycle graph. Some wheel graphs are shown in Figure 3. 
Theorem 3.5 Forn > 4, the minimum dominating energy of a wheel graph (W,,) is > V/4n — 3. 


Proof We can find the characteristic polynomial of W,, for n > 4 by calculation directly. 


2) 


U3 U4 


Figure 4 W, 


For n = 4, W4 is a wheel graph with 4 vertices. The minimum dominating sets are S = {v1} 
or S = {v2} or S = {v3}. 

For S = {v,} the domination matrix and the characteristic polynomial of W4 are respec- 
tively calculated by 


A,(G) me 


es — ee me ey 


1 
0 
1 
1 


Oo FP Re 


1 
1 
0 
1 


and K+ — 6° — 6K? — 5K —1 = (Kk? —3K—1)(K2 +2441). The characteristic polynomial is found 


to be same when S = {v2} or S = {vs}. 


For n = 5, Ws is a wheel graph with 5 vertices. The minimum dominating sets is S = {v;}. 
Calculation shows that the domination matrix and the characteristic polynomial of Ws; are 
respectively given by 


a 

2 

— 

Q 

Nees 

II 
Le er ee oe 
So FPF FF CO Fe 
KE Oo Oo ee 
eS Sy Oe ce oe 
OS cee He oS oe 


and 4° — K* — 8x3 — 4K? = (K? — 3K — 2)(K3 + 2x7). 
Similarly the characteristic polynomial of Wg, Wz and Wg are given by (Kk? — 3K — 3)(K?2 + 
Kk —1)?, (Kk? — 3K —4)(K—1)?(K +1)?(K + 2) and (Kk? — 3K — 5)(K? + K? — 2k — k)”, respectively. 


Generally, the characteristic polynomial of W,, for n > 4 using domination matrix is 


[K? — 3K — (n— 3)] P(x). 
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Solving the equation (K? — 3« — (n — 3) = 0 we get m1 = $(3- V4n—3) and kp = 
2 
$ (3+ V4n—3). Ey-min = Ey—-min(G) > > |ki|, Ey-min(G) > /4n —3. This completes the 


proof. 


Theorem 3.6 For n > 4, the minimum distance dominating energy of a wheel graph (W,,) is 


> V4n? — 24n + 45. 


Proof The characteristic polynomial of W,, for n > 4 can be obtained by calculation 
directly. 


For n = 4, W4 is a wheel graph with 4 vertices. The minimum dominating sets are S = {v;} 
or S = {vo} or S = {v3}. For S = {v1}, Calculation shows that the distance domination matrix 
and the characteristic polynomial of W4 are respectively given by 


S 

2 

— 

Q 

~~" 

II 
ee Be 
pa en SS 


Fe CO fF Fe 
Oc. Ae ce 


and o* — 0° — 607 — 50 — 1 = (0? — 30 — 1)(0? + 20 +1). 
The characteristic polynomial is found to be same when S = {v2} or S = {vs}. 


For n = 5, Ws is a wheel graph with 5 vertices. The minimum dominating sets is S = {v,}. 
Calculation shows that the distance domination matrix and the characteristic polynomial of Ws 
are respectively given by 


=) 

oS 

Q 

Se 

II 
Le ee ee OE oe 
NO FF CO 
Fe Oo NY FY Fe 
oO FP FP Ye 


1 
1 
0 
2 
1 


and o° — o* — 1603 — 200? = o(0? — 50 + 0)(o 4+ 2)?. 


Similarly, the characteristic polynomial of We, W7 and Wg are given by (0? — 7¢ +1)(07 + 
30 + 1), 0(07 — 90 + 2)(0 + 1)?(0 + 3)? and (0? — 110 + 3)(o? + 50? + 60 + 1)?, respectively. 


Generally, the characteristic polynomial of W,, for n > 4 using distance domination matrix 


Solving the equation (0? — (2n — 5)a + (n— 5)) = 0 we get 


O1 = 


02 = 
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2 
and Epy—min = Epy—min(G) > > |oi|, Ep y—min(G) > V4n2 — 24n + 45. Hence, we complete 
i=l 


the proof. 


Theorem 3.7 For n > 4, the minimum Laplacian domination energy of a wheel graph (W,,) 


is > Vn? —2n+5. 


Proof Calculation enables one to find the characteristic polynomial of W,, for n > 4 
directly. 

For n = 4, W4 is a wheel graph with 4 vertices. The minimum dominating sets are S = {v1} 
or S = {vo} or S = {v3}. For S = {v1}, Calculation shows that the Laplacian domination 


matrix and the characteristic polynomial of W4 are respectively given by 


2 -1 -1 -l 
-1 3 -1 -1 
-1l -1 3 -1 
ae is CR 


and at — 1la? + 39a? — 40a — 16 = (a? — 3a — 1) (a — 4)’. 

The characteristic polynomial is found to be same when S = {v2} or S = {vs}. 

For n = 5, Ws is a wheel graph with 5 vertices. The minimum dominating sets is S = {v,}. 
The Laplacian domination matrix and the characteristic polynomial of Ws are respectively 


calculated by 


-l1 -1 0 3 -1 
-1 0 -1 -1 3 


and a® — 15a* + 8203 — 190a? + 141a + 45 = (a2 — 4a — 1) (a — 3)? (a — 5). 

Similarly, the characteristic polynomial of We, Wz and Wg are given by (a? — 5a — 
1) (a? — 7a +11)” , (a?—6a—1) (a — 2)? (a — 4)? (a—5) and (a?—7a—1) (a® — 10a? + 31a — 29)”, 
respectively. 

Generally, the characteristic polynomial of W,, for n > 4 using Laplacian domination 


matrix is 


[a? — (n—1)a—1] P(a). 


Solving the equation (a? — (n — 1)a— 1) = 0 we get 


1 
ar=5(n ‘i n? 2n +5) 


and 
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2 
ELy—min = Ery—min(G) > y |ax;| =V n2 — 2n +5. 
w=1 


Hence the proof is completed. 


Theorem 3.8 Forn > 4, the minimum Laplacian distance dominating energy of a wheel graph 
(W,,) is > V9n? — 62n + 117. 


Proof The characteristic polynomial of W,, for n > 4 can be obtained by calculation 
directly. 


For n = 4, W, is a wheel graph with 4 vertices. The minimum dominating sets are 
S = {ui} or S = {vg} or S = {ug}. For S = {v,}, Calculation shows that the Laplacian 
distance domination matrix and the characteristic polynomial of W4 are respectively given by 


2. hm al, pak 
-l1 3 -1 -1 
-1 -1 3 -1 
ede Sh. SL. 43 


and G4 — 1163 + 3962 — 406 — 16 = (62 — 36 — 1) (6 — 4)’. 
The characteristic polynomial is found to be same when S = {v2} or S = {vs}. 


For n = 5, Ws is a wheel graph with 5 vertices. The minimum dominating sets is S$ = 
{v1}. The Laplacian distance domination matrix and the characteristic polynomial of Ws are 
respectively given by 


ED AG) SCAN sai 3s Op St 


and 


BP 15b* 4: Tap? — 048? 9458 595 = (02 98 —7)(8 = 5) (8 3): 


Similarly, the characteristic polynomial of We, W7 and Wg are given respectively by 


(6? — 8-17) (62 -98 +19)’, 
(8? + 06 — 31) (6 — 6)* (6 — 4)" (6-3) 


and 
(6? + B — 49) (63 — 146? + 638 - 91)”. 


Generally, the characteristic polynomial of W, for n > 4 using Laplacian distance domi- 


nation matrix is 


[B? + (n — 7)B — (2n? — 12n + 17)} p(8). 
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and 
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Solving the equation 3? + (n — 7) — (2n? — 12n + 17) = 0 we get 


fi. = 5(-v9n®—6antii7—n+7), 
1 
B. = 5 (Von? =62n+117—n +7) 


2 
Ep py min = Ei pymin(G) > Da |6;| = V/9n2 — 62n + 117. 


i=l 


Hence the proof is completed. 
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Abstract: A function f is called a C-geometric mean labeling of a graph G(V, EF) if 
f : V(G) — {1,2,3,---,q@ + 1} is injective and the induced function f* : E(G) — 
{2,3,4,--+,q+1} defined by f*(uv) = | FOFO)| for all wv € E(G) is bijective. A 
graph that admits a C-geometric mean labeling is called a C-geometric mean graph. In this 


paper, we have discussed the C-geometric meanness of some ladder graphs. 


Key Words: Labeling, C-geometric mean labeling, C-Geometric mean graph, Smaran- 


dache k-mean graph. 
AMS(2010): 05C78. 


§1. Introduction 


Throughout this paper, by a graph we mean a finite, undirected and simple graph. Let G(V, EF) 
be a graph with p vertices and q edges. For notations and terminology, we follow [5]. For a 
detailed survey on graph labeling we refer to [4]. 

Path on n vertices is denoted by P,. G© S», is the graph obtained from G by attaching m 
pendant vertices at each vertex of G. Let G; and G2 be any two graphs with p; and p2 vertices 
respectively. Then the cartesian product G1 x G2 has pip2 vertices which are {(u,v) : u € 
Gi,u € Go}. The edges are obtained as follows: (ui, v1) and (ug, v2) are adjacent in G, x G2 
if either uy = ug and v; and v2 are adjacent in G2 or wu, and ug are adjacent in G; and 
UV, = vg. The ladder graph L, is a graph obtained from the cartesian product of P2 and Py. 
The triangular ladder TL,,n > 2 is a graph obtained by completing the ladder L,, by the edges 
ujvi41 for 1 <i <n—1, where L,, is the graph P, x P,. The slanting ladder SL, is a graph 
that consists of two copies of P,, having vertex set {uj:1<i<n}U{uj:1<i< n} and edge 
set is formed by adjoining uj+1 and v; for all 1 <i<n-—1 ([2]). 

Let P, be a path on n vertices denoted by u11, U1,2, U1,3,°°* ,Ui,n and with n— 1 edges 
denoted by €1,€2,:-+ ,€n—1 where e; is the edge joining the vertices u1,; and ui,41. On each 
edge e;, erect a ladder with n — (i — 1) steps including the edge e;, for i = 1,2,3,...,n—1. 
The graph thus obtained is called a one sided step graph and it is denoted by ST}. Let Pon be 
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a path on 2n vertices u1,1, U1,2, U1,3,°°* ,U1,2n and with 2n — 1 edges e1, €2,...,€2n—1 where e; 
is the edge joining the vertices u;,; and u1,;41. On each edge e;, we erect a ladder with ‘7 + 1’ 
steps including the edge e;, for i = 1,2,3,...,n and on each e; erect a ladder with 2n + 1-7 
steps including e;, for? =n+1,n+2,--- ,2n—1. The graph thus obtained is called a double 
sided step graph and it is denoted by 2S7),. 


The study of graceful graphs and graceful labeling methods was first introduced by Rosa 
[7]. The concept of mean labeling was first introduced by S. Somasundaram and R. Ponraj [8] 
and it was developed in [6] and [9]. In [11], R. Vasuki et al. discussed the mean labeling of 
cyclic snake and armed crown. In [1, 3], some graph labelings of step graphs were analyzed. 


In a study of traffic, the labeling problems in graph theory can be used by considering the 
crowd at every junctions as the weights of a vertex and expected average traffic in each street 
as the weight of the corresponding edge. If we assume the expected traffic at each street as the 
arithmetic mean of the weight of the end vertices, that eases mean labeling of the graph. When 
we consider a geometric mean instead of arithmetic mean in a large population of a city, the 
rate of growth of traffic in each street will be more accurated. Motivated by this, we establish 
the geometric mean labeling on graphs. 


Motivated by the works of so many authors in the area of graph labeling, we introduced 
a new type of labeling called C-geometric mean labeling. A function f is called a C-geometric 
mean labeling of a graph G if f : V(G) > {1,2,3,---,q+ 1} is injective and the induced 
function f* : E(G) — {2,3,4,---,q+1} defined as 


f*(uv) = | F(u) F)| for all uv € E(G) 


is bijective. A graph that admits a C-geometric mean labeling is called a C-geometric mean 
graph. 


In [10], S. Somasundaram et al. defined the geometric mean labeling as follows. 


A graph G = (V, E) with p vertices and q edges is said to be a geometric mean graph if it 
is possible to label the vertices « € V with distinct labels f(a) from 1,2,--- ,qg+1 in such way 
that when each edge e = wv is labeled with f(uv) = | F(a) F@)| or | F(u) Fe)| then the 
edge labels are distinct. 


In the above definition, the readers will get some confusion in finding the edge labels which 
edge is assigned by flooring function and which edge is assigned by ceiling function. Generally, 
a graph G = (V, F) with p vertices and q edges is said to be a Smarandache k-mean graph for 
an integer k > 2 if it is labeled vertices x € V with distinct labels f(x) from 1,2,---,q+ 1 in 
such way that edge e = wv is labeled with f(uv) = | MTOOFFOF | or | FOF | then the 
edge labels are distinct. Clearly, a Smarandache 2-mean graph is nothing else but a geometric 
mean labeling graph. 


n [10], S. Somasundaram et al. have given the geometric mean labeling of the graph 
Cs UC, as shown in Figure 1. 
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3 4 4 11 11 40 10 9 


Figure 1 A geometric mean labeling of C5 U C7. 


From the above figure, for the edge uv, they have used flooring function | f(iwf @)| and 
for the edge vw, they have used ceiling function [Vv fiwf @)| for fulfilling their requirement. 


To avoid the confusion of assigning the edge labels in their definition, we just consider the 


ceiling function [Vv f(u)f @)| for our discussion. Based on our definition, the C-geometric 
mean labeling of the same graph Cs U C7 is given in Figure 2. 


Figure 2 A C-geometric mean labeling of Cs U C7 


In this paper, we have discussed the C-geometric mean labeling of the ladder graphs L, 
forn > 2, Ln © Sm forn > 2 andm < 2, TL, forn > 2, TL, © Sm forn > 2 andm < 2, SL, 
for n > 2, SL; © Sp for n > 2 and m < 2, step graph ST), and double sided step graph 2ST ,. 


§2. Main Results 


Theorem 2.1 The graph Ly is a C-geometric mean graph for n > 2. 


Proof Let uy, u2,--+,Un and v1,V2,+-++ ,Un be the vertices of L, = P, x Pj. Then the 
ladder graph L,, having 2n vertices and 3n — 2 edges. 


Define f : V(Ln) > {1,2,3,...,3n— 1} as follows: 


f(u) = 1, 

f(u;) =3i-1, for2<i<n, 
f(v1) =3 and 

f(u) = 3i- 2, for2<i<n. 
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Then the induced edge labeling is obtained as follows: 


Hence, f is a C-geometric mean labeling of the ladder P,, x P2. Thus the ladder P,, x P 
is a C-geometric mean graph for n > 2. 


Theorem 2.2 The graph Ly, © Sm is a C-geometric mean graph for n> 2 andm < 2. 


Proof Let uy, u2,++* ,Un and v1, V2,°++ , Un be the vertices of Ly, = P, x P2. Let wl, wh, 


. we and al, ao, ee 


ladder L,,, for 1 <i<n. 


as) be the pendent vertices attached at each vertex u; and v,; of the 


Case l. m=1. 


Define f : V(Ln © S1) > {1,2,3,--- ,5n — 1} as follows: 


f(u) = 3, 

f(u;) =5i-—3, for2<i<n, 

f(v1) =4, 

flu) =5i-2, for2<i<n, 
f(w®) =5i-4, forl<i<n, 
f(a) = 2 and 

f(a) = 5i—1, for2<i<n. 


f* (uti) = 5t, forl <i<n-1, 
f* (wivigi) = 5i+1, forl<i<n-1, 
f*(ui1) = 4, 
f (ui) = 5¢-— 2, for2<i<n, 
f* (usw?) =5i-3, forl<i<n, 
f (ve) = 3 and 
fi (vie) =5i-1, for2<i<n. 


Case 2. m=2. 
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Define f : V(Ln © S2) > {1,2,3,--+ ,7n — 1} as follows: 


3 i=l 
flu)=4 71-2 2<i<nand i is even 
7i—5 2<i<nandiis odd, 

5 i=l 
f(vi) = Ti —4 2<i<n and iis even 
7-1 2<i<nandiis odd, 

1 i=l 
f(wl?) = Ti —3 2<i<n and? is even 
7i—6 2<i<nandiis odd, 


sa eS 
f(e®) = 7i-6 3<i<nand iis even 
7i—3 3 <i<7nand 7 is odd 
8 i=1 
and f(a) = Ti —5 2<i<n andi is even 
7i—2 2<i<nandiis odd. 


Then the induced edge labeling is obtained as follows: 


6 t=1 
Ti 2<i<n-l, 


f*(uitizi) = 


f* (wivigi) = Ti +1, forl<i<n-1, 
f* (uv) = 71-3, for l<i<n, 
2 t=1 
f*(ujw®) = Ti —2 2<i<nand iis even 


71 —5 2<i<nandiis odd, 


c (i) 7i—1 1<i<nand iis even 
f*(uiwy") = 

71-4 1<i<nand iis odd, 

(i) 71 —5 1<i<nand iis even 
f(vity’) = 


7i—2 1<i<nand iis odd 
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7 z=1 
and f (vic) =4 71-4 2<%i<nand i is even 


Ti-1 2<it<nandiisodd. 


Hence, f is a C-geometric mean labeling of the graph L,, © S,,. Thus the graph Ly, © Sm 


is a C-geometric mean graph for n > 2 and m < 2. 


Theorem 2.3 The graph TL, is a C-Geometric mean graph for n > 2. 


Proof Let the vertex set of TL, be {ui, U2,++* , Un; V1, V2,°** ;Un} and the edge set of TL, 
be {ujuizi; 1 <i< n—-1U {vviga; 1 <i<n—-1U ful <i < n}u {oii 1 <i<n-1}. 
Then TL, has 2n vertices and 4n — 3 edges. 


Define f : V(TLn) > {1, 2,3,--- ,4n — 2} as follows: 


f(ui) = 4¢ - 3, for 1 <i<n, 
f(u) =4¢-—1, for 1 <i<n-—J1and 
f(un) = 4n — 2. 


Then the induced edge labeling is obtained as follows: 


f  (uimigi) = 40-1, for 1 <i<n-1, 

f* (uivi) = 40-2, for 1 <i<n, 
f* (vivigi) = 41 +1, for 1 <i<n-1 and 
f* (viuigi) = 4¢, for 1 <i<n-1. 


Hence f is a C-geometric mean labeling of TZ,. Thus the triangular ladder TL, is a 


C-geometric mean graph for n > 2. 


Theorem 2.4 The graph TLy, © Sm is a C-geometric mean graph for n > 2 andm < 2. 


(4), (@) (2) 


Proof Let ui, t2,++* , Un and v1, v2,°-+ , Un be the vertices of TL,,. Let wi’, w5’,-++ , wm and 
v2)... ,2) be the pendent vertices attached at each vertex u; and v; of the ladder Ln, 
forl<i<n. 


Case 1. m=1. 
Define f : V(TLn © 51) > {1,2,3,--- ,6n — 2} as follows: 


f(u1) = 3, 


f(u;) =6i-—4, for2<i<n, 
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f(u) =6i-2, forl<i<n, 


6 
f(w) =6i—5, forl<i<n, 

2 

67 —3, for2<i<n. 
7-1, forl<i<n-1, 
t+1, forl<i<n-1, 


2, forl<i<n-1, 


i—3, for2<i<n 


’ 


Ss 
* 
— 
= 
I 
D> w D> Da F&F D@D ®D 


f (uw? =6i-—4, forl<i<n, 
f* (vy? = 3 and 
fi (vj0 = 61-2, for2<i<n. 


Case 2. m=2. 


Define f : V(TLn © S2) > {1,2,3,--- ,8n — 2} as follows: 


f(u1) = 3, 

f(ui) = 8-3, for2<i<n, 

f(v1) = 5, 

f(u;) = 8i—5, for2<i<n, 
fey) =t 
f(w\”) =8i —4, for2<i<n, 
f(w)) = 2, 
f(wS?) = 81-2, for2<i<n, 
fat?) =4, 

f(a) = 81-7, for2<i<n, 
f(a?) = 6 and 
f(c$?) = 8i —6, for2<i<n. 


Then the induced edge labeling is obtained as follows: 


f*(uiu2) = 7, 
f* (uimi41) = 8i+1, for2<i<n-1, 
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f* (wrv2) = 8, 

f* (wivigi) = 81-1, for2<i<n-1, 
f* (uivi) = 8t-— 4, for l<i<n, 
f*(viu2) = 9, 

f* (viii) = 8i, for2<icn-1, 

FP (urwy") = 2, 

f* (uw?) = 81-3, for2<i<n, 

f*(urwy”) = 3, 

f* (uws”) = 81-2, for2<i<n, 

F(z’) =5, 

f* (ujay’) = 8i —6, for2<i<n, 

f* (irs) = 6 and 
f* (uivs’) = 81-5, for2<i<n. 


Hence, f is a C-geometric mean labeling of the graph TL, © S,,. Thus the graph TL),© Sm 


is a C-geometric mean graph for n > 2 and m < 2. 


Theorem 2.5 The graph SL, is a C-geometric mean graph for n > 2. 


Proof Let the vertex set of SL, be {u1, W2,++* , Un, V1, V2,°°° 


, Un} and the edge set of SL, 


be {ujtiga;3 l <i<n—1SU {oyi41;1 <i <n —-—1}U {ujuigaj;1 <i<n-—1}. Then SL, has 


2n vertices and 3n — 3 edges. 


Define f : V(SL,) > {1,2,3,--- ,3n — 2} as follows: 


fw) = 1, 

f(uis) = 3t—4, for 2<i<n, 
f(uj) = 3t, for 1 <i <n-—1 and 
f(vn) = 3n—-2 


Then the induced edge labeling is obtained as follows: 


f* (uru2) = 2 
f* (uimigi) = 3t-— 2, for2<i<n-1, 
f* (wivigi) = 314+ 2, for l<i<n—-2, 
f* (Un-1Un) = 38n — 2 and 
f* (viuigi) = 3t, for 1 <i<n-1. 


Hence f is a C-geometric mean labeling of SL,,. Thus the slanting ladder SL, is a C- 


geometric mean graph for n > 2. 
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Theorem 2.6 The graph SLy,© Sm is a C-geometric mean graph for n > 2 andm < 2. 


Proof Let uz, u2,°++ ,Un and v1, V2,°++ , Un be the vertices of SL,,. Let w, wh, ee wand 
oe), ... ¢ be the pendent vertices attached at each vertex u; and v; of the ladder Ln, 


forl<i<n. 


Case 1. m=1andn>3. 


Define f : V(SLy, © S1) > {1,2,3,--- ,5n — 2} as follows: 


f(w?) =1, 

f(w?) =5i-7, for2<i<n, 

fay = 

fic) =bi+1, for2<i<n-—land 
5 


Then the induced edge labeling is obtained as follows: 


3t 1<i<2 
f*(uitizi) = 

5i- 3 3<i<n-Il, 
f* (wivigi) = 5i +3, for l<i<n-2, 
f*(n—10n) = 5n — 3, 
f* (viuigi) = 5i, for 1 <i<n-1, 
F*(uwl”?) =2, 
f*(ujw?) = 51-6, for2<i<n, 
fret?) =7, 
f*(vi2) = 5i+1, for2<i<n-—1and 
f* (vax ) =5n —2 


Case 2. m=2andn>3. 
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2i+1 
flu) =4 71-6 
71-9 
7n — 10 
f (Un) = 
7n—9 
9 
fw) =) Ti+2 
Ti-1 
7n— 13 
f (vn—2) = 
7n—15 
f(Un-1) = 7n— 5, 
f(vn) = 7n—- 3, 
1 
fw) =4 8 
7i—7 
7i— 10 
7n — 10 
41 —2 
fw) =. 5 
7i-8 
euga| oT 
7n—8 
8 
FGP)=4 7% 
7i-3 
(ee ya 
7n—17 
ae eae 
in—7 


Define f : V(SLy © S2) > {1,2,3,--+ , 7 — 2} as follows: 


1<i<2 
3<i<n-—1 and iis even 
3<i<n-—1andiisodd, 
n is even 

n is odd , 

i=1 
2<i<n-—83 andi is even 
2<i<n-—3 andi is odd, 
n is even 


n is odd , 


t=1 

2<1<3 
4<i<n-—1 andi is even 
4<i<n-—landiisodd, 
n is even 

n is odd , 

1<i<2 
3<i<n-—1 and iis even 
3<i<n-—1andiisodd, 
n is even 

n is odd , 

i=l 
2<i<n-—83 andi is even 
2<i<n-—3 andi isodd, 
nm is even 

n is odd , 

nm is even 


n is odd , 
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11 
fe?) =< 41 
7i—-2 
Ce ae 
7n — 16 
f(a") = 7-6 
and f(a”) =7n—2. 
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i=1 
2<i<n-—83 andi is even 
2<i<n-—83 andi is odd, 
n is even 


n is odd , 


Then the induced edge labeling is obtained as follows: 


di 
Putin) = 9 
7i-4 
7n —13 
a (Un—1Un) = 
7n—11 
: 12 
f° (viviei) = 
7i+4 
7n —9 
fe (Un—2Un—1) = 
7n — 10 


f* (Wn-1Un) = Tn — 4, 


nm is even 
n is odd , 
1=1 
2<i<n-3 


9 


nm is even 


n is odd , 


f* (viui4i) = 71, for 1 < a <n- 1 ; 


2 
‘ 671 — 7 
Paw) = 2 
71—6 
71-9 
7n — 10 
f* (up) = 2°” 
7n—9 
37 
fr(usw$?) =) 7-5 
71-8 
f*(unw”) = Tn — 8, 
9 
F(a?) =) 41 


i=1 

2<1<3 
4<i<n-—1 andi is even 
4<i<n-—1andiisodd, 
n is even 

n is odd , 

1<i<2 
3<t<n-—1 and? is even 


3<i<n-—1andiisodd, 


i=1 
2<i<n-—83 andi is even 


2<i<n-—3andiisodd, 
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7n— 12 n is even 


7n — 16 n is odd , 


f* (vary) = Tn —3, 
10 G1 
f (vic?) = Ti +2 2<i<n-—83 andi is even 
Ti-1 2<i<n-—3 andi is odd, 
r (vm 206") _ 7n— 11 n is even 
7n— 15 n is odd , 
i (vps?) =7n—5 
and f* (vnx$”) = 7n — 2. 


Case 3. m=1,2 andn=2. 


The C-geometric mean labeling of SL2z © S$; and SLz © So is given in Figure 3. 


Figure 3 The C-geometric mean labeling of SL2 © S; and SL2 © S3. 


Hence, f is a C-geometric mean labeling of the graph SL, © S,,. Thus the graph SL, © Sm 
is a C-geometric mean graph for n > 2 and m < 2. 


Theorem 2.7 The graph ST, is a C-geometric mean graph for n > 2. 
Proof Let 1,1, U1,2, U1,35°°* ; Ui,n; U2,1, U2,2, U2,3,°°* » U2,n, U3,1, U3,2, U3,3,°°* > U3,n—1; U4,1, 
U4,2, U4,3;°°* 5 U4n—2,°°* 5 Un,1,Un,2 be the vertices of the step graph ST). 


In u;,;,7 denotes the row (counted from bottom to top) and j denotes the column (counted 
from left to right) in which the vertex occurs. 
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Define f : V(ST,) — {1,2,3,--- ,n? +n—1} as follows: For 1 <i<n-1, 


(n+ 1-1)? +2(j-1) Legs || 
f (wig) = en ae - 
(n+1-—i)(n+3-—%)-274+1 [pet | +1 <j<n+1-i, 


f(tinto-1) = (n+1-i)(n+3—-i1), for2<i<n-1, 
f (un) = 38 and 
f (un,2) =); 


Then the induced edge labeling is obtained as follows: 


For 1<i<n-—- 2, 


3 
+ 
ee 
| 
~. 
Ny 
| 
T 
i) 
QQ. 
| 
eS 


Soa et 
j= | >= | and 7 is odd 
l 


3 
ran 


ate | and i is even 


( 
FP (wij 4ij41) = | 
( 


n+1—i)(n+3-i%)-2j nte=tl|+1<j<n-i, 


f* (Un—1,1Un—1,2) = 5 
f* (UnUn,2 =2 
(n+1-—i)(n+2-1%), for2<i<n-1, 


=(n+1-i)(n—-i), forl<i<n-2, 
=4 


) 

) 
ik (Uint1—iVitin+2—-) i 
ia (ui 1Uiti,1) 
) 


Te (Un—1,1Un,1 


For 1<i<n-8, 


f* (ui guia) = 


nee | and i is even 


(nt+1—-i)(n-i)+2j-1 2<j<|Se"/-1 

(n+ 1-)(n-—7 427-1 j = | 2" | and i is odd 
(n+1—4)( F=| 

(n+1—4)( 


nteil 41 <j<n-i, 
f* (Un—2,2Un—1,2 =o, 
ss (Un—1,2Un,2 


f (Uint1—iUi+inti—i =(n+1—- i)?, for l<i<n-2 


) 
) 
) 
and f*(Un—1,2Un,2) = 3. 


Hence, f is a C-geometric mean labeling of ST;,. Thus the step graph ST}, is a C-geometric 


mean graph, for n > 2. 


Theorem 2.8 The graph 2ST2,, is a C-geometric mean graph, for n > 2. 


Proof Let uy ,1, 01,2, U1,3,°°* , Un, 2,1, U2,2, U2,3,°°* , U2,2n, U3,1, U3,2, U3,35°°* » U3, 2n—2, Ud,1, 


U4,25 U4,3,°°* 5 U4, 2n—4,°** » Un41,1,Un+41,2 be the vertices of the double sided step graph 2ST. 
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In u;,;,7 denotes the row (counted from bottom to top) and j denotes the column (counted 


from left to right) in which the vertex occurs. 


Define f : V(2STon) > {1,2,3,--+ ,2n? + 3n} as follows: 


2n?+n+1+2(7-1) 1<j<n 


fur) = . 
2n? + 38n—-2(j—n—-1) n+1<j<2n, 


for2<i<nand2<j<n+2-1i, 


f(usj) =2n+1—-1)? + (n+2—-%) +29 — 2), 


for2<i<nandn+3-—i<j7 < 2n+3- 21, 


f(uij) = 2(n + 1-1)? +3(n+1-%) - 206+ 5 -—n—-83), 
f(u21) = 2n? +n —-2, 
f(u1,1) = 3, 
f(ui2) = 1, 
f(ui1) =2(n+2-i)P +n—-i, for3<i<nand 
f (Uizn+4—21) = 2(n + 2 i)? n+1-1, for2<i<n. 


Then the induced edge labeling is obtained as follows: 


2n?+n4+2+2(7-1) 1<j<n 
Qn? + 3n+1-2j—n) nt+1<j<2n-1, 


f* (ui ju1,341) = 


for2<i<n—land2<j<n+2-i, 


f* (wigs gi) = 2(n +11)? + (n+3—-%) +2(j - 2), 


for2<i<n—-—landn+3-i< 7 < 2n4+2-—- 21, 


(n+1—1)? +3(n4+1-—2)4+1-20+3—n—2), 
n+2-—i), for2<i<n-1, 


f (u11U2 j= 4r +n, 
f* (U1,2nU2,2n) = n+n+1, 
f (Ui 2Uida1 = (n 1 i)? tn+1—i, for2<i<n-l, 


f* (Un 2Un41,1 
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2n? —n+2+2(j — 2) 2<j<n 


f* (ur,ju2,5) = 
- n2+n—-1-Aj—-n—-1) n+1<j<2n-1, 


for2<i<n—land3<j<n+2-i, 


FP (ij viti,j—1) = 2(n+1 i)? (n t 1 i) t 2(7 2), 


for2<i<n—-—landn+3—-i< 7 < 2n4+2-—- 21, 


is a C-geometric mean graph, for n > 2. 


f* (ui juitij—1) = 2(n +1- i)? Tr (n +4—- i) _ 2(i +7 Te 1), 
f*(us1us2) = 2(n + 1-1)? +3(n4+1—i4) +1, for2<i<nand 
= 2(n+1-i)? +3(n+1-—i)+2, for2<i<n. 


f* (Ui2n4+3—21Ui,2n44—21) 


Hence, f is a C-geometric mean labeling of 2S'T>,,. Thus the double sided step graph 2S7>,, 
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Abstract: The maximum order of partition of the edge set E(G) into edge hub sets is 
called edge hubtic number of G' and denoted by €-(G). In this paper, we determine the edge 
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§1. Introduction 


By a graph G = (V,F), we mean a finite and undirected graph without loops and multiple 
edges. A graph G with p vertices and q edges is called a (p,q) graph, the number p is referred 
to as the order of a graph G and q is referred to as the size of a graph G. In general, the degree 
of a vertex v in a graph G denoted by deg(v) is the number of edges of G incident with v. The 
degree of an edge uv is defined to be deg(w) + deg(v) — 2. Also A‘(G) denotes the maximum 
degree among the edges of G, and 6‘(G) denotes the minimum degree among the edges of G. 
|a| is the greatest integer less than or equal to x. In a tree, a leaf is a vertex of degree one, a 
leaf edge is an edge incident to a leaf. We refer to [6] for terminology and notations not defined 
here. 

Introduced by Walsh [13], a hub set in a graph G is a set H of vertices in G such that 
any two vertices outside H are connected by a path whose internal vertices lie in H. The hub 
number of G’, denoted by h(G), is the minimum size of a hub set in G. A connected hub set in 
G is a vertex hub set F' such that the subgraph of G induced by F' (denoted G[F]) is connected. 

Let G be a graph, let e = (u,v) and f = (u1,v1), a path between two edges e and f 
is a path between one end vertex from e and another end vertex from f such that d(e, f) = 
min{d(u, ur), (u, v1), (v, ur), (v,v1)}. Internal edges of a path between two edges e and f are 
all the edges of the path except e and f [11]. A subset H, C E(G) is called an edge hub set 
of G if every pair of edges e, f € E\ H. are connected by a path where all internal edges are 
from H.. The minimum cardinality of an edge hub set is called edge hub number of G, and 
is denoted by h.(G) [11]. An edge hub set H. C E(G) is called a connected edge hub set, if 
the subgraph [H-] is connected. The minimum cardinality of a connected edge hub set of G 
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is called a connected edge hub number and is denoted by hee(G) [1]. For more details on the 
hub studies we refer to [10]. Graphs G1, and G2 have disjoint vertex sets V; and V2 and edge 
sets E, and E> respectively. Their union, G = G; U G2 has, as expected, V = V; UV2 and 
E = E, UE, [6]. 


A set D of vertices in a graph G is called dominating set of G if every vertex in V \ D is 
adjacent to some vertex in D, the minimum cardinality of a dominating set in G is called the 


domination number 7(G) of a graph G ([7]. 


A set B of edges in a graph G is called an edge dominating set of G if every edge in E\ B 
is adjacent to some edge in B, the minimum cardinality of an edge dominating set in G is called 
the edge domination number 7/(G) of a graph G ([7]). An edge-domatic partition of G is a 
partition of E(G), all of whose classes are edge-dominating sets in G. The maximum number of 
classes of an edge-domatic partition of G is called the edge-domatic number of G and denoted 
by ed(G) ([1]). 


A double star S,,m is the tree obtained from two disjoint stars Ky,,-1 and Ky,m—1 by 
connecting their centers [5]. The line graph L(G) of G has the edges of G as it is vertices which 
are adjacent in L(G) if and only if the corresponding edges are adjacent in G [6]. A friendship 
graph, is the graph obtained by taking m copies of the cycle graph C3 with a vertex in common 
and denoted by F;,. The following results will be useful in the proof of our results. 


Theorem 1.1([10]) For any graph G, he(G) < q— A’(G), and the inequality is sharp for any 
path Py, p> 4. 


Proposition 1.1({10]) For any graph G, he(G) < p— 3. 


Theorem 1.2([10]) For any tree T with p > 3 vertices and I leaves, 
A(T’) = hee(T) =p — (0+ 1). 


Proposition 1.2({9]) For any graph G, €(G) < 6(G) + 2. 


§2. Main Results 


Definition 2.1 The maximum order of partition of the edge set E(G) into edge hub sets is 
called edge hubtic number of G and denoted by €-(G). The maximum order of partition of the 
edge set E(G) into connected edge hub sets is called connected edge hubtic number of G and 
denoted by Ece(G). 


It is obvious that €.(G) > &.(G), since he(G) < hee(G). We first determine the edge hubtic 
number of some standard graphs. 
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Observation 2.1 (1) For any cycle Cp, 


3, ifp=3; 
4, ifp=4; 
Ee(Cp) = ; 
2, ifp=5,6; 
1 


, ifp> 7. 


(2) For any path P,, 
3, ifp=4; 
£e(Py) = 2, if p= 3,5; 
, ifp>6. 


(3) For the wheel graph Wj p-1, p > 4, 


6, ifp=4; 
&e(Wi»-1) = 4 4, ifp=5; 
3, ifp>6. 


(4) For the star IY p-1, €e (B51 p-1) =p— 1. 

(5) For the double star Sim, €e(Snjm) = 3. 

(6) For the complete bipartite graph Knim, $e(Knm) = max{n, m}. 

We will check that if the edge hubtic number is a suitable measure of stability?. Now we 
ask, does the edge hubtic number discriminate between graphs. There are many examples of 


graphs which propose that €.(G) is a suitable measure of stability which is able to discriminate 
between graphs. For example, consider the graphs G , Gz and G3 in Figure 1. 


YY 


Gy Go G3 
Figure 1: G , Go, and G3. 


It is clear from Figure 1, that ed(G1) = ed(G2) = ed(G3) = 3, the edge domatic number 
does not discriminate between graphs G1, Gz and G3, but €.(G1) = 3, €-(G2) = 4 and €.(G3) = 
5, therefore €.(G1) # (G2) # (G3). So the edge hubtic number discriminates between 
graphs G, Gz and G3. 
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Observation 2.2 For any graph G, 0 < &(G) < q. 


Theorem 2.1 If a graph G is a tree with at least 3 non-leaf edges and the induced sub graph 
G|(E \ L)] is not a star where L is the set of all leaf edges in G, then €.(G) = 1. 


Proof Let a graph G be a tree with at least 3 non-leaf edges and the induced sub graph 
G|(E \ L)] is not a star, we discuss the following cases: 


Case 1. Suppose that H, is a set of all non-leaf edges, clearly any path between two leaf edges 
does not pass through another leaf edge. So, H, is an edge hub set of G, and by Theorem 1.2 
it is minimum edge hub set. Now, suppose Z C EF \ H. be an edge hub set of G. Since G is a 
tree with at least 3 non-leaf edges and the induced sub graph G[(E \ L)] is not a star, then the 
induced subgraph G[E \ Z,] is not complete. Also any path in a tree never passes through a 
leaf edge. Therefore there are at least two non adjacent edges e, f € FE \ Z. such that no path 
between them is in Z,, this is a contradiction. Hence H, is the only edge hub set. 


Case 2. Suppose that H, is an edge hub set of G but not containing all non-leaf edges. Since G 
has at least three non-leaf edges, let {e1, e2, e3} be non-leaf edges where e; and e3 not adjacent, 
let 11,13 be two leaf edges adjacent to e; and es, respectively. Clearly, G[{l1, e1, e2, 3, l3}] is a 
path Ps. As h.(Ps) = 3, then H, contains at least three edges from Ps. Therefore any other 
edge hub set of G must intersects H, since size of Pg is 5. Then €.(G) = 1. 


Proposition 2.1 For any (p,q)-graph G, &(G) < my? where h.(G) £0. 


Proof Let H = {H,, Ho, H3,--- , Hz}, be the edge hubtic partition of G and €.(G) = t. 
Clearly |H;| > he(G),i = 1,2,--- ,t and we get ¢ = S~\_, |H;| > th-(G), hence the result. 


Observation 2.4 Let G’ be a subgraph of G, then is not necessary €_(G’) < €-(G). 
For example, G = Ki + Pu, and G’ = Ki + Ps, &(G’) =5 £3 =€(G). 
Proposition 2.2 For any (p,q)-graph G of order p > 5, 


&e(G) < 0 (G) +2. 


Proof By the definition of edge hub number it is obvious that h.(G) = h(L(G)), so 
&-(G) = €(L(G)). By Proposition 1.2, €.(G) = €(L(G)) < 6(L(G)) + 2, since 6’(G) = 6(L(G)), 
the result follows. 


Corollary 2.1 For any (p,q)-graph G of order p> 5, 


&e(G) +he(G) < 5(G)+p-1. 


Proof By Proposition 1.1 and Proposition 2.2, we get the result. 
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Theorem 2.2 For any (p,q)-graph G of order p, €¢(G) +&-(G) < be) , and the inequality is 
sharp for stars K13, and Ky,. 


Proof By Observation 2.2, €.(G) <q and &( 


G 
A 
aI 
FS 
ap 
oO 
=) 


€e(G) + &(G) < q+ 


OI 
II 


Theorem 2.3 Let G be a (p,q)-graph. Then 


Ee(G) +he(G) <q4+2. 


Proof By Theorem 1.1, he(G) < q— A’(G). Hence he(G) < q— 6'(G). Proposition 2.2, 
completes the proof. 


Observation 2.5 If € (G1) = &(G2),then not necessary he(G1) = he (G2). 


For example, G; = K1,3, and Gz = F such that €.(Gi) = (G2) = 3, and h.(G1) =0 F 
3 = h, (G2). 


Theorem 2.4 Let G be a graph of size q. Then €.(G) = q if and only if G with 0’ > q—2. 


Proof Assume that €.(G) = q, then there is a q partition of E(G) into edge hub sets and 
every partite set consists of one edge, we have the following cases: 


Case 1. All edges of G are adjacent, so any edge of G is an edge hub set of G. So 6’ =q-1. 


Case 2. Any edge of degree q — 1, is adjacent to all edges and hence it constitute an edge 
hub set of G, and since any edge of degree q — 2, is adjacent to all edges of G except one, so 
every edge of them must be an edge hub set for G, hence 6'(G) = q — 2, if we consider any 
edge f such that deg(f) < q—2, in this case let deg(f) = q —3, so there is two edges e1, e2 not 
adjacent to f, now if the set {f} is an edge hub set for G then e; must be adjacent to e2, but 
by this assumption {e,} is not edge hub set for G, since eg not adjacent to f and e; not a path 


between them. So €.(G) = q only if the graph G satisfies 6’(G) > q — 2. Converse is obvious. 


Proposition 2.3 For any two connected graphs Gy and G2, 


1, of Gy or Go is with 6’ <q-1; 


Ee(G1 UG2) = 
2, if Gy and G2 are with 6’ =q—-1. 


Proof Let G1, Gg be two graphs both with 6’ = q—1, clearly E(G1) is an edge hub set for 
G1 UG» and E(G2) is an edge hub set of the same graph, therefore €.(G1 U G2) = 2. Suppose 
that G, or G2 is with 6’ < q—1, then any edge hub set of Gj UG»2 must contain all of the edges 
of G,; and any edge hub set of Gg, therefore €.(G1 U G2) = 1. 


Corollary 2.2 For any disconnected graph G with m > 3 components, &.(G) = 1. 
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Abstract: There are 3 crises in the development of mathematics from its internal, and 
particularly, the 3th crisis extensively made it to be consistency in logic, which finally led 
to its more and more abstract, but getting away the reality of things. It should be noted 
that the original intention of mathematics is servicing other sciences to hold on the reality 
of things but today’s mathematics is no longer adequate for the needs of other sciences such 
as those of theoretical physics, complex system and network, cytology, biology and economy 
developments change rapidly as the time enters the 21st century. Whence, a new crisis 
appears in front of mathematicians, i.e., how to keep up mathematics with the developments 
of other sciences? I call it the 4th crisis of mathematics from the external, i.e., the original 


intention of mathematics because it is the main topic of human beings. 


Key Words: Mathematical crisis, reality, contradiction, TAO TEH KING, mathematical 
universe hypothesis, Smarandachely denied axiom, Smarandache multispace, mathematical 


combinatorics, traditional Chinese medicine. 
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§1. Introduction 


As we known, one or the main function of mathematics in science is it can establish exact math- 
ematical expressions for scientific models on things. Certainly, a theory can not be without the 
practice, and it can be only from the practice. By this view, the creating source of mathematics 
can be only from solving problems appeared in practice of human beings, and then move its 
method and technique upward a mathematica theory for understanding the reality of things in 

the world. 
Usually, a thing is complex, even hybrid with other things sometimes. Then, what is the 
reality of a thing? The reality of a thing is its state of existed, exists, or will exist in the world, 
1 An extended version of the preface of my book Mathematical Reality - My Philosophy on Mathematics with 


Reality, Published in USA, August, 2018. 
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independent on the understanding of human beings, which implies that the reality holds on by 
human beings maybe local or gradual, not the reality of a thing. Hence, to hold on the reality 
of things is the main objective of science in the history of human development. 

But, a mathematical conclusion really reflects the reality of a thing? The answer is not cer- 
tain because the practice of human beings shows the mathematical conclusion do not correspond 
to the reality of a thing sometimes, for instance the Ames Room. Usually, the understanding 
of a thing is by observation of human beings, which is dependent on the observable model, 
data collection by scientific instruments with data processing by mathematics. Such an ob- 
servation brings about a unilateral, or an incomplete knowledge on a thing. In this case, the 
mathematical conclusion reflects partial datum, not all the collection, and in fact, all collection 
data (by different observers with different model) with data processed is not a mathematical 
system, even with contradictions in usual unless a data set, which implies that there are no 
mathematical subfields applicable. 

We all know that it appeared 3 crises in mathematical development. In each time, mathe- 
matics itself was enriched, improved and completed. However, along with the solving process of 
the 3th mathematical crisis, the trend of mathematical developing in 19th and 20th centuries 
shows that a mathematical system is more concise, and its conclusion is more extended, then 
farther to the reality of things because it abandons more and more characters of things. Besides, 
more and more researchers only pay attention on questions or problems in himself branch along 
with the mathematical branch divided, and few peoples consider his research whether is or not 
valuable for developing the whole mathematics, for understanding the nature and beneficial to 
human progress, which finally results in mathematics father to the practice of human beings, 
weaker for hold on the true face of things in the world. 

As the time enters the 21st century, science developments change rapidly, and meanwhile, 
a few global questions constantly emerge, such as those of local war, food safety, epidemic 
spreading network, environmental protection, multilateral trade dispute, more and more ques- 
tions accompanied with the overdevelopment and applying the internet, ---, etc. Clearly, 
today’s mathematics is no longer adequate for the needs of other sciences. It is far falling be- 
hind the development of society. A new crisis appears in front of mathematicians, i.e., how to 
keep up mathematics with the developments of other sciences for hold on the reality of things? 
I think this is a big and more important problem in the development of mathematics in the 
21th century, and call it the mathematical 4th crisis because holding on the reality of things is 
the central objective of human beings. 

The main purpose of the review is analyzing this crisis and points out the way of one how 
to out this crisis by establishing new mathematical theory, also provides an envelope theory, 


i.e., mathematical combinatorics as the candidate for the way. 


§2. Be Understood or Not 


For reality of things, an elementary but fundamental question should be answered first. That 
is, can one really holds on the reality of things? For this question, there are two but quite 


opposite answers. One is the reality of things can not completely understanding, i.e., one can 
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only holds on the approximate reality of things. Another is one can finally understanding the 
reality of all things, i.e., Theory of Everything. We respectively discuss them following. 


Not Understood. There is a well-known philosophical book: TAO TEH KING written 
by an ideologist Lao Zi in ancient China. In this book, it discussed extensively on the relation 
of TAO, a more general object than the reality with name and things, and shown in its first 
but central chapter ([8]): 


The Tao that experienced is not the eternal Tao; 

The Name named is not the eternal Name; 

The unnamable is the eternally real and naming is the origin of all particular things; 

Freely desire, you realize the mystery but caught in desire, you see only the manifestations; 

The mystery and manifestations arise from the same source called darkness; 

The darkness within darkness, the gateway to all understanding. 

For explaining the relation of Tao with knowing ability of human beings respectively in his 
Chapter 42: 

Tao gives birth to One, One gives birth to Two, Two gives birth to Three and Three gives birth to all things; 

All things have their backs to the female and stand facing the male. When male and female combine, all things 


achieve harmony. 


and also in Chapter 25: 
Human beings follow the earth, Earth follows the universe, 


The universe follows the Tao and the Tao follows only itself. 


By the view of Lao Zi, the reality of things is not understood because the Tao that expe- 
rienced is not the eternal Tao, the Name named is not the eternal Name, and human beings 
is born after Three along with Three gives birth to all things, particularly, the reality. I agree 
Lao Zi’s notion, i.e., it is difficult to know the reality of all things, and all mathematical reality 
is only approximate reality, not the reality. For Tao, One and Two before Three, we can only 
analyze their various possibility by science, can not really hold on their true faces. 


Be Understood. The notion is the supporting and main trending in scientific community 
today, i.e., the reality of all things can be understood by human beings and one can finally holds 
on and become the dominate of the world. Particularly, the physical world is nothing else but 
a mathematical structure ([12], [13]) by Tegmark Max, a famous Swedish-American physicist 
and cosmologist in MIT now. 

Here, I would like to analyze 2 hypotheses, i.e., the Big Bang and mathematical universe 
hypothesis on the physical world. 


1. Big Bang Hypothesis. The Big Bang model states that the earliest state of the 
Universe was an extremely hot and dense one, and that the Universe subsequently expanded 
and cooled, which is based on general relativity following: 


Applying his principle of general relativity, i.e. all the laws of physics take the same form 
in any reference system and the equivalence principle, i.e., there are no difference for physical 
effects of the inertial force and the gravitation in a field small enough, Einstein got the equation 
of gravitational field 


1 
Ruv — 5 hou + Agu = —87GT yp, 
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where Ruy = Rup = Ro 


piv? 
Ore. 4or? 
R® = be pee ie ae T2.T? = To T? ; 
piv Ox” Oxi pir av pv ai? 
rs -_ 1 pa(O9mv “a Ognp _ 29mn 


Ou” Ou™ OuP 
and R= g’" Ryy. 


Combining the Einstein’s equation of gravitational field with the cosmological principle, 
i.e., there are no difference at different points and different orientations at a point of a universe 
on the metric 104l.y., Friedmann got a standard model of universe. The metrics of the standard 


universe are 


dr? : 
ds? — — dt? + ela + r? (dé? + sin? Ody")| 
and 2( 
R(t 2 p2 2239 
Git = 1, Orr = [oR 908 = —r*R*(t)sin* 0. 


The standard model of cosmos finally enables the birth of Big Bang model for the cosmos 
in thirties of the 20th century, and finally, the NASA’s explorer mission WMAP(Wilkinson 
Microwave Anisotropy Probe) determined the radius of the universe was 13.7).l.y on Big Bang 
hypothesis. 


Mathematical Universe Hypothesis. The mathematical universe hypothesis proposed 
by Max Tegmark, is a speculative Theory of Everything, which claims that our external physical 
reality is a mathematical structure([12], [13]), i-e., the physical universe is not merely described 
by mathematics, but is mathematics (specifically, a mathematical structure), which implies 
the mathematical existence equals to that the physical existence, and all structures that exist 
mathematically exist physically as well. And observers, including humans ourself, are self- 
aware substructures (SASs), and in any mathematical structure complex enough to contain 
such a substructure, it will subjectively perceive itself as existing in a physically real’ world. 


4 All Mathematical Structure 


Non-Visible Universe Non-Visible Universe 


Big Bang Universe 


(Visible Universe) 


Non-Visible Universe Non- Visible Universe 
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According to Lao Zi’s birth ruler, the WMAP is essentially determined the radius of visible 
universe by human beings is 13.7b.l.y, but we can not claim the radius of universe is 13.7b.l.y 
just by the Big Bang hypothesis. Otherwise, we are in an awkward situation and can not answer 
what is the outer of the sphere of radius 13.70./.y unless its radius is finite. The advantage of 
Max Tegmark’s hypothesis is it avoids the finite or not of the universe but claims its physical 
reality is a mathematical structure. These 2 hypotheses are simply shown in Fig.1. 

Therefore, the Big Bang hypothesis is only a notion locally on the universe. But why various 
experimental of human beings verifies it maybe right just because our human beings are after 
Three, i.e., after the Big Bang by Lao Zi, and the Friedmann’s standard model of universe is a 
special solution of Einstein’s gravitational equations, which is essentially to explain the general 
by special cases. However, we have many solutions on Einstein’s gravitational equations, even 
with constant \ = 0 ((2]). Certainly, the Max Tegmark’s hypothesis is on the whole universe 
but it also contains lethal deficiency following: 


(1) If the Big Bang hypothesis is right, i.e., we can only hold on the reality of the visible 
universe, how can we verify the external universe, i.e., non-visible unverse is mathematics or 
not; 

(2) Is our mathematical theory can already be used for understanding the reality of all 
things in the world? The answer is not certain because mathematics is homogenous without 
contradictions, i.e. a compatible one in logic but contradictions exist everywhere in the world 
by philosophy. Thus, the reality known by mathematics on things can be only a subset of the 
reality set ([4], [5]), ie., the mathematical structure is not equal to the physical reality. 


All of these show that even if the Big Bang and Max Tegmark’s hypotheses are both right, 
we also need to establish a new mathematical theory so that the mathematical structure is 


equal to the physical reality, i.e., a mathematical crisis is confronted with mathematicians. 


§3. Mathematical Crisis in 21th Century 


3.1 Brief Review 3 Crises of Mathematics in History 


As we known, there are 3 crises in the development of mathematics following, each of them 


motivates mathematics itself constantly enriched, improved and completed.. 


First Crisis. The early Pythagorean mathematics was based on the so-called commen- 
surability principle, i.e., Pythagorean’s assertion: “everything is a number”. According to this 
principle two geometric values Q and V have common measure, divisible by it, i.e., their ratio 
can be expressed as the ratio of the relative prime numbers m and n. However, Hippasu, a 
member of Pythagorean’s found the length of the diagonal of a unit square is 2, which can 
not be as a ratio of two relative prime numbers, i.e., it is an irrational number. This discov- 
ery became a turning point in mathematics development, which ruined the former system of 
Pythagoreans, extended the rational to real numbers and finally resulted in new mathematical 
theories. 


Second Crisis. Even at present, calculus is a subject with the most widely applying 
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to other science for hold on reality of things. However, its foundations refers to the rigorous 
development of the subject in its early time. The cause was the unrigorous use of infinitesimal 
quantities in that time, which resulted in the second crisis of mathematics, i.e., the foundation of 
calculus. Certainly, there are many mathematicians work hard for going out this crisis, formed 
new mathematical theories. For example, the limitation of Weierstrass eventually became 
common of calculus base, instead of infinitesimal quantities as the rigorous approach, and 
established real analysis which included full definitions, theorems with rigorous proof of calculus. 


Third Crisis. The third crisis of mathematics came from the foundation of mathematics, 
i.e., set theory by Russell paradox following: 


Let R be the set of all sets that are not members of themselves. If R is not a member 
of itself, then its definition dictates that it must contain itself, and if it contains itself, then it 
contradicts its own definition as the set of all sets that are not members of themselves, i.e., 


R={ale¢gar}, thn RERSRER. 


Russell paradox finally resulted in the establishing of axiomatic set theory, i.e., Russell’s 
type theory and the Zermelo set theory in 1908. 


3.2 Mathematical Crisis in 21th Century 


The mathematical crisis, or the 4th crisis in 21th century does not come from its internal but in 
the external needs or in its original intension. As we discussed, the axiomatic and abstract on 
mathematics in the 19th and 20th centuries finally results in mathematics away from practice. 
This trend also found by physicists in 20th century. Einstein once complain mathematics: as far 
as the laws of mathematics refer to reality, they are not certain; and as far as they are certain, 
they do not refer to reality. Besides, more and more problems appeared in the practice can not 
find an applicable mathematics and don’t know how to hold on their characters. In fact, there 
are more examples supporting this claim with social development in 21th century. 


Elementary Particle. We have known matters consist of two classes particles, i.e., 
bosons with integer spin n, fermions with fractional spin n/2, n = 1(mod2), and by a widely 
held view, the elementary particles consists of quarks, leptons with interaction quanta including 
photons and other particles of mediated interactions ([6], [7]), which constitute hadrons, i.e., 
mesons, baryons and their antiparticles. 


Fig.2 
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Although quark model is a formal classifying scheme for hadrons, i.e., the quarks and 
antiquarks of Sakata, or Gell-Mann and Ne’eman, it appeared subconscious in the multiverse 
interpretation of H.Everett on the superposition of particle. It should be noted that the mul- 
tiverse interpretation or quark is proposed by physicist for explaining behavior of particles 
without an applicable mathematics. However, it completely changed the usual notion that a 
particle is an geometrical point or a subset of space, and opened a new way for understanding 
the reality of a hadron in notion, i.e., we are not need to insist again that a hadron is a geo- 
metrical point or a subset of space such as those of assumptions in determinable science. For 
example, a baryon is predominantly formed by three quarks, and a meson is mainly composed 
of a quark and an antiquark in quark models, such as those shown in the right of Fig.2, where 
a particle consists of 5 quarks can be also found on the left. 


Biological Population. The biological populations are dependent each other by food 
web, i.e., a natural interconnection of food chains and a graphical representation of what- 
eats-what in an ecological community on the earth. For example, a food chain starts from 
producer organisms (such as grass or trees which use radiation from the sun to make their 
food) and end at apex predator species (like grizzly bears or killer whales), detritivores (like 
earthworms or woodlice), or decomposer species (such as fungi or bacteria). Usually, a model 
of a biological system is converted into a system of equations. The solution of the equations, by 
either analytical or numerical means, describes how the biological system behaves either over 
time or at equilibrium. In fact, a food web is an interaction system in physics which can be 
mathematically characterized by the strength of what action on what. For a biological 2-system, 
let xz, y be the two species with the action strength F’(a — y), F(y — x) of x to y and y to z 
on their growth rate, ([1]). Such a biological 2-system can be quantitatively characterized by 
differential equations 

t= Fly 2) 
y=F'(a>y) 


on the populations of species x and y. However, this method can be only applied to the small 
number (< 3) of populations in this web. If the number m of populations> 4, such as those 
shown birds in Fig.3, 


Fig.3 
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a natural way for characterizing the behavior of m birds is to collect all dynamic equations of 


cells, i.e., 
4 = Fy(t,%1) + ) A; (@j > 71) 
Al 
v2 = F(t, F2) + > Hi; (; —_— Fo) 
jP2 


where F; : R° — R? is generally a nonlinear function characterizing the external appearance of 
ith cell and H;(x; — 2;) is the action strength of the jth cell to the ith cell in this system for 
integers 1 <i,7 <™m. 

However, this system maybe non-solvable, i.e., one can not characterize the behavior of 
birds in the hope that of solutions. Thus, even in mathematical biology one has no a math- 
ematical branch applicable for the reality of biology unless partially by differential equations 
and statistical analysis. 


Regional Economy. Today, the regional trade regulation enables each one of his member 
develops extensively on other members. Achieving mutual benefit, and finally striving for trade 
balance is the purpose of the regional trade organization. This situation appears both in the 
global or area economy because there are few countries or areas still in self-sufficient today. 
The trade surplus and deficit usually result in the trade disputing in members, processes the 
multilateral negotiations, and then reachs a new regulation for members in the international 
trade. Usually, one can obtains statistical data published by customs or statistical services in 
a country or an area, but there are no a mathematical subfield for characterizing the global or 


local changing in economy. 


Furthermore, we can easily get other fields that there are no mathematical subfields appli- 
cable. For example, the complex network, including community network, epidemic spreading 


network with their behaviors. 


Why I say it to be the 4th crisis of mathematics? The reason is that mathematics is 
gradually ignored and replaced by computer simulation in society development. It already 
become a subject of indulge in self-admiration by mathematician ourselves. If we do not turn 
it around, it will be abandoned sooner or later by the society finally. 


§4. Out of the Crisis 


The essence of out of the 4th crisis of mathematics is return to the original intension, i.e., 
unveiling the reality of things. For such a objective, two main things should be processed 
first on today’s mathematics. One is the contradiction between things such that different 
things should be in equal rights, and another is the dependence of things because there is a 
universal connection between things by philosophy, which implies one should redefine elements 
in mathematics. 


Let us discuss the two things flowing. 
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1. Smarandache Multispace. Today, we have known a kind of geometry breaking 
through the non-contradiction in classical mathematics, i.e., Smarandache geometry (1969) by 
introducing a new type axiom for space. An axiom is said to be Smarandachely denied if the 
axiom behaves in at least two different ways within the same space, i.e., validated and invalided, 
or only invalided but in multiple distinct ways. A Smarandache geometry [10] is a geometry 
which has at least one Smarandachely denied axiom (1969). If & is a Smarandache denied 
axiom on space 7, then all points in Y with & validated or invalided consist of points sets 
TH” and TN™), and if it is in multiple distinct ways invalided, without loss of generality, 
let s be its multiplicity. Then all points of F are classified into T/”,T3”,--- ,T”. Hence, we 
get a partition on points of space J as follow: 


eT Pree: on. ea Tl Jes Pe 


This shows that .7 should be a Smarandace multispace. 


Generally, let (51;R1i), (H2;Re), +++, (Um; Rm) be m mathematical spaces, different two 
by two. A Smarandache multispace ¥ is a union U %; with rules R= U R; on ¥, ie., the 


w=1 i=1 
rule R; on b; for integers 1 < i < m ([3],[9]-[10]). Thus, the reality of things, whatever its 


accurate or approximate should be characterized or found out on Smarandache multispaces. 
Whence, the Smarandache multispace solved better the contradiction in classical mathematics. 
However, an abstract Smarandache multispace is nothing else but an algebraic or set problem 
({11]), which worked out finely the equal rights, but 


(1) To be also new conceptions accumulation; 
(2) Not solve the universal connection of things; 
(3) Can not extensively applies achievements in today’s mathematics, ---, etc.. 


Thus, for understanding the reality of things, a new envelope theory should be established 


on Smarandache multispace, i.e., mathematical combinatorics. 


U1 fi v2 
iz fa 
fi fs 
Un Vi+1 Us V3 
Fig.4 


2. Mathematical Combinatorics. What is mathematical combinatorics? The math- 
ematical combinatorics is such a mathematics over topological graphs G; i.e., establish an 
envelope mathematics on the elements of universal connection of things, which worked out 
finely both the equal rights and universal connection of things. And how to combine classical 
mathematics with topological graphs G? I found a typical set of labeled graphs Gt , called 
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continuity flows can be viewed as mathematical elements, i.e., labeling their edges by elements 
in a Banach space 4 with two end-operators on & and holding on the continuity equation 
=> => 
on each vertex in G. For example, such a continuity flow over C’, is shown in Fig.4, where, 
Aton, = 1 Aty,_, = 2 and . 
fi + (21 - 1) F(t,z) 


fi gi-1 


for integers 1 <i <n Then, such a set of labeled graphs G_* inherits the character of today’s 
mathematics, i.e., if ene fees vee mer are oriented topological graphs and & a Banach space, 
then all such labeled graphs Gt with linear end-operators is also a Banach space, and further- 
more, if @ is a Hilbert space, all such labeled graphs GE with linear end-operators is a Hilbert 
space too. 

Now, there are 2 kinds of problems on continuity flows GL: 

(1) Globally, given a graph family eee cneee aon > 1 and a Banach space &, 
whether there exists continuity flows over graphs GE, GE, a GL to be elements form a math- 
ematical space; 

(2) Locally, for a continuity flow GL, if some vertices are no longer conserved by outside 
interference, how to make it conserved again such that it is still a continuity flow. 


The first problem has been solved by a series papers of mine (See references of [5] in 
details), but for the second problems, there are only a few local or partially results. In fact, 
an independent energy system, including computer, car and human body, cell tissue, biological 
populations, ---, etc. adaptive system is nothing else but a continuity flow, and furthermore, 
conservation flow. Thus, we can use continuity flows to characterize behavior of these systems 
for reality. 

Here, I would like to introduce the twelve meridians theory in traditional Chinese medicine 
({14]), which can be viewed as a typical example of continuity flows, particularly, in treating 
an illness. It is in fact to make the patient balance in Yin and Yang on acupoints of meridians, 
i.e., conservation, where Yin (Y~) or Yang (YT) can be viewed as negative or positive energy, 
tendency, ---, etc. are basic conceptions in traditional Chinese culture, i.c., YT and Y~ are 


everywhere with that Yt in Y~ and Y~ in Y™, such as those shown in Fig.5, where the black 
and white areas respectively represent Y~ and YT. 


Fig.5 


According to the characteristics of human body, the traditional Chinese medicine proposed 
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12 meridian theory, i.e., there 12 meridians in human body completely reflects the physical 
condition. They are respectively Hand Tai Yang small intestine meridian (H,), Hand Shao 
Yang Tri-Jiao meridian (Hz), Hand Yang Ming large intestine meridian (H3), Hand Tai Yin 
lungs meridian (H4), Hand Shao Yin heart meridian (Hs), Hand Jue Yin pericardium meridian 
(Hg), Foot Yang Ming stomach meridian (F,), Foot Jue Yin liver meridian (F2), Foot Tai 
Yin spleen meridian (F3), Foot Shao Yin kidney meridian (F4), Foot Shao Yang gallbladder 
meridian (Fs), Foot Tai Yang bladder meridian (F¢), such as those shown in Fig.6(these red 
lines in human bodies without acupoint). 


Fig.6 


The balance of {Y~, Y*} at points on the 12 meridians is the basic ruler for human body 
in traditional Chinese medicine. If there exists a point in one of the 12 meridians in which 
{Y—,Y*} is imbalance, this person must be ill, and in turn, for a patient there are must be 
points on the 12 meridians in which {Y~,Y*} are imbalance. This is the healing theory of 
traditional Chinese medicine, and by thousands of years of testing, there are no counterexamples 
appeared in China. 

Certainly, the healing theory of traditional Chinese medicine is nothing else but continuity 
flows. Notice that the 12 meridians are in fact 12 directed pathes H, H2, H3, H4, H5, He, Fi, Fo, 
F3, Fy, Fs, Fg with vertices of acupoints. Define 


6 6 
#-(U»)U(Us) 
i=1 i=1 
with L: V(G) — {Y-,Y7}, then, G* should be conserved on its vertices in {¥", YT} fora 
person, i.e., a continuity flow. 

For a patient, i.e., there are points to be imbalance on the 12 meridians, the doctor detects 
the points on which meridians, at which acupoints and the imbalance is Y~ more than Y*, or 
Y* more than Y~, and then by a natural ruler of the universe in traditional Chinese culture, 
i.e., reducing the excess with supply the insufficient, the doctor regulates these related acupoints 
by acupuncture or drugs so that the acupoints balance in {Y~ , Y} again. Clearly, this implies 
a mathematical process for a continuity flow Gt again. 
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Certainly, there are no specific amount for the action strength H(a; — xo), where zo is the 
acupoint with {Y—,Y*} imbalance, x; is the related acupoints, 1 < i < s, which completely 
depends on the judgement of the doctor, and continuous regulation based on the actual situation 
of the patient, i.e., a process of response. This also implies that getting a continuity flow Gt 
again maybe by repeatedly regulation of the flows on conditions. 
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we know nothing of what will happen in future, but by the analogy of past 
experience. 


By Abraham Lincoln, an American president. 
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